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ABSTRACT 


Level Mountain is a Late Cenozoic shield volcano 
located near 131°20'W,58°31'N in the Stikine Volcanic Belt 
of northern British Columbia. The volcanic plateau with an 
average elevation of 4500',(1372m), is younger than Upper 
Miocene in age and is comprised of up to four sequences of 
alkali basalt and ankaramite flows and tuffs. The central 
region of peaks and ridges rises to a maximum elevation of 
7200',(2195m), and has a repetitive bimodal distribution of 
alkali basalt and peralkaline salic lavas and tuffs that 
spans the period from 4.5 million years b.p. to recent 
times. Feeble basaltic vents with spatter, bombs and scoria 
apparently postdate continental glaciation. A paleomagnetic 
study on two stratigraphic sections that span the entire 
range of volcanism samples the earth's major polarity 
reversals for the past 6 MY. 

Petrochemically the lavas belong to the sodic alkali 
basalt series and are of continental affinity. Oxygen 
isotope values near 5.6°/,,.0'*O SMOW have been measured for 
unaltered rocks including basalts, salics, and ultramafic 
inclusions. This would seem to indicate that the mantle 
under northwestern British Columbia is normal with respect 
to oxygen and that the Level Mountain magmas are mantle 
derived. Corroborating evidence is available from the whole 
rock strontium isotope studies which show a total range of 


*7Sr/**Sr from .7025 to .7071. The presence of basic and 
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Salic, undersSaturated and ovesaturated, peralkaline and 
metaluminous rocks attests to the complexity of petrogenetic 
processes at Level Mountain. Although the two volumetrically 
important lava types appear to be primary and mantle 
derived, the presence of some rhyolites enriched in ‘0, 
granitic gneiss incluSions in basalts and tristanites and 
great variation in whole rock lead isotope ratios are 
probably all indicative of some degree of interaction 
between magmas and crustal rocks. Major and trace element 
geochemical variations and mineralogy indicate an upper 
mantle origin at shallower than 15 kbar pressure from an 
undersaturated ultramafic source that possessed an abundance 
of alkalis and incompatible elements. The mineralogy and 
chemistry of the major rock types indicate conditions of low 
oxygen fugacity and a dry gas phase composition for both 
peralkaline flows and basalts. The field relations of some 
comendite flows implies fluid behavior and seems to indicate 
a viscosity two to three orders of magnitude lower than 
would be predicted on the basis of silica content and whole 
rock chemistry. This could be partially reconciled if the 
gas phase had been halogen rich and dry. 

Calculations of the volume and energetics of the 
volcanism of the Stikine Belt indicate that less than 15% of 
the modern heat flow is necessary to maintain the levels of 


volcanicity since Pliocene times. 
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CHAPTER 1. INTRODUCTION 


LEVEL MOUNTAIN VOLCANIC CENTRE 

Level Mountain is a major, alkaline and peralkaline 
continental shield volcano in the Intermontane Belt of 
British Columbia (Souther, 1977a). The description and study 
of the Late Cenozoic volcanics is of considerable importance 
to the further understanding of regional tectonics and 
volcanism in Western Canada. Investigations into the 
tectonics of this region have been published by Souther 
(1970, 1977b), Monger and Souther (1972), Gabrielse and 
Souther (1962), White (1959), and St. Armand (1957). 
Contemporary and recent studies on the volcanism of the 
Intermontane Belt include the work of Casey (1980), Bevier 
(1977, 1978), Fiesinger and Nicholls (1977), and Souther 
hone. oft v3 

The petrogenesis of alkaline volcanoes of bimodal 
character has been an open question and a topic for 
investigation by petrologists for many years, e.g. Daly 
(1933), Bowen (1938), Shand (1927), Chayes (1963, 1964), 
Gass and Mallick (1968), Gibson (1974), Bryan (1964, 1966, 
1976), and Clague (1978). The controversy over the origin of 
peralkaline rocks is a long standing one whose beginnings 
stem from investigations of the volcanics of the Isle of 
Pantelleria Washington (1913, 1914). The conditions for the 
appearance of peralkaline ferromagnesian minerals have been 


studied experimentally by Ernst (1962), Marsh (1975), and 
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Edgar and Parker (1974) with many questions still left to be 
answered. The mineralogy and petrochemistry of natural and 
synthetic peralkaline systems has been treated by Nicholls 
and Carmichael (1969), Bailey (1974) and Bailey and 


MacDonald (1970, 1975). 


OBJECTIVES 

The present study was begun in the Fall of 1974. The 
intent of the study was to geologically map and sample the 
Level Mountain Volcanics and to relate them through 
petrological and geochemical studies to Mount Edziza and 
Late Cenozoic volcanism of the Stikine Volcanic Belt as a 
whole. In addition to the original aspects of the study, 
which included mapping, petrography, petrochemistry, and 
mineralogy, the problem of accurately defining the age and 
Span of volcanism led to a paleomagnetic study and to 
attempts at radiometric dating. The questions on 
petrogenesis led to the acquisition of trace element and 
rare earth data in addition to analyses for oxygen, 
Strontium and lead isotopes. Other questions regarding the 
petrogenesis of the peralkaline salic lavas prompted a 
Series of one-atmosphere melting and crystallization 
experiments on natural compositions under conditions of 
controlled oxygen fugacity. Finally, attempts to relate the 
volcanism at Level Mountain to the regional picture of Late 
Cenozoic volcanism in Western Canada led to: a review of the 


existing literature on the geophysics and tectonics of the 


Intermontane Belt, field trips to other areas of the 
Intermontane, and some physical calculations regarding the 


nature of the volcanism. 


LOCATION, ACCESS, PHYSIOGRAPHY, AND CLIMATE 

The Level Mountain Range, centered near 
131°20'W-58°31'N, is the prominent physiographic feature of 
the Nahlin Plateau which is located in the Cassiar land 
district of northwestern British Columbia. The division 
between the Cassiar and Atlin mining districts crosses the 
plateau. The nearest road approach is to Telegraph Creek 
from the Stewart CassSiar Highway. Access to the area is best 
provided by helicopter from Dease Lake, 45 miles, (72 km), 
to the East. , 

The plateau is bounded by steep cliffs of columnar 
jointed basalt flows on the southern, western, and 
northeastern sides. The slopes to the northwest and 
southeast fall away more gently due to the influence of 
continental glaciation. The drainage from Level Mountain has 
a basically radial pattern, following the U-shaped glacial 
valleys that dissect the central region. The principal 
drainages are: Beatty Creek, Little Tahltan River, Egnell 
Creek, Dudidontu River, Kakuchuya Creek, Matsatu Creek, 
Megatushon Creek, Little Tuya River and Mansfield Creek. The 
southern and eastern watershed is part of the Stikine 
drainage system. The northern and western watershed is to 


the Nahlin River and the southwestern watershed is to the 
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Sheslay River. 

The climate of the region is rather severe due to the 
elevation and almost constant winds of greater then 15mph 
(24kph). Precipitation in the general area, (Atlin, 
Telegraph Creek, Dease Lake), is 10-15 in. (25-38 cm) per 
year. Summer weather is usally brisk, less than 13°C, and 
partly cloudy to overcast. Breakup does not usually happen 
until June, and during this time the plateau becomes a mess 
of deranged drainage, bogs and fens for up to 3 weeks. Snow 
fields usually remain throughout the summer on north facing 
slopes and may not ablate or melt for several years. The 
rather severe conditions on the plateau are conducive to the 
formation of desert varnish, which is seen on most outcrops 
and joint faces. The general geomorphology of the plateau is 
characteristic of low precipitation subalpine areas. Soils 
on Level Mountain are thin and generally show poor 
development and limited lateral extent. The two dominant 
types are acid well drained regosols on the steep slopes and 
brunisolic grey luvisols on the plateau surface and in 
alpine valley bottoms. Cryoturbation generally obscures any 
horizon development (Fengor, 1981) 

Tertiary Geomorphology 

The Tertiary Geomorphology of the Stikine region was 
characterized by mature drainage systems and gentle 
topography with a maximum of 1000m of relief; more typical 
variation was on the order of 150m. Most of the major 


lowland valleys such as the Tahltan River and the Dudidontu 
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River are either exumed or resequent valleys. This is in 
marked contrast to the young consequent streams of the lava 
plateau. The major river systems of the area such as the 
Nahlin-Inklin, the Shesley and the Stikine are resequent but 
rejuvenated by Late Tertiary uplift of the entire region. 
Investigations by Johnson (1926) and Kerr (1936) indicate a 
net uplift on the order of 1000m since Mid-Tertiary time. 
Even more extensive uplift has been reported for the 
adjacent coast range region (Parrish, 1981). This uplift has 
not only affected the drainage and geomorphology but is 
probably also related to the major tectonic processes 
responsible for the widespread alkaline volcanism and normal 
faulting of Late Tertiary times. There is a contrast in the 
geomorphology of the Intermontane Belt from north to south 
that could be explained by different rates of uplift. In the 
Stikine region the Intermontane Plateau has been incised by 
youthful rivers. Tertiary valleys that had recently held 
Pleistocene glacial deposits have been exumed. In contrast 
to this, the Tertiary topography of the Cariboo region is 
Still buried beneath glacial deposits and youthful lava 
plateaux. This is evident near Anahim Lake (Tipper 1969). In 
the Cariboo, the headward erosion of major rivers has not 
yet crossed the western boundary of the Intermontane Belt. 
The geomorphology of the southern Intermontane Belt has been 
more strongly influenced by the paths of Quaternary 
glaciation and subsequent streams than any of the area to 


the north or west (Fulton, 1975). 
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PREVIOUS WORK 

There are few previous geological investigations in the 
Level Mountain Range. Single traverses over or around the 
plateau margin were made by Cockfield (1925) and Mathews 
(Watson and Mathews 1944). During "Project Stikine" and the 
preliminary mapping of the Dease Lake Sheet (Gabrielse and 
Souther, 1962) a long traverse was made by Root and Ostensoe 
in the summer of 1958 (Ostensoe, 1960). In these early 
investigations of the Level Mountain Volcanics, field names 
for rock types include basalts as well as andesite, 
rhyolite, and quartz latite giving a misleading impression 
toward calc-alkaline character. Ostensoe (1960) was the 
first to penetrate the central region, collect rock samples 
and perform petrographic analyses. His preliminary findings, 
including the report of titaniferous augite in the basalts 
and aenigmatite from a comendite demonstrate the alkaline 
nature of the Level Mountain Lavas. 

Souther, of the Geological Survey of Canada, has done 
considerable regional mapping of the Late Tertiary and 
Quaternary volcanics of the Intermontane Belt. His studies 
have been highlighted by the work on Mount Edziza, the other 
major centre of the Stikine Volcanic Belt (Souther, 1966). 
His studies include detailed petrological (Souther, 1974), 
mineralogical (Yagi and Souther, 1974), and paleomagnetic 
work (Souther and Symons, 1974). His pioneering 
investigations into the volcanism and tectonics (Souther, 


1977a, 1977b) have revealed the presence of peralkaline 
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lavas and the possibility that the Northern Stikine Volcanic 
Belt may represent a situation of incipient continental 
EittanorASoutherny 1977b); 

The type locality of the Level Mountain formation 
(Souther,1971) is a continental shield volcano of uppermost 
Tertiary and Quaternary age centered near 131°20'W, 58°31'N. 
Alkalic lavas outcrop over an area of about 700 square miles 
(1812 square kilometers) in the form of a basaltic plateau 
with average elevation of 4500 feet (1372m). Capping the 
plateau is a geologically complex central region of 
glacially dissected ridges and vents with elevations to 7200 
feet (2195m). The distal portions of the shield consist of 
flat lying flows, with primary dips less than 2°, of alkali 
basalt, ankaramite and tuffs. Four mappable stratigraphic 
units have been recognized on the plateau (Hamilton and 
Scarfe, 1977). The central region with its dissected 
stratocone and domes has a bimodal distribution, with 
peralkaline trachytes and comendites predominating over 


basalts. 
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Chapter 2. REGIONAL SETTING 


HISTORICAL GEOLOGY OF THE INTERMONTANE BELT 

Reviews of the historical geology of the Cordillera 
have been given by White (1959) and Monger et al (1972). The 
focus of those studies was the geological evolution of 
British Columbia since Paleozoic time. Souther (1970) 
divided B.C. into five major northwest trending geological 
evolutionary belts. From west to east these are: The Insular 
Belt, The Coast Crystalline Belt, The Intermontane Belt, The 
Omineca Belt, and the Eastern Marginal Belt. These 
physiographic/tectonic regions and the distribution of 
Mid-Tertiary and older volcanic rocks is given in Figure 
2-1. The tectonic development of British Columbia and 
portions of Alaska has been discussed by St. Armand (1957). 
The Intermontane Belt (IMB) lies between 180 and 420 miles 
inland and parallel to the present day Pacific Coast. 
Intermittent but extensive volcanic activity since Paleozoic 
time has characterized the IMB. This distinction as a 
volcanic region continues on through Late Cenozoic and to 
Recent times (Souther, 1977a). 

The IBM is characterized by subdued topography and a 
multitude of small isolated Tertiary basins, each with its 
own sequence of continental sediments, lavas, and a veneer 
of Quaternary glacial deposits. The basement is 
characteristically comprised of Mesozoic calc alkaline 


plutonic rocks as major batholiths, and structurally 
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Figure 2-1, 
Physiographic belts of the Cordillera and the distribution 
of older volcanic rocks, adapted from Souther, 1970. 
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deformed marine sediments of Upper Paleozoic and Mesozoic 
age. Regional metamorphic rocks are rare and usually 
associated with the distinct margins of the IMB. 

The volcanism of the IMB from Paleozoic to Middle 
Cenozoic times has been dominated by calc-alkaline 
intermediate magmas. The sedimentary section from Permian 
time onwards has been strongly influenced by local sources 
and basin development has been dominated by local tectonic 
movements (Monger et al, 1972). These two observations tend 
to indicate that for most of the past 250MY the IMB has been 
located near a major plate boundary, in the modern tectonic 
sense, or at least the IMB has been coincident with an 
active continental margin. Since the onset of the Mesozoic 
orogenies with their widespread plutonism, the western part 
of British Columbia has been locally emergent (Souther, 
1967). Most of the IMB has been emergent since the close of 
Jurassic time. Volcanism of the Late Cretaceous appears to 
have been relatively extensive but erosion has taken its 
toll and even the major volcanic formations such as Ootsa 
Lake or Sloko are now represented by remnants of small areal 
extent that have been preserved in down-faulted regions 
(Souther, 1970). The volcanism of the Eocene and Oligocene 
was Subaerial and calc-alkaline with dominant control by 
local eruptive centres and shallow crustal magma chambers 
with compositions ranging from basaltic andesite through 
andesite, dacite, and rhyolite. This period of volcanic 


activity was definitively more extensive in the central and 
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southern portions of the IMB than in the north. A hiatus in 
volcanic activity occurred between Oligocene and Middle 
Miocene time. This period has been interpreted as 
corresponding to the cessation of subduction along the 
western Canadian continental margin with the disappearance 
of the Kula plate (DeLong et al, 1978). This was the time of 
most extensive Tertiary basin development. The resumption of 
volcanism in Miocene time brought a pronounced change in the 
chemical nature of the lavas. From Miocene to Recent time 
the volcanism has been predominantly alkaline and 
characterized by both major lava plateaux and isolated 
centres. Repeated glaciation of both continental and alpine 
type has occurred since Upper Miocene time. 

As one of the regions of most varied igneous and 
tectonic processes since the Paleozoic, the IMB is an 
important place to study the complex geological history of 
the continental margin and the evolution of the western 
North American crust and upper mantle. The traditional 
geosynclinal models, with their strong reliance on localized 
Or regional vertical crustal movements, are being replaced 
by a patchwork tectonics of Cordilleran microplates (Yorath 
and Chase, 1981; Riddihough, 1977, 1981). It is currently 
being proposed that many of the Cordilleran terranes, 
including the Stikine, are in fact allochthonous and derived 
from more southerly latitudes, with docking times ranging 
from the Mesozoic to the Mid-Tertiary (Monger and Irving, 


1980). With the possibility of a separate geological 
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evolution for adjacent areas of the Cordillera the value of 
a single historical geology and tectonic development is cast 


in doubt. 


LEVEL MOUNTAIN AND LATE CENOZOIC VOLCANISM IN BRITISH 
COLUMBIA 

The distribution of Miocene to Recent volcanism of the 
Intermontane Belt is presented in Figure 2-2. The Middle to 
Upper Miocene onset of this latest volcanic episode is 
preceeded by a hiatus of 10 to 15MY duration, and a 
fundamental change in the character of the lavas from calc 
alkaline to continental alkaline. Regional mapping by the 
G.S.C. (Operation Stikine) and reconnaissance by the author 
on the following sheets: Kettle River, Princeton, Nicola, 
Bonaparte Lake, Quesnel, Anahim Lake, Prince George, and 
Stikine River, demonstrate the two stage nature of Tertiary 
volcanism. This generalization holds from the Okanagan to 
the Cariboo and from the Omineca to the Stikine. The 
alkaline nature of the Miocene to Recent lavas of the 
Intermontane Belt is supported by whole rock chemical 
analyses presented in this study, as well as by independent 
conclusions by Souther (1977a) and Fiesinger and Nicholls 
(1977b). Souther (1977) has subdivided the Late-Cenozoic 
alkaline lavas of the Intermontane Belt into three belts: 
(i) Plateau Belt, central B.C., (ii) Anahim Volcanic Belt, 
western B.C., (iii) Stikine Volcanic Belt, northwestern B.C. 


Other volcanic belts (calc-alkaline) such as the Garibaldi 
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Figure 2-2. 
Miocene to Recent volcanism of the Intermontane Belt, 
adapted from Souther, 1970, 1977. 
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Belt and Wrangell Volcanic Belt are demonstrably related to 
active or fossil convergent plate boundaries and slabs of 
subducted crustal lithosphere (Green, 1977; Armstrong et al, 
1977). It has been Souther's contention that the alkaline 
lavas and their distribution are similarly associated with 
major lithospheric tectonic features. For the Anahim 
volcanic belt this could be a tensional region associated 
with the tear fault edge of the subducted slab which was 
responsible for the Garibaldi Volcanic Belt (Bevier et al, 
1979). For the Stikine Belt the feature might be related to 
a subducted ridge or incipient continental rifting (Souther, 
1977b). Typically the alkaline volcanics of the world have 
one of two principle settings. The most abundant are 
permitted, discordant types associated with lithosphere in 
tension along ridge-spreading centres - rifts and transform 
faults. The other type is tectonically uncorrelated such as 
on oceanic islands, and isolated alkaline piles that are set 
in tectonically complex areas (Sorensen, 1970; MacDonald, 
1974: Bailey, 1977). Unlike the second cycle calc-alkaline 
magmas that in some fashion involve remelting of crust 
(McBirney, 1969), the alkaline magmas are primary and 
reflect only the physical conditions and chemistry of the 
upper mantle (Sorensen, 1970; Wyllie, 1971). It is possible 
that the Stikine and Anahim Volcanic Belts are related to 
crustal tension associated with the shear geometry of the 
Pacific continental margin and the Queen Charlotte Fault 


system (Souther, 1977b). Another possibility is that the 


of betale> ¢Mastaiomst ots fisgipiasiter th 
to sdele Sra aeitabhnces S7a6G sthineuade’ San 
te 39 phosdemia {7TR! .eeT7e! eisteodtti Cateeae 

eee ne ientL0e nsed aad at 
apfsaud in 2a'h aise 


— P ~ 
cd 


efiitsila si> Fees | 


diiw betel ove2e 1£%6 j “° 6 w 
Pape 
a .@e sat oifdtnes stredees ‘4 
bejul>oses :« isactans? 2 6d Bigoo efae tied 
ze i daie bs - sit te sobs 2663 oated 
L@ Fs Siveao? J.is 3 ible i:1ad 209 125 
93 Beiala ‘ia and dé PSe one 
radauc of issnsi ,> peigiver 2e°-se96is om 
aa 
eved SLrow ts ebios olay «hi tamhh smo Garena 


#18 tacknwe sop 21? . epoziges ehgianiag oe 

gi erence rdttl fous 6) wee easy | +h$B rewetb a ) 
euolarand Bis egii5° - S22 So; ¢ TSS prgoa-sphia ene : 
26 AQwe= OS75.5 ohh)? L casi to dors $7 Ad 7) 2 (onze ot 
jea ot6 Jenmszt 26.4 - feds 3 CoJ Glee: eG . 2hwadei 
,5la0eGs6m ‘G.e wenstod), 2a97%6 *224q008 inate > 
etifatis-olso s.o¢5 Dhotet Sra sXiind PPE 
tauwio to gaizfomes evisveal sponded one 5 


baa “Pr eMisG. 414 shea =n Laws eds, 


ed 16. y2ade limes Ses <i le? 


eidiavoq at st ofh¥O) reer SOT a 


15 


Late Cenozoic volcanism of the Intermontane Belt is 
primarily related to thermal and compositional aspects of 
the upper mantle and only influenced to a lesser degree by 
variations in crustal fabric, thickness and composition. 

The alkaline lavas of the IMB are predominantly fluid 
basalts that fill local Tertiary basins to create extensive 
plateaux and local shield edifices. Examples of the basaltic 
plateaux are the upper Kamloops Group (Bonaparte Sheet), the 
Quesnel Group (Campbell and Tipper, 1971) and the Stikine 
Group (Souther, 1970, 1977a). Examples of more localized 
shields and composite volcanoes, which include peralkaline 
trachytes and other salic lavas, are the Itchas, Ilgachuz, 
and Rainbow Mountain Ranges of the Anahim Volcanic Belt and 
Edziza and Level Mountain of the Stikine Volcanic Belt. 

The distribution of lavas in the Stikine Volcanic Belt 
is given in Figure 2-3 which has been adapted from the 
1":4mile geologic map of Project Stikine (map 1957b) and 
Subsequent mapping by: Souther and Symons (1974), Hamilton 
and Scarfe (1977) and Casey (1980). Other than the spatial 
juxtaposition of the Edziza Volcanism with the Mess Creek 
Fault and the Level Mountain Volcanism with the Nahlin 
Fault, there is no apparent relation of Late Cenozoic 
Volcanism to Pre-Tertiary crustal structures. The wide 
disparity in age between the alkaline volcanics and the last 
activity on these old crustal structures (Gabrielse and 
Wheeler, 1961; Wheeler and Gabrielse, 1972; Souther and 


Armstrong, 1966) may demonstrate their spatial association 
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Figure 2-3. 
Geological sketch map of the Stikine Volcanic Belt showing 
the distribution of Late Cenozoic alkaline lavas and ancient 
faults. Adapted in part from: Stikine River Map 9-1959, 
Souther, 1970; Casey, 1980. LEGEND:(a) - Kawdy Plateau, (b) 
- Heart Peaks, (c) - Nahlin Fault, (d) - Level Mountain, (e) 
- Tanzilla Plateau, (f) - Tahltan Plateau, (g) - Stikine 
River, (h) - Mount Edziza, (i) - Castle Rock, (j) - Klappan 
Range, (k) - Mess Creek Fault, (1) - Hoodoo Mountain, (m) - 
Iskut River. 
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to be of no genetic significance. However, the Nahlin Fault 
may be a pervasive crustal structure (suture zone), formed 
as a result of the Mesozoic emplacement of the allochthonous 
Stikine terrane (Monger and Price, 1979). In that alkaline 
magmas are typically associated with extensional tectonics 
and fracture zones (MacDonald, 1974), the lavas of the 
northern Stikine Volcanic Belt (Heart Peaks, Kawdy Plateau, 
Level Mountain Range and Tanzilla Plateau) could be related 
in part to a Late Cenozoic reactivation of this old crustal 


weakness. 


CHAPTER 3. GEOPHYSICS OF THE INTERMONTANE BELT AND THE 


TECTONIC SETTING OF THE LEVEL MOUNTAIN RANGE 


GEOPHYSICS AND REGIONAL TECTONICS OF THE INTERMONTANE BELT 
In the past two decades over twenty papers of major 
consequence have been written on various geophysical studies 

in the Cordillera. The results of these studies will be 
presented as a series of four discussions, according to the 
type of geophysical information. The discussions will be (i) 
seismology, (ii) potential fields, (iii) time variant 
electromagnetism, and (iv) heat flow. A final summary will 
treat the overall interpretation and its shortcomings. Due 
to limited accessibility, particularily for any long term 
Surveys or data acquisition, most of the geophysical 
information is concentrated in the southern Cordillera, 
(Clowes, 1981). The extension of these interpretations into 
the northern portion of the Intermontane Belt is made purely 
on geological and physiographic grounds. As there are no 
extensive geophysical surveys in the northern Cordillera 
extrapolations of geophysical data or conclusions from such 
Studies to the crust and upper mantle of the Stikine should 
be regarded as highly speculative. 

The seismicity of Western Canada has been treated by 
Milne (1963) and Milne et al (1978). The distribution of 
earthquakes greater than magnitude five is coincident with 
the Insular Belt. These major events give fault plane 


solutions consistent with compression in a NW-SE sense, 
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parallel to the Queen Charlotte - Fairweather transcurrent 
fault system. Milne et al (1978) have mapped the strain 
release in units of quantity of magnitude five events per 
10,000km* per hundred years. On this plot, the zero contour 
of strain release lies mainly within the Coast Crystalline 
Belt. The smaller magnitude events show a pronounced 
regional bias to southern B.C. This reflects only the 
distribution of seismic observatories and the strong need 
for geophysical research in northern B.C. Presumably the 
true distribution of events of magnitude four and smaller is 
more uniform than mapped. Due to the transform fault nature 
of the continental margin, one might expect some dilational 
features on the continental side as a shear-gash fracture 
effect. The likely orientation for such features would 
depend on the local crustal fabric with respect to the 
general NE striking tensional axis. It is possible that the 
Stikine Volcanic Belt and Anahim Volcanic Belts of Souther 
(1977b) are just such crustal features. The presence of hot 
shallow upper mantle and warmer than average continental 
crust would tend to reduce earthquake-style strain release 
in the IMB, (Milne et al, 1978). Regardless of the 
explanation, the IMB appears to be a relatively aseismic 
region. 

The best information on lithospheric structure comes 
from seismology. A single-end refraction profile across the 
southern Cordillera was performed by White et al (1968). 


Wiggens (1976) used refracted body waves from natural 
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seismic events and records from the Canadian Seismic Network 
to examine the upper mantle beneath the IMB. Roebroek and 
Nyland (1975) also looked at Vp data for natural events. 
Surface wave dispersion studies were performed by Wickens 
(1977) for eight stations in south and central B.C. to give 
a regional model for the crust and upper mantle. A more 
recent model incorporating S-wave, P-wave and surface wave 
data is presented by Wickens and Buchbinder (1980). The 
model obtained by White and Savage (1965) is given in Figure 
3-1. Three important statements can be made about the 
Intermontane Belt: (i) the crust is thin (approximately 
30km), (ii) the compressional wave velocity (5.9km/s) is 
rather low by contrast with adjacent crust and the Conrad 
discontinuity is shallower or not seen beneath the IMB, 
(iii) the upper mantle is shallow and has a lower 
compressional wave velocity (7.6km/s) than the mantle 
beneath the adjacent Pacific or the main North American 
continent to the East and South. The data of Dey-Sarkar and 
Wiggens (1976) indicate a low compressional wave velocity, 
approximately 7.8km/sec, and a high Poisson ratio for the 
upper mantle beneath the northern IMB. Roebrock and Nyland 
(1975) examined P wave residuals from over 1500 natural 
events for ray paths that bottomed beneath Western Canada. 
Their travel time residual versus epicentral angle (delta) 
plots for three sub-regions of the Cordillera show a large 
variation in the delta range from 14° to 17°. This would be 


expected for the presence of a 30 to 50km low velocity zone 
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Figure3-1. 
Geophysical cross section of Cordillera taken from Lambert 
and Caner (1965) originally adapted from White and Savage 
(1965). 
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(LVZ); however, their data do not warrant so singular an 
explanation. They find no significant regionalisation of P 
wave velocity structure from north to south. However, their 
data does support some lateral variation including the 
evidence for modern subduction of the Juan de Fuca Plate. 
Wickens (1971) presents surface wave velocities for Canada, 
Figure 3-2. The regular westward decrease in Love wave phase 
velocity at 50km depth is probably indicative of a regional 
thermal or compositional gradient. Such a feature is in 
agreement with the White and Savage (1965) model. More 
detailed work by Wickens (1976) in the Cordillera 
demonstrates the following characteristics of the IMB: (i) 
thin crust and shallow Moho, (ii) no high velocity lid on 
upper mantle such as is found beneath the Alberta Plains or 
the Insular Belt. Nicholls et al (1981) have argued that a 
high velocity lid may still exist. If the difference between 
high velocity lid and low velocity zone were as small as 
0.05 km/sec, and if the lid is sufficiently thin (<1.0km), 
it may escape resolution by low frequency data. 

Regional gravity studies in the Cordillera have been 
made by Garland and Tanner (1957) and Stacey (1973, 1974). 
Kanasewich (1966) has used joint inversions of gravity and 
seismic data to generate consistent models for the crust and 
upper mantle of Western Canada. There is a widespread Bouger 
minimum over central B.C. and the IMB. From crustal 
seismology there is sufficient evidence to rule out the 


normal explanation of crustal thickening, Airy root, as the 
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oe on Figure 3-2, 
ve wav : 
Ie 4) © Velocity Scricture at a depth of 50km from Wickens 
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Figure 3-3. 
Density variation in the Cordillera computed from joint 
inversion of gravity and seismic data; taken from Berry and 
Forsyth (1975). Contoured density in units of kg per cubic 
meter. Negative sign denotes low density regions for crust 
and possibly upper mantle. Roman numerals denote 
physiographic belts: 1. InSular belt, II. Coast crystalline 
belt, III. Intermontane belt, IV. Hinterland belt, V. 
Omineca belt. 
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cause. A model for density variation in the crust and upper 
mantle is given in Figure 3-3 (taken from Berry and Forsyth, 
1975), which shows the Bouger effect as a large low density 
region. Because this coincides with thinner crust, it 
implies a low density upper mantle which could have at least 
three explanations: (i) the absence of a high density lid on 
the upper mantle beneath the IMB, (ii) the presence of a 
thick LVZ with partial melting, or (iii) a less refractory 
upper mantle. This type of potential field information 
cannot distinguish among the three cases. Kanasewich 
presents a Pratt - type model which attempts to account for 
gravity anomalies with gross density variations. This model, 
in Figure 3-4, integrates all three possibilities mentioned 
above. 

Caner (1969a,b) examined aeromagnetic profiles across 
the Cordillera and found long wavelength (greater than 
150km) anomalies west of the Rocky Mountain Trench. He 
presents three possible explanations: (i) hot crust and 
consequently a shallow Curie-isotherm depth, (ii) small 
scale break up of magnetic (and geological) structures, and 
(iii) more silicic crust. It would be quite reasonable in 
the light of the present geological knowledge of the 
Cordillera to accept all three hypotheses. 

Caner (1969a,b; 1970) has also performed time variant 
electromagnetic studies of the Canadian Cordillera by both 
geomagnetic depth sounding and magnetotellurics. The 


inversion of his electromagnetic data as a layered 
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Figure 3-4, 
Model Cross section of crust and upper mantle of the 
Cordillera based on joint inversion of seismic and gravity 
data; taken from Kanasewich (1966). 
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Pgupe: 3-5. 
Comparison of geoelectric sections for world rift areas with 
B. C. Data for Baikal from Gornostayev (1972), for East 
Africa from Banks and Ottey (1974), for Rhine Graben from 
Reitmayr (1975), for B. C. from Caner (1969, 1970, et al 
1969). Data base includes magnetotellurics, geomagnetic 
Gepth sounding and aeromagnetics. Error bars indicate range 
of reported values. World rift areas have a zone of 
anomalously low resistivity (<50ohm-meters) at shallower 
than 30 km depth. 
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resisivity structure predicted lower crustal apparent 
resistivities of about 5ohm-meters, which led Caner to 
speculate that the lower crust is probably volatile rich 
(hydrated ?) and/or partially molten, and that the Moho is 
hotter than 750° with a moderately conductive upper mantle 
which could be partially molten. An alternate explanation 
for such a highly conductive lower crust would require 
temperatures above 900°C, which is considered less likely. A 
comparison of Caner's magnetotelluric sounding data for the 
Cordillera with rift areas of the world is given in figure 
3-5. The geoelectric signature for the IMB is sufficiently 
Similar to established rift areas to speculate that the IMB 
may be an incipient continental rift. 

The western Cordillera of Canada and the IMB very 
probably has heat flow values on the order of 1.8 to 2.5 
H.F.U. (ucal/cm?/sec), (Jessop and Judge, 1971; Jessop and 
Souther, 1980 ), which would compare with the Cordilleran 
geothermal anomaly of the western U.S. (Roy et al, 1968). A 
value of 2.0 was given (Jessop and Judge, 1971) for 
Penticton with values of 2.0 and 2.3 in adjacent Washington. 
The prairies are characterized by values of 1.0 to 1.2 and 
the shield has values like 0.8. The more recently active 
volcanic regions of the IMB no doubt exhibit the highest 
heat flow. The cause of the high western heat flow values 
must lie within the mantle. Geologically unreasonable levels 
of radioactive elements would be required for crustal rocks 


to explain the observed heat flow. To some extent, higher 
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temperatures in the crust and mantle would be consistent 
with all of the previously discussed geophysical data. 
Increased temperatures have a lowering effect on 
compressional wave velocities (Anderson and Grew, 1977). The 
range of laboratory determinations for mafic rocks and 
minerals is reported to be -2.5x10°* to -38.0x10-% 
km/sec/°C. For garnet and forsterite this lowering effect 
shows moderate increases with temperature. Soga (1966) 
derived polynomial expressions for extending these data to 
higher temperatures which show the effect to be non-linear. 
If we use the value of -5x10-* km/sec°C for the temperature 
dependence of velocity, the normal U.M. velocity of 8.5 
km/sec and the anomalous U.M. velocity of 7.8 km/sec, we can 
calculate a hypothetical lateral temperature difference of 
600°C. Although this sort of difference could exist on a 
lateral scale of hundreds of kilometers between normal and 
hot mantle, the actual difference is likely to be more on 
the order of 400°C or less. Roy et al (1968) report 35km 
isotherms of 860°C under the western U.S. and 460°C under 
the eastern U.S. Using 800°C as a temperature estimate for 
the Moho of the IMB, the 8.1km/sec normal mantle under the 
continent could be 600°C, but certainly not as low as 200°C. 
If instead of the -5x10-* km/sec/°C velocity - temperature 
dependence given above, Hughes and Maurette's value for 
Dunite of 38x10-* km/sec/°C is used, the model temperature 
difference would only be 78°C. With data of this quality one 


cannot distinguish between temperature and compositional 
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effects. As was suggested in the discussion of geomagnetic 
induction data, both temperature and composition effects are 
likely to play a part. 

A single cross-section has been constructed to 
Summarize the general geophysical characteristics of the 
Intermontane Belt and its setting within the Cordillera. 
This cross-section is shown in Figure 3-6. In summary, the 
geophysical data for the IMB indicates the following 
generalizations. The crust may have a shallow 0 to 3km upper 
layer of variable low velocity (3 to 4km/sec) sediments. 
Under the IMB the Conrad is shallower and discontinuous 
(poorly resolved). This coincides with possible upwelling of 
the Curie Isotherm (to i8km ?) and presumably higher than 
average continental heat flow values (2.0 H.F.U.). The 
compressional wave velocity of the IMB crust is only 
5.9km/sec by contrast with the 6.1 km/sec of the Coast 
Crystalline Belt or the Eastern Marginal Belt. There is a 
lower crustal (below 10km) zone of high conductivity which 
could indicate volatile enrichment and/or partial melting. 
All of these observations are consistent with a hotter more 
Salic crust for the IMB. The upper mantle of the IMB lacks a 
distinct high velocity lid such as is present 300km to the 
east or west. The upper mantle is characterized by lower 
compressional and Love wave velocities, but the occurrences 
of lateral-compositional variations cannot be distinguished 
from lateral-temperature variations. The most important 


generalization of these physical data is that the IMB is a 
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Figure 3-6. 
Hypothetical cross-section for the crust and upper mant 
beneath northern B. C. from seismic velocities and gr 
Gata. The coast crystalline belt after Johnson et al 
IMB region from Wickens and Buchbinder (1980), Rocky 
Mountains from Merevu et al (1977), density variations after 
Stacey (1973) and Kanasewich (1966). 
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distinct geological region for the entire thickness of the 
continental crust and the upper 40km of the mantle. The 
compositional and thermal causes of volcanism in the IMB are 
certainly related to these pervasive geophysical 
peculiarities. Any realistic tectonic or geological theory 
for volcanism in the IMB must take into account at least 


several of these physical features. 


INTERPRETATIONS OF AEROMAGNETIC DATA FOR THE NAHLIN PLATEAU 

The Nahlin Plateau, including Level Mountain, has been 
mapped aeromagnetically on a 1:50,000 scale by the 
Department of Energy, Mines and Resources. The spatial 
resolution of the maps is sufficient to discern features of 
wavelength greater than 1.6km. The resolving capability of 
the maps for magnetic features on the basaltic plateau can 
be expressed as the average and maximum lateral gradients, 
of 250nT/km and 700nT/km, respectively. Preliminary 
geological interpretations have been made of prominent 
magnetic features, particularly where they relate to 
basement geology, structure and tectonics, and to Late 
Cenozoic basalts. Several anomalous features and general 
trends are portrayed in sketch-map form in Figure 3-7. 

In these high latitudes, above 58°N, a geometrically 
regular vertically magnetized body of rock, having a 
susceptibility contrast to surrounding rocks, should display 
a vertically symmetric anomaly that is either a high or a 


low depending on the sense of the magnetization (up or 
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Figure 3-7, 
Schematic interpretation of aeromagnetic maps for the Nahlin 
Plateau. 
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down). Prior to the Cenozoic there have been structural 
deformations and tectonic movements such that the present 
Orientation of the rocks, with respect to the earth's main 
field, has changed with respect to the original position. 
There are three major types of apparent susceptibility 
Signatures on the Nahlin Plateau: (i) highly magnetic, Late 
Cenozoic basaltic lavas, (ii) moderately magnetic 
cryStalline rocks that make up the local basement 
(plutonics, older volcanics and metamorphics), and (iii) 
weakly magnetic sediments. 

The dominant magnetic strike in the area is WNW-ESE 
which is in good general agreement with broad geological and 
structural trends mapped by the GSC (Gabrielse and Souther, 
1962). Major basement faults such as the Nahlin Fault and 
the Snow Peak Fault are discernable in that they separate 
geological terranes of different magnetic expression. 
Recognition of these faults is based on alignment of paired 
high-low features (second vertical derivative zero 
crossings) and to a lesser extent on horizontal inflection 
of magnetic trends. 

The marine sediments of the Nahlin Trough and the 
Dudidontu Lowlands, northwest of Level Mountain, are 
typified by horizontal gradients below 20 nT/km. Approximate 
calculations (by the methods of Peters, 1949) on the depth 
to the basement below these sedimentary areas indicate a 
sedimentary thickness on the order of 2km. There is also an 


indication of a different susceptibility for the basement 
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rocks on opposite sides of the Nahlin Fault, lending support 
to regional geological mapping (Gabrielse and Souther, 
1962), which indicates the Nahlin Fault to be a major 
crustal structure. 

The magnetic fabric of the igneous and metamorphic 
basement is dominated by subsymmetric features with 
wavelengths less than a few kilometers and intensities in 
the range 250 to 1500nT . The prominent igneous basement 
features appear to be non-vertically magnetized. At least 
five of these anomalies could be explained by bodies that 
were initially-vertically magnetized (in normal field) and 
were subsequently rotated to the southwest. This poses a 
tectonic question as to whether the bodies were originally 
vertically-magnetized and subsequently rotated down to the 
southwest, or whether there was a major shift in latitude 
(more than 15° to the north), of this block of crust afer 
the Jurassic. Symons (1978), on the basis of paleomagnetic 
evidence, has suggested crustal tilting 30° down to the west 
for a section of the Coast Range Batholith a few hundred 
kilometers to the south. This sort of simple tectonic 
rotation is by far the least demanding hypothesis 
energetically. Because the Cretaceous and Cenozoic outcrops 
of the Level Mountain Region are essentially horizontal, any 
hypothetical crustal tilting must have occured between 
Jurassic and Middle Cretaceous time. The possibility that 
the Stikine terrane is allochthonous and has been 


tectonically displaced prior to the Cenozoic cannot be ruled 
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out. 

In addition to the east-west belt of lows discussed 
with respect to the continuation of the Ketchum stock under 
the Level Mountain Plateau, there is another north-south 
belt of lows parallel to Beatty Creek. This trend is neither 
parallel to regional structure nor to the Mesozoic 
intrusions. There exist two possible explanations of this 
feature: (i) this is a basement topographic high and the 
lavas are thin or (ii) a major portion of the plateau 
section here could be reversely magnetized. The second idea 
is unlikely in the light of the development of the lava 
plateau, discussed later, and in terms of the measured 
paleomagnetic section at the Little Tahltan Canyon, which 
shows most of the thickness of the plateau to be normally 
magnetized. The first explanation would make Beatty Creek a 
resequent stream. The Dudidontu to the west and the Tuya to 
the east are exumed or resequent. If the historical drainage 
density was similar to that of today, Beatty Creek is 
probably a resequent stream as well. 

The remaining two magnetic features on the Level 
Mountain Plateau are magnetic highs with amplitude greater 
than 2000nT and wavelength 4 to 8km. These anomalies have 
been profiled, digitized and interpreted by the geometric 
techniques of Peters (1949), Vacquier et al (1951). The 
domical high indicated by a 1 on the aeromagnetic sketch map 
is the most prominent magnetic feature on Level Mountain as 


a whole. Its amplitude against the plateau expression is 
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2200nT and against true regional background is greater than 
2500nT . It is elongate NNW-SSE with a maximum expression of 
8km long by 4km wide. Depth rule estimate derivative methods 
(Peters, 1949), in addition to assumed susceptibility 
contrasts and 2-D modelling by use of a graticule after 
Grant and West (1965) and Peters (1949), indicate a depth to 
the top of the body of 3.46km (2.15miles) below ground 
Surface (ie. 1630m below sea level). These features relate 
to hypabyssal mafic intruSions which are probably intimately 
related to Level Mountain volcanism. The two features occur 
on the Beatty Creek Sheet near the southern end of the 
Central Chain and on the Granite Lake Sheet on the Northern 
Level Mountain Plateau. Both assumed bodies are coincident 
with regions of dyke outcrop, and the first is possibly 
associated with a high level keystone graben which is 
indicated on the cross-section of the geological map. In 
surface outcrop the dykes are quite thin, meters to a few 
tens of meters, and the susceptibility contrast is not 
sufficient to account for even a few percent of the observed 
anomaly. If the initial discussion in this section (on the 
map resolution of magnetic features) is taken into account, 
estimates of the lateral extent of the more prominent 
anomaly indicate an upper surface area of 20km’. If its 
thickness is assumed to be roughly the same as its lateral 
dimensions, the volume can be estimated as 100 to 160km’, 
which is about the volume of an average pulse of surface 


volcanism for the plateau events. In a similar fashion, the 
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2000nT anomaly in the southwest corner of the Granite Lake 
Sheet has been interpreted to be due to two adjacent 
vertically-magnetized basaltic plugs that reach the plateau 
base (957m below plateau surface or 719m above sea level). 
The northwesterly of the two bodies has been estimated to 
have a volume of 80km* and the southeasterly a volume of 
55km*. In terms of petrogenesis, these hypabyssal stocks 
most likely represent high level magma chambers that fed the 
plateau eruptions. The estimated susceptibilities of these 
bodies, based on measurements of Level Mountain rocks, would 
indicate them to be basalt (gabbro) in composition. Neither 
of these bodies can be realistically considered to be the 
residue of protracted crystal fractionation, such as might 
be required to generate the observed salic volcanism. Their 
calculated volumes are too small by at least an order of 
magnitude to account for the observed volumes of salic magma 
and their location is substantially removed from the 


extensive salic piles. 
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CHAPTER 4. GENERAL GEOLOGY OF THE LEVEL MOUNTAIN RANGE 


HISTORICAL GEOLOGY AND FIELD MAPPING 

Field mapping and sampling of the Level Mountain 
Formation was performed by the author during the 1975, 1976, 
1977 field seasons. Stratigraphic sections were examined, 
described, and correlated between the central glacially 
dissected ridges and also between the southern and western 
distal regions of the volcanic shield. Stratigraphic control 
and sampling elevations were determined by a hand-held 
altimeter (Thommens c), accurate to within 15 feet. Mapping 
was done on one inch to one mile air photos with 
approximately 400% ground coverage and on 1:12,500 scale 
topographic base maps. The 1:12,500 scale topographic maps 
were enlarged mosaics made from the 1:50,000 NTS sheets. The 
geologic cross-sections have been constructed along 
composite traverse lines with combined topographic and 
altimetric control. The geologic map of the Level Mountain 
Volcanic Centre, Map 1, was constructed from field maps and 
an air photo mosaic and transferred to a 1:63,000 scale 
topographic base. The actual traverse control on the mapping 
can be ascertained from the traverse and sampling map, Map 
2. The numbered traverses indicate the sampled and measured 
sections. The northern and eastern portions of the shield 
have been interpreted by helicopter traverse, photo geologic 
methods and with aeromagnetic maps. The air photos (NAPL) 


were taken from an altitude of approximately 20,000 feet 
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with a six inch lens taken over the 19 year period from 1948 
to 1967. The aeromagnetic maps are from EMR, 1979. 
Basement Geology 

Traverse and mapping in the basement rocks 
(Pre-Miocene) have been concentrated in three areas: (i) 
Dudidontu River Valley - Ketchum Lake, (ii) Egnel Creek, and 
(iii) Little Tahltan River Valley. The geology varies ona 
local scale of 2km or less with examples of marine sediments 
and metasediments, acid plutonic rocks, porphyritic andesine 
andesites and conglomerates derived therefrom, subaerially 
erupted pyroxene andesites and continental sediments. These 
units will be briefly descibed in the hope of making a 
contribution to the understanding of the complex geologic 
history of the area. The names and ages ascribed to the 
various basement units do not always coincide with the 
preliminary mapping by the GSC (Gabrielse and Souther, 
1962). This reinterpretation of the geology is intended to 
Synthesize existing information with field work done by the 
author. The description of units are presented in order of 
increasing age. 
Continental Sediments - Sustut Group 

Isolated sections of Early Tertiary fluvial and 
lacustrine clastic sediments occur along Classy Creek, 
Mansfield Creek, Tuya River and Nahlin River. Silts and 
shales predominate over sands and thin lenticular coal 
deposits are present along with some occurrences of 


mineralized and petrified wood (G. Sivertz, pers.comm.; T. 
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Mould, pers.comm.). The GSC reports fossils of Paleocene age 
from this unit in the Tuya River (Gabrielse and Souther, 
1962). The unit was likely deposited along a series of 
separated Early Tertiary Basins in respose to some regional 
change in base level, similar to the Deadman River Formation 
on the Bonaparte Lake sheet (Campbell and Tipper, 1971) or 
the Hotnarko Sediments on the Anahim Lake Sheet (Tipper, 
1969). 
Sloko Volcanics 

In the valley of Egnell Creek and in places on the 
Wrathall Plateau are outcrops of subaerially erupted 
calc-alkaline volcanics, see selected analyses in Table 4-1. 
The entire formation is up to 500 feet thick (152m) and 
formed predominantly of green and purple flows of pyroxene - 
andesite and pyroxene - hornblende andesites. The package of 
andesite flows is obliquely cut by dykes from 1 to 3m in 
width of the same material. Capping this pyroxene andesite 
package on the surface of the Wrathall Plateau are local 
occurrences of fine grained basaltic andesite, rhyodacite 
and rhyolite breccia. The rhyolite breccia occurs near the 
northern edge of the Wrathall plateau and appears to 
continue into a plug-like body which probably represents a 
vent area. This entire sequence of calc-alkaline volcanics 
has been subjected to zeolite facies metamorphism and to 
considerable alteration with fracture and joint 
mineralization and infilled vesicles. The groundmass of most 


rocks is devitrified and recrystallized to a matte of 
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Table 4-1. Sloko Volcanics Chemical 
Norms and Trace Elements. 
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$102 52.71 
Ti02 1.52 
A1203 15.30 


Fe203 4.26 
Fed 6.22 
MnO OO. 17 
MgO 4.51 
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Na20 3.43 
K20 1.80 
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Trace element analyses ppm (XRF) 
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zeolites and clays or calcite and quartz. This unit has 
exposures along Egnell Creek, the Sheslay River and to the 
north and west of Ketchum Lake. Near the Dudidontu River, 
aphyric basaltic andesites outcrop with predominant colours 
of grey, red, and brown, with the brown material being very 
intensely altered and friable. Chemical analyses on 
specimens from Egnell, Wrathall and Dudidontu are presented 
in Table 4-1. Structural deformation of these volcanics 
appears to be mainly tilting 10° to 15° to the northeast and 
some block faulting. Some sulfide mineralization is locally 
present near the vent area and in exposures to the west and 
north of Ketchum Lake. This unit (Sloko) does not appear on 
the preliminary mapping of the Dease Lake sheet (Gabrielse 
and Souther, 1962). Instead, the Level Mountain Formation is 
mapped in its place. These rocks are older than the Level 
Mountain formation by tens of millions of years Level 
Mountain Formation and they are distinguishable on the basis 
of their field appearance of weathering in pastel colours, 
calc-alkaline mineralogy, major and trace element chemistry, 
their zeolite facies metamorphism and mild state of 
structural deformation. Judging from the appearance of the 
Eocene and Oligocene volcanism in the central and southern 
part of the Intermontane Belt, these Sloko rocks is older 
Still. It is tentatively assumed to be Uppermost Cretaceous 
or Paleocene in age and roughly correlative with the Sloko 


Volcanics which have been discussed by Souther (1967, 1970) 


and Monger (1968). 
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Tahltan Group (Lower Cretaceous) 

Outcrops of dense green, porphyric andesine andesite 
and indurated conglomerate, with the andesite as the 
predominant clasts, occur in the Little Tahltan River Canyon 
at about 3600 feet elevation. Outcrop of this material is 
also present on Kakesta Mountain to the southeast. The name 
given to this rock by exploration geologists in the area is 
Turkey Track Porphyry due to cumulophyric aggregates of 
andesine. These volcanic rocks and closely associated 
sediments are distinct from both the earlier amphibolite 
facies meta-andesites (Triassic) that outcrop on the 
southeastern side of the Heart Peaks and from the Sloko 
Volcanics (Upper Cretaceous-Paleocene) described above. 
Presumably this unit is representative of the lower Jurassic 
Stuhuni Group, although very similar material from the 
Tulsequah sheet is mapped as Cretaceous (Souther, 1971). 
Sheslay Batholith (Coast Plutonic Complex, Triassic or 
Jurassic) 

Outcrops of multiphase, intermediate to acid, 
calc-alkaline plutonic rocks are found in the vicinity of 
Ketchum Lake and to the southeast of Heart Peaks. Similar 
plutonics outcrop adjacent to the Sheslay River. Predominant 
lithologies are hornblende biotite granodiorite and quartz 
monzonite, but locally pods of diorite and pegmatite dykes 
are exposed east of Ketchum Lake. These plutonic rocks occur 
as a series of nested and cross-cutting stocks with average 


dimensions of about 0.4km. These intrusive rocks have 
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experienced hydrothermal alteration and potassium 
feldspathization in addition to both syn- and post-orogenic 
Structural deformation. The most intense deformation is 
represented by NNW-SSE shear zones to the northeast of 
Ketchum Lake, with intensive chloritization and some 
assorted supergene copper mineralization. The general 
appearance of this plutonic complex is reminiscent of the 
earliest Mesozoic plutonism seen throughout the Intermontane 
Belt, such as the Thuya Batholith on the Bonaparte Lake 
sheet (Campbell and Tipper, 1971). These intrusive rocks 
have definite cross-cutting relationships with greenstones 
and low-grade regionally metamorphosed pelitic sediments to 
the east of Ketchum Lake. They also have cross-cutting 
relationships with amphibolite facies metavolcanics (Stuhini 
Group) along the east side of the Heart Peaks. 
Marine Sediments and Submarine Volcanics 

A sequence of shales and siltstones with green volcanic 
flows outcrops in the Little Tahltan Valley. This unit is 
folded about a NW - SE axis with high angle dips (greater 
than 50°). The northern and southern limits of this map unit 
are high angle faults. Other than the expected effects of 
burial diagenesis and minor calcite veining, the sediments 
are mineralogically unaltered. They could be Triassic in age 
and belong to the Stuhini Group, which would put them in map 
unit 4 as mapped by Gabrielse and Souther (1962), or they 
could belong to map unit 7 and be part of the deepwater 


facies of the Lower Middle Jurassic Inklin Formation. 
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Regional Metamorphics (Paleozoic) 

Low-grade, regionally metamorphosed, pelitic sediments, 
chlorite schists, and greenstones outcrop east of Ketchum 
Lake. These rocks are crosscut by the Ketchum Stock as well 
as having fault contacts with the stock. These rocks are 
isoclinally folded and faulted about a NNW - SSE structural 
trend. The principal metamorphism appears to predate the 
structural deformation, while the structural deformation 
postdates the emplacement of the Ketchum Stock. This does 
not preclude that the metamorphic recrystrallization was 
related to the emplacement; however, the metamorphism 
appears to be regional rather than contact. These rocks do 
not fit easily into the map units of the Dease Lake sheet. 
It is the opinion of the author that these regional 
metamorphics, along with the amphibolite facies basic 
volcanics that outcrop the SE of Heart Peaks, represent 
metamorphosed paleozoics. They could be the same as map unit 
13 (Gabrielse and Souther, 1962), which occurs an roadcuts 
along the Dease River in the northeast part of the Dease 
Lake sheet. The other possibility, although not a strong one 
in the light of contrasting compositions and lithologies, is 
that these are the metamorphosed equivalents of the 
Permo-Triassic Cache Creek-type rocks that outcrop to the 
north of Level Mountain across the Nahlin fault. 

Within the basement geology of the Nahlin Plateau 
Region there is evidence for five or more episodes and 


directions of structural deformation, two distinct periods 
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and grades of metamorphism and five or more episodes of 
igneous activity. All the events took place in the past 
three hundred million years. 

The Nahlin Fault is a major tectonic feature both in 
terms of surface outcrops (Gabrielse and Souther, 1962) and 
crustal structure, as can be seen from the 1:50,000 
aeromagnetic maps. This major fault probably dates from 
Permo-Triassic time with movement continuing through the 
Late Jurassic. The Snow Peak Fault, also WNW in strike, is 
probably related to the Nahlin Fault, although it may be 
younger in part. The major WNW valley of the Tahltan River, 
which has a coincident aeromagnetic expression similar to 
the faults mentioned above, may also be a major upper level 
crustal structure. It has been hypothesized by Gabrielse, 
Souther and Monger (pers. comm.) that the Nahlin fault is a 
fossil transform relating to the accretion of a Stikine 
plate to the Cordillera in Mesozoic time. This could account 
for the pronounced differences in stratigraphy, structure, 
lithology, metamorphic grade and magnetic expression of the 
two sides of the fault. The location of the Level Mountain 
volcanics alongside this fault may indicate a zone of 
pervasive crustal weakness. 

In three places minor faulting was observed, each with 
its own orientation. In the Little Tahltan River Valley, two 
high angle reverse faults striking N100°E were observed to 
cut marine sediments which are either Triassic (Stuhuni 


Group) or Lower Middle Jurassic (Inklin Formation). The 
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faulting is probably Jurassic, although the sense is not 
clearly related to the major Triassic-Jurassic structures 
previously discussed. The Ketchum Stock and the pelites to 
the east are both faulted by a set of thrusts or high angle 
reverse faults with strike N165°E. These are probably 
Cretaceous in age. At Egnell Creek, high angle normal 
faults, of Eocene age (map unit 7 of Gabrielse and Souther, 
1962) bound a downdropped block of Upper 
Cretaceous-Paleocene Sloko volcanics. None of these minor 
Structural trends fits easily into a single crustal fabric 
or deformational scheme. 

The only straightforward point that can be made in 
review of the Pre-Tertiary geology of the Level Mountain 
Region is that it has been an area of repeated tectonism, 
volcanism, and sedimentation, which is to say that the 
Stikine has long been associated with an active continental 
margin. The outcrops in the area attest to the complexity of 
Shallow crustal geology as do the aeromagnetic maps. The 
basement rocks are predominantly calc-alkaline plutonics and 
clastic marine sediments, and aside from the few major 
faults, blocks of similar age and structural fabric probably 


do not extend more than a few kilometers. 


PHYSICAL VOLCANOLOGY 
The Level Mountain range is a low profile, continental 
shield volcano with two principle parts: a flat lying 


basaltic plateau and a composite stratocone cap of bimodal 
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composition. The plateau has pronounced cliff margins due to 
eroSive action of continental ice. Alpine glaciation has 
deeply incised the central statocone to form a series of 
radial U shaped valleys with separating ridges (plate 4-1). 

The Level Mountain shield is made up of sequences of 
thin fluid alkali basalt flows with occasional tuffs and 
localized fluvial or glacial sediments. The basalts 
generally display 2 to 3m thicknesses with prominent 
columnar jointing. These columnar jointed basalts are 
usually made up of 3 to 8 flows so closely spaced in time 
(hours to days) that they cooled as a single unit. Most of 
these cooling units have intermittent basal agglomerates or 
fire-red clays consisting of quenched basalt fragments, 
carbonized plant materials and wood at lower elevations, and 
baked soils and clays. These basal agglomerates have 
pervasive hydrothermal alteration with localized veins of 
calcite and zeolites. Collapsed lava tubes can be seen on 
the present day surface of the plateau and central ridges 
and occasionally in cliff exposures. All tubes in basalts 
had inside diameters less than 2m, see plate 4-2. 
Vesiculation occurs in the upper meter or so of most basalt 
flows and bubble trains or pipe vesicles often extend 
through the flow. Blocky aa basalt flows were only observed 
in two localties: due east of Ketchum Lake on the western 
plateau margin and in the Little Tahltan River valley. In 
both cases the blocky lava flows overlie thick 


volcanoglacial or fluvioglacial units, plate 4-3. Basaltic 
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Plate 4-1. 
Glacially dissected stratocone cap. View to the SE. 
Subhorizontal trachyte flows and tuffs (units 5 to 7) 
dominate the stratigraphic succession of the 
Kakuchuya-Dudidontu ridge. Alpine valley of Kakuchuya creek 


in foreground is downcut to the level of the basaltic 
shield. 
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Plate 4-2. 
Cross section of 2 meter high lava tube in a fine-grained 
alkali basalt-hawaiite. Map unit 4, Little Tahltan Canyon. 
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Plate 4-3. 
150 meter thick pile of unconsolidated fluvioglacial and 
volcanofluvioglacial sediments exposed at the forks of the 
Little Tahltan River on the southern plateau margin. Steep 
dips of up to 70° and folding may have resulted from 
deposition on top of ice, with later slump to their present 


position. 


jnseorq iste 


| choot in yeaa 
| rs 


: eh « 


51 


dykes tend to be linear features. Too little is known about 
basement structures to say in each case whether basaltic 
dyke orientation is the result of basement fabric or high 
level volcano-tectonic control. Some of the more prominent 
dyke features, however, can be explained. The dykes in the 
region of Egnell Creek are parallel to Early Cenozoic 
faulting, where Sloko volcanics are preserved ina 
down-dropped block. The dykes on the northern plateau are 
coincident with a major aeromagnetic anomaly that is 
interpreted to reflect a hypabyssal gabbro body, possibly a 
crystallized shallow magma chamber. Tristanite dykes (1-3m 
wide) in the central chain make a linear feature from below 
the forks of Beatty Creek to Meszah Peak. This i0km trend is 
not parallel to folds or basement faults in the region 
(Gabrielse and Souther, 1962), yet it is difficult to 
imagine a single high level magma chamber of this dimension. 
Three of the thicker flows of the upper plateau unit in 
the westen region formed tumulus structures. Reports by 
Watson and Mathews (1944) and a helicopter survey indicate 
several tuyas on the northeastern plateau surface. Tuyas are 
attributed to be the result of eruption into the base of an 
overlying ice sheet. Other volcano-glacial deposits include: 
(i) palagonite tuff breccia with basaltic balloons and baked 
lacustrine silts, plates 4-4 and 4-5 , (ii) fresh water 
pillow basalts, (iii) cold lahars (till and agglomerate), 
(iv) till cemented by pisolitic siliceous sinter, and (v) 


glacial erratics at the base of flows and tuffs. The most 
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Plate 4-4. 
Palagonite tuff breccia of basaltic composition outcropping 
on the western plateau margin. Vitroclastic bombs and 
glass-rimmed balloons (to right of hammer) stand out in a 
poorly-bedded tuff. 
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Plate 4-5. 
Basaltic volcano-glacial tuff contains boulders of granite 
not indigenous to the volcano, fragments of older 
fine-grained basalt (both to right of hammer), and baked 
lacustrine silts (to left). 
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recent volcanism appears to be feeble bombs (fusiform) and 
Spatter at localized vents in the central region. 

Several physical and textural peculiarities are present 
in the salic flows: 

(1) peralkaline trachyte flows tend to be 2 to 4m thick 
with relatively high phenocryst content, yet judging from 
flow morphology, they are as fluid as basalts in the same 
section, 

(ii) thicker ‘units of comendite and trachyte have 
columnar jointing and in some horizons abundant lava tubes 
(see plate 4-6), 

(iii) aphyric and vitrophyric units often show peculiar 
flow-banding and mottling. This is generally a purple and 
gold blebby banding parallel to the original flow surface. 
Glassy comendite dykes often have fresh cores and white 
crystalline edges. The glass is peralkaline, whereas the 
white crystalline margins are now metaluminous. These 
changes in texture and chemistry are probably due to 
crystallization, devitrification or alteration (Noble, 
1967). 

(iv) aphyric and fine-grained trachytes often show 
onion-skin exfoliation features. In some massive units this 
weathering feature is so pervasive that it masks the true 
flow bedding. It can be mistaken for structural deformation 
and steep dips. The appearance of the exfoliated lavas is 
like a metamorphic schist. The key to understanding these 


features are outcrops on the crest of alpine ridges where 
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Plate 4-6. 
Infilled lava tube from a comendite trachyte of unit 7 on 


the SW face of Meszah Peak. Fluid features of these salic 
flows attest to their low viscosity upon eruption. 
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the exfoliation relates to periglacial and frost features, 

(v) four rock bosses occur in the basal salic unit. All 
of these are thick, columnar-jointed, aphyric masses that 
presently sit in the bottom of alpine glacial valleys. They 
may be neck plugs or other hypabyssal volcanic features or 
they might be the result of eruption of salic magma into 
ice. On the whole, the salic rocks are physically much more 
variable in appearance, thickness, distribution, phenocryst 
content and chemistry than the basalts. 

(vi) there is a high level trachytic stock on the north 
Side of Meszah Peak. This stock has sSub-vertical miarolitic 


Cavities along two zones and well developed joint sets. 


GEOMORPHOLOGY AND GLACIAL GEOLOGY 

After the volcanic landforms, the most striking 
geomorphological features visible on Level Mountain are 
glacial or periglacial in origin. According to Bostock 
(1948) the Stikine plateau was covered by Pleistocene ice to 
an elevation of 6500 feet. The central, volcanologically 
complex, region has been dissected by alpine glaciation to 
form a series of radial U-shaped valleys with separating 
ridges (see figure 4-1). The ridges are sculptured by 
cirques which may have fed alpine glaciers or may have 
entirely postdated them. Many of the cirques have small 
tarns and moraines which adjoin or overlie the lateral 
moraines of alpine valley glaciation. The radial valley 


glaciers generally cut a sharper gradient through the Level 
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Figure 4-1, 
Interpretation of glacial geomorphology for Level Mountain. 


@LACIAL GEOLOGY OF THE LEVEL MOUNTAIN RANGE 
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Mountain Formation than the small primary dips of the lavas. 
The differential hardness of basalts and trachyte flows as 
compared to alternating tuffs and agglomerates gave rise to 
a riegel and basin profile down the valleys. The basins 
often contain shallow, bedrock floored lakes: Dudidontu 
Valley, two lakes; Kakuchuya Valley, two lakes; North Beatty 
Valley, one lake. The remainder of the stepped valley 
bottoms are generally covered with till, alluvial and/or 
mass-wasting deposits that usually mask the geology (for 
example, the South Beatty Valley and the upper valley of the 
Little Tahltan). Both the alpine glacial valleys and the 
plateau margins display rock bosses, whalebacks, and 
drumlin-like hills. The aphyric and relatively homogeneous 
high level intrusSives are more resistant to glacial erosion 
than the surrounding well jointed thin flows and tuffs. A 
hypabyssal trachyte stock, northeast of Meszah Peak, forms a 
whaleback with glacial striae and chatter marks. Massive 
Sills or laccolithic type plugs in the Kakuchuya and Upper 
Beatty Valleys form rock bosses with asymmetric profiles, 
that step down the valley in the direction of ice flow, 
plate 4-7. 

The upper regions of glacial valleys and the ridge 
crests have abundant examples of felsenmeer, sorted frost 
polygons, stony solifluction and patterned ground. On slopes 
the stone polygons are elongate down slope with vertical 
stone bands on the perimeters and interior regions of fines. 


Steep slopes on the central region, as well as on the 
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Plate 4-7. 
Rock boss, sill of columnar jointed aphyric peralkaline 
phonolites in lower (western) end of Kakuchuya alpine 
valley. Valley wall to north comprises comendite, phonolite, 
hawaiite and alkali basalt flows of unit 5. 
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plateau margins, have sliplenses and terracettes. A 
landslide of modest proportion occurred during July, 1977, 
in the Little Tahltan Valley. 

Stratified screes flank most of the domes along the 
central range, see plate 4-8. High gradient snow melt 
Streams have cut through the scree to reveal steep initial 
dips and layering with different degrees of inclination. 
Most ridge crests and cliffs have frost wedging and 
crevassing. In thick columnar units near cliff faces these 
crevasses may be 2 to 3m wide and more than 10m deep. Frost 
shattering and exfoliation of salic flow tuffs gives rise to 
onion-skin textures and apparent cross-bedding. Sometimes 
exfoliation is so severe as to mask the horizontal nature of 
Salic flows. 

Small moraines or rock glaciers occur below the 
cirques. Lateral and end moraines are found at the mouths of 
most glacial valleys (Dudidontu and Little Tahltan). The 
lower reaches of the alpine valleys and the major extent of 
the plateau surface have mudfields and thufur (mud 
hummocks). Frost boils are abundant in mudfields and in 
seasonally dry ephemeral lakes (meltwater ponds). One small 
(4m high x 10m x 15m) pingo-like feature was found near 
Ketchum Lake in a bog-filled, shear feature of the Ketchum 
Stock. The entire southeastern plateau surface (south of 
Nalachaga Mountain) is covered with string bogs, with the 
pushed up ridges forming terraces of elongate bogs across 


the low sloping ground surface, plate 4-9. Steep slopes on 
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Plate 4-8. 
Statified screes, cupola and headed dykes (top right) on 
phonolite dome in Dudidontu alpine valley. Dissected 
stratocone flank in near background. Coast range in far 
background. 
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} Plate 4-9. 
String bogs cover the low gradient plateau surface south of 
Nalachaga mountain, between Beatty and Mansfield creeks. 
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the plateau margin give rise to mudflows and landslides. A 
pronounced variability of dissection and patterning of the 
plateau surface is evident on the air photos. The plateau 
surface has oriented ridges of till, usually a meter or less 
in height, and parallel bands of abundant lakes and elongate 
ponds. This upper plateau unit also has basalts which 
contain erratics, or formed as tuyas, pillow basalts, 
palagonite tuff breccias, and lava flows that overly two 
till units (Pliocene). These features may be found at the 
canyon of the Little Tahltan Valley, and on the northeast 
plateau, where the Dudidontu River crosses the western 
plateau margin, on the western plateau margin above Ketchum 
Lake, and in the lower region of the Little Tahltan Valley. 
This episode almost certainly marks a major continental 
Glaciation. Dating by fission track and Rb/Sr would indicate 
that this continental glacial episode occured around 
4.5MYBP. The striae, till ridges and ablation till fields 
all indicate a southeastern direction of ice flow. This ties 
in nicely with major outwash and eolian deposits to the 
south and southeast of the Level Mountain Plateau. In 
contrast to this, a northerly ice movement is indicated by 
bedrock features in the lowlands west of the Level Mountain 
Plateau. Perhaps these features represent an earlier 
continental ice flow, or even a later one that was thin 
enough to have bypassed Level Mountain without again 
Scouring the plateau surface. If this latter interpretation 


is correct, the time may have been correlative with the 
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valley glaciers which dissected the stratocone. This 
interpretation would make the penultimate continental 
glaciation stronger than the ultimate in the Level Mountain 


aré€a. 


STRATIGRAPHIC DEVELOPMENT OF THE LEVEL MOUNTAIN FORMATION 


Plateau Building Stage 

Volcanism at Level Mountain began in Upper Miocene time 
between six and seven million years ago. The topographic 
relief in the area was mild and the region presently covered 
by the lava plateau probably held two or more mature 
drainage systems Similar to the lowland portion of the 
present day Dudidontu River. The first flows were rather 
volumetric and served to infill those Mid-Tertiary valleys. 
The two major ancestral drainages probably ran north-south, 
one between the Little Tahltan River Canyon and Elm Lake and 
the other on the east side of Beatty Creek through Nalachaga 
Mountain toward Megatushon Creek. Regional drainage for 
these streams was probably to the north-northeast similar to 
the Dease system (to the East) and the exhumed Dudidontu 
Valley (to the west). 

The basal flows rest uncomformably on three extensive 
Pre-Tertiary rock units: (i) marine sediments of Mesozoic 
and possibly some of Paleozoic age, (ii) intermediate to 
acid calc alkaline plutonics, and (iii) upper Cretaceous 


andesitic volcanics. Thick agglomerate sequences and 
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basaltic stream conglomerates are exposed in several places. 
Always low in the Level Mountain section, these represent 
the basal infilling of previous topographic lowlands. Along 
the bottom of the Little Tahltan Canyon a thick basaltic 
stream boulder conglomerate apparently underlies the lowest 
columnar ankaramite flows exposed laterally in the cliff 
section. This Late Tertiary conglomerate rests on a highly 
indurated Jurassic conglomerate whose predominant clasts are 
subangular fragments of prophyritic andesine andesite 
Similar to those exposed on Kaketsa Mountain to the 
southwest. The boulders in the basal Level Mountain 
conglomerate are coarse-grained ankaramites, commonly 
Sub-rounded and up to a few tens of centimeters in size. The 
rounding could easily be accomplished because the coarse 
ankaramites typically display spheroidal weathering or 
decompose to a pyroxene grus; however, clasts of this size, 
Making up a unit 7 to 10m in thickness, could only be 
Carried and deposited in a high energy and fairly major 
drainage system. Similar basal conglomerate/agglomerate 
Outcrops occur where the north fork of the Beatty Creek 
incises the north-east flank of Wolf Bones Ridge and where 
the Megatushon downcuts through the plateau to the northeast 
(Ostensoe, 1960). The character and placement of these 
volcanically derived sediments and certain linear features 
on the 1:50,000 N.T.S. aeromagnetic maps (1979) have 
Suggested the ancestral Tertiary north - south drainage 


reconstruction as discussed above. 
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After the existing topography was infilled, 
construction of the plateau edifice could commence. The 
lowest flows (plateau unit 1, Hamilton and Scarfe, 1977) of 
the plateau are massive ankaramites. These, as well as 
subsequent plateau events, were probably major fissure 
eruptions. The dyke features mapped on the Northern Plateau 
represent such fissures for unit 4, earlier vents have 
probably been covered. Each of the plateau units (I-IV, 
Hamilton and Scarfe, 1977) has up to eight major 
equivolumnal cooling units for an average total volume of 
165km° per eruptive pulse. The flows and flow packages that 
compromise the plateau are essentially horizontal. Lateral 
correlations of 10km or more suggest dip rates of less than 
1° and maximum uncertainin stratigraphy ties of a few 
hundred feet. Because the contacts of the first three 
plateau units are essentially exposed only in steep cliff 
sections, such as the Little Tahltan or Egnell Canyons, 
these three units have been lumped for display purposes on 
the geologic map. 

Each plateau unit is composed of several cooling units 
and each cooling unit is in turn made up of one or more 
flows. These cooling units generally have a thin basal 
agglomerate or "fireclay" and a closely timed series of 
flows which may still retain their own basal crystal-rich 
horizons or upper vesicular regions, yet still cooled and 
Crystallized as a single unit and show pervasive columnar 


joints, plate 4-10. The spacing on the columnar joints is 
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Plate 4-10. 
Alkali basalts and hawaiites of map unit 4, southern plateau 
margin. Ledges develop in basal fireclays and rubbly 
vesicular flow tops between more resistant flows. 7m thick 
unit at base of cliff is comprised of 4 closely spaced flows 
which cooled as a single unit with pervasive columnar 
jointing. 
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greatest for the thickest units and generally decreases 
upwards. These polygonal joint sets are probably tensional 
cooling cracks that propagate downwards as the flows 
crystallize. The time span associated with the flows of a 
Single cooling unit is likely to be on the order of days or 
weeks. Nowhere in the well exposed canyon sections were 
paleosols observed between flows of a major plateau unit, so 
that it can probably be safely inferred that the hiatus 
between succesSive eruptions was a maximum of about a 
thousand years and probably more like a few hundred years. 
In contrast to the equi-volumal, equi-temporal, chemical and 
petrographic similarities of the lavas that comprise each 
plateau unit, pronounced changes occur between one eruptive 
period and the next. The basal flows of each subsequent 
sequence may infill channels or rest on fluviatile or 
lacustrine sediments in areas where the stratigraphic break 
is most pronounced. The lowest flows more typically display 
basal agglomerates and laharic breccias which overly baked 
boggy soils and peats. The charred branches and roots 
exposed in the basal contact of plateau unit II (in Little 
Tahltan and Egnell Canyons) are probably willows and birch. 
This would indicate a Similar environment to the present day 
environment seen on the southern and western plateau 
margins. 

The four plateau units alternate between ankaramite 
flows and/or tuffs and alkali basalt flows. This alternation 


is probably indicative of some cyclical processes in the 
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source region or in the transport mechanism. The plateau 
building period took approximately 2MY using the timing 
correlation obtained from paleomagnetic studies. Evidence 
for an episode of continental glaciation is found during the 
time span of plateau unit 4. The volcanological and some of 
the glacial geomorphology discussed above mostly belong to 
this episode. Souther and Symons (1974) report a till from 
about this time in their paleomagnetic section of Edziza. 
Stratocone Stage 

After the completion of plateau building a fundamental 
change occurred in the type and volume of Level Mountain 
volcanism. The next three eruptive cycles, map units 5 
through 7, were predominantly salic lavas and tuffs with 
only a subordinate volume of basalts (salic to basic volume 
ratio of greater than 4.0). Each of these three salic 
composite pulses had about one-third the volume of the 
preceeding plateau events. The eruptive activity during this 
period seems to be related to a few central orifices, 
particularily for the salic lavas. The longest lived of 
these composite vent areasS waS a Stratocone located to the 
south of Meszah Peak centered in the present alpine glacial 
valley occupied by Kakuchuya Creek. The flows and tuffs of 
Kakuchuya Creek were erupted over a period of about 2MY and 
at itsS maximum during eruptive stage 7, this cone probably 
reached an elevation of about 8000 feet above sea level. 
Contemporaneous with the eruptive activity at Kakuchuya 


cone, there were also major eruptive centers in the west 
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fork of Beatty Creek and near the divide between the Tahltan 
and Dudidontu alpine valleys. During stage 7 volcanism at 
Meszah Peak and in the valley of the north fork of Beatty 
Creek began to take prominence. The shift in eruptive 
activity from Kakuchuya Cone to Meszah Peak, plate 4-11, and 
Gravel Top (another vent 1.5km southeast from Meszah along 
the same ridge), plate 4-12, was coincident with the 
appearance of widespread volcanotectonic features. Faulting 
was observed at two places in the West Beatty Valley and 
both examples were subvertical with a strike parallel to the 
trend of the valley (northwest). The eastern most of these 
two faults occurs on the southeast nose of Wolf Bones Ridge, 
where the plateau basalts on the southwest side of the fault 
are apparently downthrown about two hundred feet and rotated 
to give the flows orientations ranging from N60°E, dip 55°NW 
to N45°E, dip 37°NW. The fault trace was a series of sag 
ponds. The second example occurs higher up the west fork of 
the Beatty on the opposite side of the valley. Here the 
plateau basalts and basaltic agglomerates have orientation 
N50°E dip 71°SE and the fault trace is a stream. These two 
examples alone might be dismissed as slump blocks related to 
glacial activity; however, there are additional examples of 
tensional features. Faulting occurs along the Egnell Creek 
which is a radially oriented graben feature. That this 
region is a zone of structural weakness dates from at least 
Cretaceous time when it was a volcanic-tectonic graben for 


calc alkaline lavas (see cross section at base of map). At 
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Plate 4-11. 
View to NE of Meszah Peak, highest point of Level Mountain 
range. Resistant outcrops (units 7, 8 and 9) of peralkaline 
trachyte, benmoreite and hawaiite flows with comendite flows 
and ash flow tuffs. Dark scoria slope to right is flank of 
Gravel Top, a trachyte-benmoreite scoria cone. 2m wide 
discontinuous tristanite dyke crossing flank of cone points 
to Meszah Peak. Plateau surface (top right) is the drainage 
of Matsatu creek. 
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Plate 4-12. 
View to SE of Gravel Top. Moderately well stratified 
trachyte-benmoreite scoria cone on crest of alpine ridge SE 
of Meszah Peak. 
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present the northwest margin of this feature has a series of 
linearly aligned post-glacial hawaiite dykes, sag ponds and 
keystone features. It is also likely that the central chain 
with its calc alkaline rhyolite domes (map unit 7c) is 
located on a subsided block as indicated in cross section. 
The most abundant salic lavas and tuffs of unit 7 are 
comendites which are markedly less peralkaline than are the 
pantellerite trachytes of units 5 and 6. Additionally, map 
unit 7 also contains domes and breccia pipes of 
calc-alkaline rhyolite, see plate 4-13, and areally limited 
flows and dykes of metaluminous tristanite. To 
Stratigraphically tie the basal units of the central ridge 
requires that the ridge be a keystone graben which has been 
down dropped 200 feet or more with respect to the outcrops 
on Meszah Ridge or Wolf Bones Ridge. The limited occurrence 
of rhyolites and other calc-alkaline lavas in this central 
chain is probably related to a high level crustal magma 
chamber and to some degree of shallow tectonism. The 
coincidence of present alpine valleys with the one time 
active centres of extrusion may simply be related to the 
predominance of tuffs, high primary dips, and pronounced 
thermal alteration at these locations. However, the 
possibility of radial volcanotectonic grabens cannot be 
totally discounted. 

The latest eruptive activity at Level Mountain has been 
a series of basalt flows, spatter and dykes that occurred in 


at least 3 separate episodes over the past million years. 
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Plate 4-13. 
View to south from summit of Kakuchuya-Beatty dome. Dark 
comendite dyke in foreground above tarn. Domes of central 
chain are dominantly comendite and rhyolite tuffs with 
Subordinate obsidian flows and ignimbrites. Low elevation 
background (top right) is southern plateau margin. Coast 
range in far background. 
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Selected hand specimens from this recent activity are shown 
in plate 4-14. These basalts are phenocryst-rich hawaiites 
with prominent plagioclase megacrysts in the dykes along the 
Egnell and Wrathall Plateau and inclusions of gneiss and 
troctolite on Meszah Peak. The volume of these is very 
subordinate by comparison with the plateau basalts and is 
more comparable to the basalts of units 5b and 6b. In that 
most of the erosion responsible for the radial glacial 
valleys is probably Pleistocene in age, up to 60% of these 
lavas may have been removed by glacial activity. 

According to the paleomagnetic evidence for timing, no 
period in the history of Level Mountain has seen so little 
eruptive activity as the last million years. Although some 
eruptive features definitively postdate glaciation, their 
volume is so small that Level Mountain should be classed as 
dormant or extinct. A review of the energetics of Level 
Mountain volcanism, given in a subsequent section, supports 


this viewpoint. 
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Plate 4-14, 
Hand specimens of recent hawaiites and basalts. Top left, 
breadcrust bomb; Top centre, vitroclastic scoria; Top right, 
polygonal "micro"-columnar jointed vesicular hawaiite; 
Bottom left, fusiform aerial bomb; Bottom right, ropy lava 
from hawaiite flow top. 
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CHAPTER 5. PETROGRAPHY OF THE LEVEL MOUNTAIN LAVAS 


INTRODUCTION 

In a discussion of the petography of the Level Mountain 
lavas it is important to describe all of the rock types, 
textures, mineralogies and crystallization sequences and to 
relate them in a reasonable genetic framework. A summary of 
petrographic data and representative modes are presented in 
Table 5-1. Modes include both phenocryst and groundmass 
phases. The petrochemistry, norms and detailed mineralogy 
for most of these samples are presented in subsequent 
sections. The petrographic discussion of major lava types 
emphasizes the flows and fine grained tuffs rather than the 
agglomerates. This primarily reflects a volume bias; 
however, very few samples were collected of the 
agglomerates, lahars and coarse grained pyroclastic rocks 
due to their dubious value as petrochemical or genetic 
indicators. The petrographic discussion is the result of 
observing over six hundred thin sections and reconnaissance 
microprobe work on over one hundred polished sections. The 
norms are more restricted as chemical analyses were limited 
to about two hundred and fifty rocks that were considered 
sufficiently fresh and representative. The important point 
to consider throughout the discussion of each lava type is 
that volumetrically there are only two major and important 
Magmas, alkali basalt and peralkaline trachyte. Every other 


lava type ultimately bears some genetic relationship to one 
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or both of these. 


ANKARAMITES 

Ankaramites and olivine-rich basalts (ANK in Table 5-1) 
are the least glassy of the Level Mountain lavas. Two 
textural varieties are common: pandiomorphic subophitic 
intersertal, see plate 5-1, and pandiomorphic protoclastic 
cumulophyric intersertal, see plate 5-2. The protoclastic 
cumulophyric variety may have phenocryst contents as high as 
45% by volume, although 15 to 20% is more typical. The 
groundmass olivines and titaniferous augites of this type 
have an intergranular relationship to the felty plagioclase 
laths. Oxides are cubic to subdendritic. Intersertal areas 
are less common than in other basaltic flows. The common 
primary infillings are devitrified brown glass. The 
anakaramites have essential phenocryst olivine which is more 
forsteritic than Fo,,;. Groundmass olivine may be as low in 
magnesium as Fo,;. The size distribution for olivines from 
phenocryst to goundmass is usually seriate with the largest 
phenocrysts being a few millimeters. The forsteritic 
phenocrysts are hypidiomorphic to allotriomorphic and may 
have concentric or patchwork optical zonation and occasional 
anhedral cores. Embayed olivines from the flows of plateau 
unit 2 frequently contain inclusions of brown spinel. 
Optical and microprobe analyses on these 
xenocrystic-appearing olivines and spinels reveal them to be 


chemically variable even within a single flow. Lamellar 
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eiace: 5-7; 
Ankaramite (H)(X.P.L., 2.97x2.07mm) with subophitic zoned 
titaniferous augite (sides) intersertal plagioclase An63, 
and subhedral olivine altered to iddingsite. 


Pine D2 . 
Ankaramite (8/23-7/3785)(X.P.L., 2.97x2.07mm) protoclastic 
pandiomorphic with polycrystalline clots of olivine having 
brown spinel inclusions (xenocrysts). Groundmass contains 
euhedral plagioclase, titanomagnetite and olivine with 
intersertal augite and brown glass. 
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exsolved diopside sometimes accompanies forsteritic olivine 
in these xenocryst clots, see photomicrograph of sample PAR, 
shown in plate 5-3. The olivines of the ankaramites are 
considered to be in part xenocrystic fragments of the 
basaltic source region and in part of the earliest 
crystallizing phase. Although most ankaramites have olivine 
as the only phenocryst phase, the following multiphase 
assemblages have been observed: (ol, opx), (ol, opx, pl) 
(ol, pl), (ol, cpxy pl), and (ol, cpx). While phenocryst 
olivine and orthopyroxene may show fracture and imperfectly 
developed crystal forms, they do not show a mutual reaction 
relationship or textures indicative of clear cut sequential 
crystallization. Neither do they show a marked reaction 
relationship with the mesostasis. By contrast, whenever 
clinopyroxene phenocrysts occur they show pronounced 
reaction coronas. Plagioclases usually show reaction rims 
and fritted cores. The crystal forms and aggregate textures 
of the plagioclase and clinopyroxene phenocryst clots 
Suggest rather that they are cognate plutonic xenocrysts. 
The clinopyroxene of the groundmass of the ankaramites is 
invariably titaniferous showing brown or purplish colours in 
plane polarized light and zonation to higher titanium 
contents at the crystal rims. The crystal form of the 
groundmass clinopyroxene within a given flow is usually 
constant, although overall the texture ranges from 
intergranular prismatic, to feathery acicular, to 


subophitic. Skeletal magnetites sometimes form as parallel 
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Biate 5-3. 
Alkali basalt (PAR)(X.P.L., 2.97x2.07mm) Protoclastic 
xenocryst or cognate xenolith clots of olivine and lamellar 
exsolved diopside in a groundmass of labradorite, augite, 
titaniferous magnetite and glass. 
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poikilitic inclusions within subophitic augites. 

Plagioclase feldspar compositions are in the 
labradorite to bytownite range with the groundmass 
plagioclase usually being An,;, to An... Phenocryst feldspars 
tend to be tabular and twinned but not markedly zoned. 
Groundmass feldspar is lathlike to anhedral, normally zoned 
and often complexely twinned. 

The predominant crystallization sequence is olivine 
(with or without diopsidic pyroxene and spinel), 
plagioclase, titanaugite, Fe-Ti spinel. This is reflected 
both in texture and serial size range. Olivine 
crystallization is pre-eruptive while the crystallization of 
other phases is usually not. Primary intersertal glass was 
10% by volume or less. This is usually altered to blue green 
chlorite or mixtures of chlorite and zeolites. Vesicles are 
rare in ankaramites and where present have been infilled by 
calcite. The ankaramites are generally more chemically and 
mineralogically altered than the alkali basalts or 
hawaiites; but not sufficiently so as to mask their 


crystallization history or chemically distinctive nature. 


ALKALI BASALTS 

Alkali basalts (AB in Table 5-1) are the most variably 
textured lava type as well as the most abundant. Flows 
predominate over tuffs and the tuffs, once hypohyaline, are 
now severely altered. Flows may be very fine grained and 


nearly aphyric or phenocryst laden panidiomorphic seriate or 
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medium to coarse grained diabasic (Sub-ophitic) in texture, 
see plate 5-4. The occurence of phenocrysts and phenocryst 
clots (xenolith clots?) is very common. For an example see 
plate 5-5. Forsteritic olivine, pale green diopside, 
bytownite to labradorite plagioclase and iron titanium 
oxides occur as phenocrysts. Of these olivine is the most 
common. Phenocryst clots of clinopyroxene and 
clinopyroxene-olivine are the most abundant aggregate types. 
Plagioclase and clinopyroxene phenocrysts usually show 
twinning and multiple concentric zonation. Chemical zonation 
is usually in the normal sense but reset and marked by a 
zone of fine inclusions. Non-zoned types in crystal 
aggregates uSually show more optical homogeneity and often 
have triple-point textures. The groundmass has essential 
plagioclase, titaniferous augite, and iron titanium oxide, 
see plate 5-6. Intersertal brown glass or the alteration 
products calcite, chlorite and zeolites are also present. 
Olivine is not always present aS a groundmass phase. When it 
does occur it shows a marked compositional break with the 
phenocryst olivine, such aS Fo,, groundmass, contrasted with 
Fo,, phenocryst. Occasionally in the alkali basalts with a 
more intergranular groundmass, melilite rods were observed. 
The clinopyroxene of the mesosStasis was invariably 
pink-brown titaniferous augite with ophitic or subophitic 
texture. Phenocryst clinopyroxene may be either titaniferous 
augite or green diopside. Plagioclase phenocrysts tend to be 


labradorite, An,; to An ;,, while that of the mesostasis 15 
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Plate 5-4. 
Fine-grained alkali basalt (16/19D)(X.P.L., 2.97x2.07mm) 
with subophitic subhyalo-ophitic texture. Zoned (An55 to 
An45) and twinned euhedral plagioclase laths, titaniferous 
augite, basaltic glass and Fe-Ti oxides. 
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mrace 5-5. 
Sample (PAL)(X.P.L., 2.97x2.07mm) Alkali basalt with 
microxenoliths of pyroxenite (diopside>hypersthene). 
Groundmass glass shows alteration to chlorite. 
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Plate 5-6. 
Groundmass of alkali basalt (PAB)(X.P.L., 1.90x1.33mm) with 
subhedral olivine microphenocrysts, magnetite cubes, and 
laths of plagioclase with intersertal titaniferous augite. 
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more sodic, in the andesine range, An;,. to as low as An;;. 
It is the author's opinion that the only truly 
equilibrium phenocryst phases are olivine and plagioclase or 
olivine and clinopyroxene. The common embayment (resorption) 
of plagioclase and clinopyroxene and reaction coronas 
indicate disequilibrium. Even when these phenocrysts are 
concentrically zoned they are probably cognate xenocrysts 
and sometimes they are genuinely foreign. The normal 
crystallization sequence is olivine, plagioclase, 
clinopyroxene, and opaques, with eruption occurring while 
plagioclase is on the liquidus. The common occurrence of 
reset oscillatory zonation, concentric zonation and 
compositional differences between phenocryst and groundmass 
plagioclase indicates that the liquidus phase and its 
composition is very sensitive to thermochemical variables. 
The changes from phenocryst diopside to mesostasis 
titaniferous augite and the presence of phenocryst ilmenite 
and groundmass titanomagnetite are particularily indicative 
of changes in oxygen fugacity and in the partitioning of 
titanium into various sites. The remelting of phenocrysts 


seems to indicate changes in magmatic heat content as well. 


HAWAILITES 

Hawaiites (HAW in Table 5-1) tend to be protoclastic 
pandiomorphic mesohyaline and are more prone to be 
xenocrystic and xenophyric than the alkali basalts. The two 


types are obviously genetically related. Hawaiites tended to 
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be vesicular pandiomorphic cumulophyric seriate with a 
pilotaxitic mesostasis, see plate 5-7. The mesostasis is 
noticeably dominated by plagioclase (andesine) with variable 
amounts of titaniferous augite, lron-titanium oxides and 
brown glass. Olivine and alkali feldspar may be present in 
the mesostasis but neither is essential. Alteration of the 
groundmass glass is less pronounced than in the more basic 
lavas. The occurrence of large plagioclase phenocrysts, and 
tabular megacrysts to a few centimeters in size was most 
common in the hawaiites. Hawaiites often contain reverse and 
complexly zoned plagioclase and sometimes plagioclase 
phenocrysts of up to three distinct geneses. An example of 
the large, plagioclase-bearing type is shown in plate 5-8. 
Although groundmass plagioclase was usually andesine, the 
phenocryst (xenocryst?) compositions had an overall range 
from andesine to anorthite. Because hawaiites have many 
xenocryst phases, they can contain such diverse phenocryst 
types as enstatite and sanidine in a single flow. While the 
groundmass may represent an important "magma type", the 
reaction relationship to the phenocryst assemblages 
indicates the petrogenetic complexity of the hawaiite. 

The occurrence of hawaiites spans the entire history of 
Level Mountain volcanism, although they are usually a 
Subordinate lava type in terms of abundance. During the 
first three plateau periods they were subordinate to alkali 
basalt and ankaramite. The upper plateau unit 4, as sampled 


in the Little Tahltan Canyon, consists primarily of 
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Plate 5-7. 
Vesicular hawaiite (8/20-2/4410)(X.P.L., 2.97x2.07mm) with 
Seriate euhedral phenocrysts of Fo82 olivine and An50 
plagioclase in mesohyaline mesostasis. 
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Plate 5-8. 
2.97x2.07mm)Plutonic textured xenocrysts of 
plagioclase (An40) and microphenocryst of titaniferous 
augite from fine-grained hawaiite R5. Groundmass is 
oxide-rich with plagioclase, titaniferous augite and brown 
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hawalites. This may or may not be true of unit 4 as a whole. 
While hawaiite is the most abundant mafic lava type during 
the stratocone formation, it was quite subordinate in volume 
and areal extent to the peralkaline salic types. 

Texturally and optically the hawaiites are more complex 
than the other mafic lavas. Judging from their variable and 
often disequilibrium phenocryst character and the relatively 
frequent occurrence of xenocrysts and xenoliths, the 
hawaliites are probably not a primary magma type. Their 
Origin most certainly involves subsequent petrogenetic 
processes such as crystal fractionation, contamination, and 


magma mixing. 


TRACHYBASALTS 

Trachybasalts are the least volumetric and the most 
variable mineralogically and chemically of any rock type at 
Level Mountain, see Table 5-1, Section D. They were usually 
pyroclastic and agglomeratic with phenocryst and fragment 
contents typically around 40% by volume. Mineralogically and 
chemically these rocks are gradational with alkali basalts 
and with the tristanites and trachytes. Trachybasalt 
agglomerates are found in the alpine glacial valley of the 
Dudidontu and at the eastern end of Wolf Bones Ridge. These 
"flows" contain about 30% altered glass and the 
phenocryst/xenocryst assemblage may include hyalosiderite 
Olivine (Fo,, to Fo;;,) hedenbergitic clinopyroxene, 


basaltic hornblende, plagioclase (labradorite to 
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oligoclase), anorthoclase, and magnetite. The occurrence of 
phenocrysts of basaltic hornblende is restricted to 
trachybasalts and tristanites. The olivine is usually in 
aggregates, sometimes with green hedenbergite, and always 
has a reaction corona of iron oxides. Large inclusions of 
apatite are common. The appearance of these mafic 
intermediate clots suggests that were derived by the 
disaggregation of an alkali ferrogabbro, conceivably a 
cognate plutonic member of the Level Mountain suite, see 
plate 5-9. Both agglomerates and flows contain resorbed 
phenocrysts of intermediate plagioclase and frequently have 
glomeroporphyritic aggregates of anorthoclase, pyrite, and 
aegerine augite or sodic ferrohedenbergite. Occasionally the 
clinopyroxene phenocrysts have cores of titaniferous augite 
and rims of green hedenbergite. The phenocrysts in 
trachybasalt are always in reaction with the groundmass. The 
mesostasis, when crystalline, is hyalopilitic to pilotaxitic 
with laths of andesine (An,, to An,,.), Sodic pyroxene and 
yellow-green glass. The chemical classification of these 
disequilibrium rocks is difficult. They may be over or 
undersaturated with respect to silica and may vary from 
anorthite to corundum normative. Following the nomenclature 
of Irvine and Baragar (1971) these rocks may have such 
diverse names as mugearite, benmoreite, tholeiitic andesite 
or peraluminous phonolite. Modal names would certainly 
include trachybasalt and the vulsinites and ciminites of 


Italy (Washington, 1896; Tilley et al, 1964). Textural and 
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Plate 5-9. 
Fragmental trachybasalt flow (LMIII30f£)(P.P.L., 
2.97x2.07mm). Fragments include euhedral plagioclase, 
clinopyroxene and dark basaltic glass. The groundmass is 
eutaxitic with light and dark areas of different 
crystallinity. Interpreted to be a product of explosive 
magma mixing. 
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mineralogical evidence suggests that these lavas are the 
result of magma mixing of basic and sSalic types with 
incorporation of some cognate and foreign rock material from 
the crust. One trachybasalt flow had glassy inclusions which 
were presumably peralkaline due to their sparse phenocrysts 
of green sodic pyroxene, ilmenite and anorthoclase. Although 
the name trachybasalt is not an accepted term in the sodic 
alkaline basalt series, it is used to emphasize their 
peculiar character. In some areas of the world such as the 
French alkaline province (Brouse and Varet, 1966) and the 
Gregory Rift (Goles, 1976) the importance of intermediate 
magmas such as mugearites and benmoreites has been pointed 
out. Goles makes a special case for the petrogenesis of 
peralkaline salic lavas to be dependent on crystal 
fractionation from mugearite-benmoreite types. For Level 
Mountain, however, the petrographic evidence is strongly in 
favor of trachybasalts being a rather minor mixing 


derivative of the dominant basic and salic magma types. 


TRACHYTES, PANTELLERITES, AND COMENDITES 

Trachytes are the dominant rock type from the 
Stratocone. In Table 5-1, P. Trachyte refers to peralkaline 
trachyte of pantelleric affinity and C. Trachyte refers to 
peralkaline trachyte of comenditic affinity. The two most 
common are porphyritic hypocrystalline pilotaxitic and 
porphyritic eutaxitic. Phenocrysts are usually zoned and 


rimmed or resorbed. Depending on the crystallinity of the 
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groundmass feldspars the textures range from orthophyric to 
felsophyric. The groundmass ferromagneSian minerals may be 
prismatic, Subhedral lanceolate or poikilitic. A few 
outcrops of hypabyssal stocks, dikes and sills of this 
composition are more correctly termed peralkaline syenites 
and have aplitic to granophyric textures. The more 
holocrystalline examples contain free quartz in the 
groundmass. 

Phenocryst types in order of abundance are alkali 
feldspar, plagioclase, sodic clinopyroxene, fayalitic 
olivine, aenigmatite, Fe-Ti oxides and rarely biotite. The 
occurrence of two to four phenocryst types iS common in a 
Single flow and usually this assemblage is prevalent for a 
sequence of flows. The alkali feldspars are the single 
essential phase of the trachytes. They may be euhedral or 
resorbed and fritted and often occur in glomeroporphyritic 
aggregates or fragmental plutonic textured clots with mafic 
phases and/or quartz, see plate 5-10. Plagioclase 
phenocrysts have a compositional range from oligoclase to 
albite and they may occur with alkali feldspars but are not 
essential by themselves. Iron rich olivine, in the range 
Fo,; (ferrohortonolite) to Fo,, fayalite, is one of the two 
most common ferromagneSian phenocrysts. It is usually 
euhedral and rimmed or altered and often in clots with 
feldspar, sodic pyroxene and iron titanium oxides. When 
altered, it commonly has rims of sodic pyroxene or 


riebeckite. Olivine is characteristically only a phenocryst 
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Place 2-10. 
Orthophyric, xenolithic pantellerite flow (25/5e)(X.P.L., 
2.97x2.07mm). Cognate xenolith of plutonic textured alkali 
Granite comprised of quartz, microcline and aegerine with an 
aenigmatite rim. 
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phase, although occasionally pure fayalite is part of the 
groundmass assemblage. Sodic clinopyroxene is the most 
common ferromagnesian mineral both as a phenocryst or part 
of the groundmass. Phenocrysts have optics in the range from 
sodic ferrohedenbergite to aegerine augite. They are usually 
pleochroic in shades of green. Phenocryst clinopyroxene may 
only be mildly sodic and be a pale clear green in plane 
polarized light. This type usually has a discrete rim with 
dusty titanomagnetite inclusions and a more acmitic 
composition. The other variation has an intense green 
colouration which dominates the high birefringence. This 
type (oikocryst) often has a peculiar leopard frog 
appearance due to the abundant rounded inclusions of 
titaniferous magnetite. This is probably a reaction texture 
where a pyroxene plus oxide assemblage is pseudomorphous 
after some prior ferromagnesSian phenocryst phase, perhaps 
amphibole or garnet. While the phenocryst clinopyroxene may 
be fresh or altered, that of the groundmass is invariably 
fresh. The groundmass clinopyroxene tends to be more sodic 
when it is the dominant mafic of the mesostasis and is 
accompanied by titanomagnetite. Conversely, when the 
groundmass assemblage includes aenigmatite in addition to a 
sodic amphibole, the pyroxene is only mildly sodic and has 
nearly pure hedenbergite composition. Iron - titanium 
oxides, either ilmenite or titanomagnetite, are present as 
inclusions in fayalite or clinopyroxene. Titaniferous 


biotite was observed in a nepheline normative trachyte 
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(phonolite) dyke and flow of map unit 5. Aenigmatite is very 
prevalent as a groundmass phase, although occasionally it 
occurs aS a phenocryst. The aenigmatite phenocrysts are 
euhedral and prismatic and may be the only ferroan phase or 
may be subsequent to sodic pyroxene. The deep red-brown 
pleochroism tends to mask the high birefringence (0.07). The 
groundmass aenigmatite may occur as subhedral prisms 
although subophitic or poikilitic masses are more common. A 
photomicrograph showing typical groundmass textures is given 
Beaplate. 5-1 1. 

The paragenesis of aenigmatite is a rather complex and 
controversial topic. The stability of aenigmatite has been 
studied experimentally by Ernst (1962), Marsh (1975), 
Lindsley et al (1971), and Thompson and Chisholm (1969). All 
of these investigations show that aenigmatite is present in 
a relatively restricted field in T-f£0, space. When 
aenigmatite occurs as the only crystallizing ferroan 
Silicate this can be interpreted to represent the joint 
conditions of peralkalinity and low f0,. However, at Level 
Mountain the aenigmatite usually occurs in a mesostasis 
assemblage with phases that are far more oxidized such as 
alkali amphiboles. Usually iron-titanium oxides are absent. 
Since the two groundmass fractions are often only separated 
by tens of microns, this leads to conceptual problems in the 
interpretation of £0, conditions. Where the peralkaline 
trachytes are somewhat deuterically altered, a golden yellow 


astrophyllite may partially or wholly replace the 
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Plate 5-17. 
(P.P.L., 0.48x0.33mm) Pilotaxitic groundmass of peralkaline 
trachyte (7/11). Dark green phenocryst (side of photo) and 
euhedral high relief prismatic green brown phase in 
groundmass are sodic pyroxenes. Deep reddish-black anhedral 
aenigmatite (centre), clear dark blue riebeickite, light 
pale green arfvedsonite, golden astrophyllite and clear 
laths and anhedra of alkali feldspar and quartz. The 
groundmass of these holocrystalline peralkaline trachytes is 
typically heterogeneous, as if there were steep diffusion 
gradients on a scale of a few tens of microns. 
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aenigmatite. The groundmass assemblage typically has two 
feldspars in sub-equal amounts. The plagioclase is 
oligoclase to albite and normally zoned. The alkali feldspar 
is Sanidine, typically occurring as Square cross sectional 
stubby prisms. Quartz may occur in rapikivi-type 
intergrowths in the eutaxitic trachytes or in the irregular 
Spaces between feldspar laths. Possible groundmass mafic 
minerals include fayalite, sodic pyroxene, aenigmatite, 
alkali amphiboles (riebeckite, arfvedsonite), and 
astrophyllite. Rosenbuschite was recognized as a groundmass 
phase in a few flows (both by petrography and microprobe 
analysis). Common accessory minerals include baddeleyite or 
zircon. Occasionally the groundmass contained a 
cathodoluminescent phase with high contents of Ce and rare 
earths. 

The pantellerites and comedites are transitional to the 
peralkline trachytes and are not so abundant. Owing to their 
Similar mineralogy and chemistry these more silica rich 
peralkaline types are probably differentiates from the more 
common peralkaline trachyte magmas. Pantellerites at Level 
Mountain are almost exclusively holocrystalline, see 
photomicrograph shown in plate 5-12, while comendites occur 
as glassy dykes, flows, and tuffs. This situation is the 
converse for other peralkaline localities such as 
Pantelleria (Villari, 1974) or Fant-ale (Gibson, 1972, 1974) 
where the more strongly peralkaline pantellerites occur as 


pitchstones. The peralkaline eruptives of the Great Basin in 
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Plate 5-12, \ 
(X.P.L., 0.60x0.42mm) Groundmass of holocrystalline 
pantelleritic trachyte (25/5b). White areas with Sutures are 
Quartz and microgranite. Sanidine shows straight edges or 
rectangular cross sections. Some riebeckite films gives 
quartz a bluish shade. High relief prisms with anomalous 
green and pink birefringence are rosenbuschite. Aenigmatite 
prisms are euhedral with deep reddish brown to black 
colouration that masks their birefringence. 
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the western U.S. tend to be comendite tuffs (Noble and 
Parker, 1974). Petrochemically the pantellerites and 
comendites are both peralkaline rhyolites. 

Pantellerites have more than 12.5% mafic minerals in 
the norm. Comendites tend to have agpaitic indices very 
close to unity and higher quartz contents in the norm. There 
is a tendency for the comendites to be preferentially 
associated with oceanic islands and pantellerites with 
continental rifts. At Level Mountain the pantellerites were 
more coarsely crystalline, more ferromagnesian-rich 
counterparts of the eutaxitic trachytes. Common phenocryst 
assemblages were anorthoclase-aenigmatite or 
anorthoclase-fayalite. When fayalite occurred, it had 
reaction rims of riebeckite. Aenigmatite phenocrysts 
typically showed patchy oxidation to iron-titanium oxides. 
Aenigmatite and riebeckite were the most common groundmass 
mafics. Comendites tended to be aphyric or sparsely 
phenocrystic pitchstones or tuffs. The phenocryst assemblage 
was either anorthoclase, fayalite, and ilmenite or 
anorthoclase, sodic ferrohedenbergite, and ilmenite. In the 
few holocrystalline types there was essential groundmass 
Quartz and opaques. The groundmass was almost entirely 
quartzofeldspathic with minor sodic pyroxene or riebeckite. 

A summary of the typical crystallization history of 


these peralkaline salic magmas is given in table 5-2. 
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‘Table 5-2. Summary of Crystallization Sequences for 
Peralkaline Lavas 


1. Fayalite + Anorthoclase on liquidus (dry) 


2. Fayalite oxidizes to magnetite + SiO,(melt) + heat 
(causing anorthoclase resorption) 


3. Fayalite + Melt go to acmite clinopyroxene 
OR 
Soda-Iron Amphibole 
(with build up of pH,O prior to eruption) 


NOTE: Eruption may proceed after 1, 2, or 3. 


4, Groundmass crystallization usually proceeds dry at 
low £0, conditions with fayalite or 
aenigmatite; or at locally determined f0, and H,0 
content, due to diffusion lag and undercooling, 
to get aenigmatite, ferrohedenbergite, alkali 
amphibole. 


5. With late deuteric alteration get astrophyllite. 
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PHONOLITES, TRISTANITES AND RHYOLITES 

The relatively uncommon silica undersaturated lavas are 
restricted to map units 5 and 7 and the metaluminous lavas 
are restricted to map unit 7. Phonolites occur both as tuffs 
and glassy to fine-grained flows. The tuffs, usually yellow, 
white, or buff in appearance contained anorthoclase and 
sparse altered phenocrysts of analcime or sodalite. The 
flows were grey in appearance and relatively aphyric, with a 
groundmass of alkali feldspar, nepheline, opaque and glass. 
These rocks have normative nepheline in the range 4 to 6% 
with plagioclase in the oligoclase - albite range as the 
dominant feldspar. Petrographic examination of phonolite 
flows from map unit 5 in the Meszah Peak section revealed 
phenocryst amphibole and biotite. The phonolites could 
represent some fractional derivative of the more common 
basic undersaturated melts, such as alkali basalt and 
hawaiite, or could be a product of magma mixing of basalt 
and trachyte (particularly for the tuffaceous varieties 
which contain xenocryst clinopyroxene and plagioclase 
derived from basalt). Similar phonolites have been reported 
from alkaline centres such as St.Helena (Baker, 1969), the 
Auvergne (Varet, 1969), the Canary Islands (Tenerife) 
(Ridley, 1970), and the Azores (White, 1979). A review of 
the origin of phonolite lavas is given by Wright (1971). 

Tristanite occurred as flows and dykes. Both quartz and 
nepheline normative varieties were noted as well as 


metaluminous and peraluminous types. The texture is 
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glomeroporphyritic pilotaxitic hypohyaline. The phenocryst 
assemblage includes andesine, anorthoclase, chrome diopside 
(confirmed by microprobe), basaltic hornblende, biotite, and 
titanomagnetite. The groundmass always shows a high degree 
of deuteric alteration. A xenolith of gneiss was found in 
one dyke, so the metaluminous nature may be in part due to 
crustal assimilation by a more typical trachytic melt. The 

- low agpaitic index is not the only petrological peculiarity 
of the tristanites. These flows and dykes are markedly 
potassic which contrasts with the remainder of the salic 
rocks (sodic). The high water content could be inherent, as 
indicated by the hydrous phenocryst amphibole and mica. The 
length and spatial distribution of the tristanite dykes 
suggests that they could be related to high level 
volcano-tectonic features as previously discussed. These 
tristanite magmas undoubtedly have a complex genesis, 
perhaps as a mantle derived trachyte melt they suffered some 
high level crustal contamination, yet were able to retain 
early phenocrysts from depth. It does not seem satisfactory 
to derive the tristanites by crystal fractionation processes 
or by contamination of hawaiites or trachytes. 

Rhyolites are mildly metaluminous and usually contain 
only phenocryst alkali feldspar and magnetite. Pyroclastics 
and glassy varieties are more common than any other salic 
lava type. Rhyolites find their greatest distribution along 
with minor occurrences of calc alkaline dacite during the 


period of unit 7. The rhyolite flows show extensive 
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hydrothermal alteration and occasionally sulfide 
mineralization. The rocks are uSually fine grained and both 
devitrified glass or aplitic textures are common. One glassy 
comendite dyke has a fresh core, yet the hydrothermally 
altered white selvedges were rhyolite, see photomicrograph 
shown in plate 5-13. The rhyolites as flows and tuffs are 
often closely associated with comendites. For the most part 
these lavas and tuffs probably were comendites that have 
lost their slight peralkalinity through crystallization or 
alteration reactions as discussed by Noble (1967). The 
remainder are probably the composite offspring of trachyte 


differentiation and high level contamination. 


INCLUSIONS: OCCURRENCE AND PETROGRAPHIC CHARACTERISTICS 

The majority of lava flows on Level Mountain, 
irrespective of chemical type, are aphyric to fine-grained 
and relatively free of phenocrysts. It is reasonable to 
assume that these fine grained flows have been extruded at 
or near their liquidus temperatures. OcaSionally inclusions 
of both cognate and xenolithic character have been observed 
for a variety of Level Mountain lavas. Their distribution is 
rather widespread but their actual occurrence is certainly 
the exception rather than the rule. Where they are found, 
with the exception of one rhyolite vent area, they never 
comprise more than a few percent by volume of the host lava. 
A summary of petrography and modes for inclusions is given 


in Table 5-3. Plate 5-14 shows a photograph of some of the 
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GreLeuo-13. 
(P.P.L., 0.78x0.55mm) Petrographically this dyke selvedge 
(LMI3) is a comendite, with white euhedral alkali feldspar, 
white anhedral quartz, patches of inky blue riebeckite and 
green arfvedsonite with a few percent pale tan devitrified 
glass. Chemically however this is a rhyolite. 
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inclusion types. 

Megacrysts, cumulophyric aggregates, clots ad 
intergrowths comprise the cognate types. In flows and dykes 
of alkali basalt and hawaiite, single plagioclase and 
pyroxene megacrysts were found that exceed five centimeters 
in size. The most typical example is an egg-shaped 
plagioclase - pyroxene intergrowth which has coarse grain 
size, interlocking plutonic texture, and a rounded shape 
with a size ranging from 2.5 to 8 centimeters. These vary 
from dominantly plagioclase to dominantly pyroxene. The 
pyroxene is usually black and of vitreous appearance with 
conchoidal fracture more frequently seen than cleavage. It 
is eaSily mistaken for basaltic glass. The density of these 
plagioclase - pyroxene inclusions varies from 2.64 to 3.24 
grams per cubic centimeter which covers the range from 
neutral buoyancy to approximately 20% greater than their 
host basalts. The peralkaline trachyte flows of map units 5a 
and 6a occasionally contain lithophysae or cognate clots of 
alkali feldspar and sodic pyroxene with a very open texture 
and acicular intersertal crystals of alkali amphibole. The 
texture of these clots is usually diktytaxitic with actual 
voids between crystal faces. While the plagioclase-pyroxene 
eggs of the basalts might be plutonic autoliths or even 
accidentally incorporated fragments of the source region, 
the open trachytic clots are obviously of a high level 
Origin. The voids and automorphic crystal shapes are 


probably indicative of crystallization in the presence of a 
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Plate 5-14. 
Hand specimens of inclusions: top left, partially melted and 
disaggregated granite gneiss (light coloured) in dark 
hawalite; top right, calc silicate nodule from trachyte 
host, with crystalloblastic calcite and epidote; centre 
left, granitic gneiss in vesicular tristanite, no reaction; 
centre, granitic gneiss (no melting) in hawaiite; right 
centre, embayed eucrite nodule with 
yellow-weathering,glassy-devitrified rim from a peralkaline 
trachyte flow; bottom, partially melted embayed and 
dissagregated granitic gneiss (white) in hawaiite(dark). 
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separate gas phase and reduced pressure. 

The occurrence of incontrovertibly foreign inclusions 
has been used to indicate upper mantle origins of lavas and 
to estimate the compositional variation of the upper mantle 
(Sobolev, 1977). At Level Mountain the occurrence of foreign 
inclusions in fluid lavas is restricted to three 
petrochemical types: basalts, peralkaline trachytes and 
tristanites. The Hawaiites of map units 8-10 on Meszah Peak 
contain xenoliths of anorthosite, gabbro, and troctolite 
(allivalite), see photo in plate 5-15, as well as partially 
fused and quenched fragments of granite gneiss, see photo in 
plate 5-16. Both basalts and peralkaline trachytes have been 
found to contain alkaline ultramafic inclusions having 
diopside and spinel as the dominant phases with variable 
amounts of olivine, amphibole and feldspathoids, see plate 
5-17. One dyke in the Egnell Canyon had xenoliths of augite 
peridotite. An altered ankaramite flow of map unit 5b, in 
the Kakuchuya Valley, was found to contain xenoliths of 
garnet wehrlite (cpx, ol, gt, and sp). Contemporaneous 
basalts from the Stikine Canyon and Castle Rock occasionally 
contain lherzolite nodules (Littlejohn and Greenwood, 1974). 
Deformation lamellae and kink bands in the olivine and 
pyroxene crystals of these ultramafic xenoliths may either 
be the result of explosive shock upon eruption (J. Krupica, 
pers.comm.) or deformation in the source region. 

Inclusions of diorite and hornfelsic sediments were 


found in an alkali basalt flow of the Upper Plateau, unit 4, 
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Plate 5-15. 
(X.P.L., 2.97x2.07mm)Plutonic textured troctolite 
(Allivalite) xenolith (MPM) sampled from hawaiite flow on 
Meszah Peak. Twinned plagioclase dominates with euhedral to 
subhedral olivine and diopside, note fractures in all 


grains. 
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Plate 5-16. 
(X.P.L., 2.97x2.07mm) Partially fused granitic gneiss 
fragment from hawaiite on Meszah Peak. Showing grain 
boundary reaction and interstitial rapikivi texture for 


Quenched partial melts, seen as white and dark areas of 
Quartz and feldspar intergrowths. Mineralogy includes 


plagioclase, alkali feldspar, quartz, biotite and glass. 
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Plate 5-17. 
(P.P.L., 2.97x2.07mm) Alkali peridotite (8/31-85), with 
granular metamorphic-textured green diopside, fractured 
high-relief olivine, with overgrowths of tan to brown 
anhedral hornblende showing good amphibole cleavage, black 
anhedral spinel, saussuritized yellow and white low relief 
plagioclase and clear-granular patches of nepheline. 
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in the canyon of the Little Tahltan River. Granite 
inclusions from basalts at MezSah Peak show petrographic 
evidence of partial melting, reaction and quench textures. 
Armoured ultramafic nodules were noted at three separate 
locations in peralkaline trachyte flows: in unit 6a on the 
east end of Wolf Bones Ridge, on the Dudidontu dome and on 
the ridge that separates the Dudidontu and Kakuchuya 
valleys. The nodules characteristically were about 1cm in 
size and rounded with a yellowish reaction rim up to 2cm in 
thickness. The texture was coarse grained and idiomorphic. 
The mineralogy was dominantly anorthitic plagioclase and 
brown titaniferous augite with inclusions of forsteritic 
olivine, brown apatite, opaques, amphibole, mica and 
accessory interstitial nepheline, see plate 5-18. These 
nodules are both alkaline and ultramafic having a 
composition intermediate between eucrite and essexite. It is 
uncertain whether these nodules are representative of source 
material, a fractionated source residua, or a reaction 
product. However their peculiar mineralogy, texture, and 
reaction rims rule them out as fractionates from the 
trachyte magma. The most abundant inclusions in the 
trachytes are phenocryst or xenocryst clots of alkali 
feldspar, with or without oligoclase and quartz. These 
frequently have sutured internal boundaries, resorption 
outlines and alkali-rich rims, all of which may be an 


indication that they are partially assimilated crust. 
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Plate 5-18. 
(8/6-12n) from 


Eucrite/essexite 
Interlocking plutonic texture 


PR Psl., 2e97X%2.07mm) 
light and dark grey 


peralkaline trachyte flow. 
with twinned cleaved brown augite, 
feldspathoid, brown mica and black magnetite. 
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Although the inclusions are not volumetrically abundant 
at Level nbtinea we they represent quite a variety of 
lithologies. They are associated with several magma types 
and the age of the flows that carry inclusions spans most of 
Level Mountain's history. Petrogenetically the inclusions 
indicate three things. Ultramafic nodules such as 
lherzolite, peridotite, wehrlite, troctolite, anorthosite, 
and eucrite-essexite are usually interpreted as being part 
of the Upper Mantle (Varne, 1970; Basu and MacGregor, 1979; 
Mercier and Nicolas, 1975). They are common inclusion types 
in alkalic magmas and the variety of inclusions in alkalic 
lavas is one of their hallmarks. That these occur in the 
Level Mountain basalts and peralkaline trachytes, in 
addition to other lavas of the Intermontane Belt, indicates 
a mantle source for these magmas. 

The accidental inclusions of older crytalline rocks 
such as diorite, granite, and gneiss indicates the 
possibilities of crustal magma chambers and some crustal 
contamination to the primary alkalic magmas. These foreign 
rocks as well as hornfelsed sediments are found both in 
basalts and also in the low volume, late occurring, 
metaluminous types such as tristanite and rhyolite. A small 
degree of crustal contamination would likely have very 
different effects on basSic and salic magmas dependent upon 
which phase(s) were on the liquidus. In a mildly 
nephelinitic alkali basalt, Si-Al contamination could 


enhance the crystallization of calcic plagioclase and 
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possibly make the magma hypersthene- or quartz-normative. 
While in the case of a comendite melt, the same Si-Al 
contamination could cause the magma to become metaluminous. 
Many of the inclusions are clearly of a different 
texture but of sufficient chemical affinity to term them 
cognate. Into this class fall the pyroxene-plagioclase 
clots, the alkali gabbros and possibly the essexite - 
eucrite inclusions. These cognate types allude to a complex 
plutonic history in addition to the already complex 


surficial manifestations of volcanism. 


SUMMARY OF PETROGRAPHY AND PETROGENESIS 

The basaltic lavas are of alkaline affinity and 
mineralogically similar to the alkaline olivine basalt and 
sodic basalt types. The dominant phenocryst minerals which 
are likely candidates for low pressure crystal fractionation 
arguments include forsteritic olivine, calcic plagioclase 
and calcic pyroxene. The ankaramites and alkali basalts 
could be related by variable olivine and clinopyroxene 
fractionation or accumulation. Fractionation of olivine, 
clinopyroxene and plagioclase could explain the derivation 
of hawaiite from an alkali basalt primary melt. However, the 
complex zonation and reaction textures of the hawaiite 
phenocrysts suggest that the actual petrological processes 
may be more complex. The pronounced zonation and great 
Variability of the plagioclase phenocrysts may partially 


account for the existence of both nepheline and hypersthene 
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normative compositions. Where calcic plagioclase and olivine 
are the dominant fractionating phases, a nepheline-normative 
hawaiite could result as has been demonstrated 
experimentally by Yoder and Tilley (1962). Variation in the 
composition of the fractionating plagioclase, plus the 
effect of clinopyroxene, could cause the trend for silica 
enrichment to dominate, producing the hypersthenic type. Due 
to the occurrence of xenoliths and xenocrysts, simple 
fractionating arguments should not be expected to account 
for all of the observed mineralogical or petrochemical 
variation. 

The peralkaline trachytes and the more silica-saturated 
pantellerites and comendites seem to fall into two groups 
characterized by different phenocryst assemblages. Alkali 
feldspar and iron titanium oxide are ubiquitous phenocryst 
phases. Sodic palgioclase may be ancillary to anorthoclase 
but it has not been observed alone. Fayalitic olivine and 
sodic clinopyroxene are the dominant mafic types. Usually 
they do not occur together. When they do, the fayalite is in 
a reaction relationship with the liquid, while the pyroxene 
is not. Fractionation of alkali feldspar plus iron titanium 
oxide plus a mafic silicate would tend to lower the colour 
index and raise the peralkalinity. When the dominant mafic 
is fayalitic olivine, residual calcium would remain and the 
Na/K ratio would tend to be enhanced. The effect of this 
could perhaps be seen as a more calcic groundmass 


Plagioclase with the occurrence of calcic groundmass mafics 
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such as clinopyroxene or amphiboles. If the dominant mafic 
was a sodic pyroxene, the more effective calcium removal 
would lead to the crystallization of groundmass aenigmatite. 
The effect of having two sodium bearing fractionating phases 
could enhance the subsequent role of potassium. In this 
case, the proportion of groundmass sanidine would likely be 
higher, and the groundmass plagioclase would likely be more 
sodic. 

The derivation of the low volume intermediate magma 
types in the trachybasalt-mugearite-benmoreite range 
probably involves magma mixing. Phenocryst types, usually 
characteristic of either the basic or salic lavas, are 
Sometimes found together in these disequilibrium textured 
lavas. Inclusions of one type of melt in the other have been 
observed in both senses. One flow from the central range was 
a mixture of basaltic and trachytic lava. In this instance 
eruption occured too swiftly for intimate admixture to 


result. 


q 7 


i{ = 


Ra , : 
| atime ant? $i .water 


emlas wutsbe 4: 


PP 42Rhane wm ic 


dn jh , 
€ Ga 740 
t.- ¥ he 
or sad 

if 22 @; 

4 wa 
‘ iy hm Ae 


] 


'S10m0s60 ~ SIT 490 


¢ jnyzbonedas 


vr 


od 
4 


tiie 1 sacaiatabh ie 


nae Gag polved @ 


i j 7 7 


> 
> 
o 


ioGm OF 9. dreary 7 


si Lievayzs sts 38 


* 1% 
~ = «a Be 
25 ove) 5 
3 


CHAPTER 6. CHEMICAL MINERALOGY OF THE LEVEL MOUNTAIN LAVAS 


OVERVIEW OF THE MAJOR MINERAL GROUPS AND MICROPROBE ANALYSIS 
All mineral analyses were performed on the University 
of Alberta ARL electron microprobe fitted with an Ortec 
system for energy dispersive analysis. Using the energy 
dispersive facility, a qualitative identification of 
component mineral phases and their constituent elements was 
made for over two hundred representative lavas. In a 
reconnaissance sense, the variations in major and minor 
element contents were compared for all mineral assemblages. 
Using anhydrous mineral standards and four hundred second 
counting times, over ninety quantitative analyses were 
produced for the major rock forming minerals and a few of 
the accessory phases. The method of data reduction was that 
of Smith and Gold (1979). All data were collected and 
interpreted quantitatively for twenty-one elements: Na, Mg, 
ReaueeD. oS. «GhoaKkecCag*SeyaTiyoV) «Ce, (MnyaFéarGo, Bi, 2Cu, 
Zn, plus Zr, and Ba. This list of elements covers the range 
of atomic numbers eleven through thirty plus forty and 
fifty-six. The wavelength dispersive crystal spectrometers 
were routinely used to verify the presence of minor elements 
and elements not on the list. Additionally, the output of 
the wavelength spectrometers were used in reconnaissance to 
examine the distribution of minor and trace elements and to 
examine the sense of chemical zonation. Oxygen, which is the 


most abundant mineral forming element, cannot be 
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quantitatively analysed by present energy or wavelength 
dispersive techniques and was calculated by difference. In 
an analysis expressed as element weight percent, the 
difference between the reported oxygen value and the 
theoretical amount of oxygen which would be associated with 
the metals in their appropriate valences is a measure of the 
reliability of the analysis. 

In the following sections on systematic mineralogy the 
tables are arranged by mineral type and presented as element 
weight percent, as advocated by Pinsent and Smith (1975), 
and structural formulae calculated on the basis of n 
oxygens. There are several justifications for presenting 
mineral compositions from microprobe analyses as element 
weight percent, rather than as equivalent oxides. The first 
lies in the iterative form of the ZAF correction procedure 
for probe analyses, with direct comparison of element counts 
for standards and samples. Since oxygen is not directly 
analysed, or compared between sample and standards, it is 
the least reliable element and should not be convolved with 
the other elements. Inclusion of oxygen results in the 
spreading out of any error in totals or error in the fit to 
average atomic number. The compositional variations for most 
minerals of interest involve substitution of different 
metals into specific coordination sites within an oxygen 
framework. The parameter of importance in mineral 
geothermometry is the variation in a specific element, not 


its oxide, and that element's concentration (activity) into 
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Specific coordination sites. The expression of mineral 
analyses as weight percent oxides iS a carryover from the 
days of gravimetric analyses. Analyses in this form must 
always be recalculated before they can be used for 
thermochemical or petrogenetic arguments. The structure 
formula calculation was iterative with partition of iron to 
Satisfy charge balance and stoichiometry. The sample numbers 
refer to rocks. The approximate rock type is given for ready 
comparison. Each mineral section closes with a comparison of 
analysed compositions to mineral compositions predicted by 
normative calculations on the host lavas, along with 
discussions of the relevant petrogenetic roles of the 


mineral. 


OLIVINES 

Representative olivine analyses are given in Table 6-1. 
Olivine compositions range from nearly pure forsterite 
(Mg,SiO,) to pure fayalite (Fe,Si0O,) with greater than 97% 
of the compositional variation occurring in this 
two-component solid solution series. Minor element contents 
always include calcium and manganese in subequal amounts 
Such that 1 to 2% molecular glaucochroite (CaMnSiO,) is 
present in all structural formulae. In the 
hortonolite-ferrohortonolite range, molecular manganese 
exceeds calcium such that molecular tephroite (Mn,SiO,) must 
be present. A component of molecular monticellite may be 


present in the chrysolite-hyalosiderite compositions but is 
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not essential. 

Phenocryst olivine for the entire range of Level 
Mountain lavas is in the chrysolite-ferrohortonolite range. 
More forsteritic olivines occasionally occur in basalts as 
phenocrysts and in ultramafic nodules. Compositions near 
pure fayalite are restricted to the groundmass of 
peralkaline salic lavas and alkaline syenite sills. Their 
presence here is considered to be the combined result of a 
peculiar crystallization history and restricted conditions 
of £0, imposed by the composition and physical conditions of 
groundmass crystallization. Olivines, when present, are 
usually the most magnesian phenocryst mineral in a given 
lava or paragenesis. In terms of major and trace element 
variation trends, olivine fractionation shonte have 
pronounced effects on Mg, Mn and ni. 

The microprobe analyses and estimates of olivine 
compositions compare favorably with petrographic 
determinations. The detailed comparison of actual olivine 
chemistry to norm calculations reveals several surprises. 
Normative olivine chemistry for nodules and basalts is 
usually reasonable and within 10% mol Fo. However while 
olivine phenocrysts are very common in the basalts, the 
chemical analyses of the basalts may be hypersthene or 
quartz normative (depending on the iron oxidation state). In 
this case pyroxenes are calculated instead of olivines. The 
comparison of olivine to whole rock chemistry and norms for 


Salic lavas is very poor. While actual fayalitic olivine 
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exists as phenocrysts, microphenocrysts, or in the 
groundmass of salic lavas their norms show iron in oxides, 


acmite, ferrosilite, or hedenbergite. 


CLINOPYROXENES 

Clinopyroxenes, although not universally present, are 
the most widespread of the ferromagnesian minerals in Level 
Mountain lavas. Diopside and titaniferous augite dominate 
the basalts, while salic and intermediate lavas have 
ferroaugite-ferrosalite-ferrohedenbergite types. A 
collection of selected clinopyroxene analyses is presented 
in Table 6-2. The elements Ca, Mg and Fe account for greater 
than 94% of the X and Y site positions. The remainder are 
predominantly accounted for by Al, Ti, Na, and Mn. The 
distribution of the Level Mountain clinopyroxenes and their 
corresponding whole rocks 1s presented for the system 
Wo-En-Fs in Figure 6-1. The range of compositions covers: 
40.3-45.5 Wo, 24.6-44.5 En, 11.9-34.1 Fs, all expressed as 
mol%. The clinopyroxene compositional variation is close to 
the diopside-hedenbergite join when represented in the 
tetrahedron diopside - hedenbergite - acmite - johannsenite. 
There is no good correlation between the normative diopside 
component and the actual Diopside content of pyroxenes in 
the rock. 

Sodium and the acmite component are always more 
abundant than Mn and the johannsennite component for the 


entire range of compositions. Other elements which may be 
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Figure 6-1. 
Level Mountain lavas and clinopyroxenes in the pyroxene 
trapezoid (Wo:En:Fs) (atomic ratio). 
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present include Al, Ti, Cr and ferric iron. Atomic aluminum 
exceeds titanium for all clinopyroxenes. The most titanium 
rich phase is always an oxide; however, the effect of 
clinopyroxene fractionation on titanium variation cannot be 
ignored. All phenocryst clinopyroxenes have more atomic 
aluminum than sodium, so clinopyroxene fractionation will 
cause an increase in the agpaitic index and the tendency 
towards peralkaline residua. Enavopprouenes from 
intermediate and salic melts have a lower aluminum content 
than those from basaltic melts. The concentration of 
aluminum in basaltic melts is equivalent to or slightly 
lower than that in intermediate and salic melts. This 
antipathetic relation between aluminum concentration in the 
melt and aluminum concentration into the clinopyroxene 
suggests a difference in the chemical potential of aluminum 
into the clinopyroxene phase between the two dominant melt 
types. This may be related to differences in the assemblage 
of aluminum-bearing phases, to composition and concentration 
Gifferences for other elements in the melts, or to 
fundamental differences in the pressure - temperature - 
oxidation conditions under which these different 
clinopyroxenes crystallized. As to the first speculation, 
clinopyroxene is never the most aluminous phenocryst phase 
in any of the Level Mountain lavas. It is always associated 
with a feldspar, plagioclase in the basalts and anorthoclase 
in the salics. As to the second speculation, the fundamental 


differences between basalts and salics are obvious but the 
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compositional effects on the changing role of aluminum in 
the melt structure and its activity in the crystallizing 
phases are not obvious. The third possibility (which implies 
a higher pressure origin for basalts than salic lavas) is 
the most likely. This possibilty is dealt with ina 
Subsequent section including calculations on geothermometry 
and geobarometry. 

An exception to the aluminum discussion above is the 
clinopyroxenes of the tristanites. These pale green, 
eight-sided, euhedral phenocrysts are rich in aluminum with 
moderate levels of chromium and titanium. The calcium and 
aluminum levels are comparable in the clinopyroxene and the 
whole rock. In a normative sense, these pyroxenes are 
dominated by hedenbergite, calcium-tschermakite and 
diopside. 

In the salic rocks, clinopyroxene is the most 
Significant calcium-bearing phase. The concentration of 
calcium in clinopyroxene is two orders of magnitude greater 
than calcium in coexisting alkali feldspar, and is also 
Significantly higher than calcium in the whole rock. There 
is no straightforward relationship between sodium content of 
the melt and sodium content of clinopyroxene. In fact, the 
most variation in molecular acmite into clinopyroxene occurs 
for the high sodium salic melts (about 6.0% Na.,0) and is 
equally unrelated to peralkalinity. The aegerine-acmite type 
pyroxenes are generally restricted to the groundmass of 


peralkaline lavas. Analysis #11 is from the groundmass of a 
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eutaxitic pantellerite flow. In a formula sense, this sample 
is 91% acmite and 6% hedenbergite with 3% johannsennite and 
jadeite. A tracing of the spectrum of this unusual acmitic 
pyroxene is given in figure 6-2. 

The typical phenocryst clinopyroxene of the peralkaline 
trachytes and related lavas is sodic ferrohedenbergite. 
These is a positive correlation between sodium content of 
the melt and manganese content of the clinopyroxene. 
However, this manganese expressed as molecular johannsennite 
does not follow the variation of hedenbergite or acmite 
components. 

Clinopyroxenes occuring as phenocrysts are almost 
always zoned. In basalts the clinopyroxene cores are more 
magnesium rich, while rims and groundmass clinopyroxene were 
higher in iron and titanium. Zoned phenocrysts from salic 
melts showed higher iron, manganese and sodium contents for 
crystal margins. In the salic rocks the groundmass 
clinopyroxene may or may not be more sodium rich than 
coexisting phenocryst clinopyroxene. Here the presence of 


other sodic ferromagnesian phases play the deciding role. 


FELDSPARS 

A collection of representative feldspar analyses is 
given in Table 6-3 for plagioclases and 6-4 for alkali 
feldspars. The sample number and rock type are given, along 
with the weight percent elemental analysis and structure 


formula calculated on the basis of eight oxygens. For ease 
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Figure 6-2. 
Energy dispersive XRay spectrum for Na-Fe-Si microphenocryst 
from comendite ignimbrite (9/5-107/6600). Note minor content 
of Zr, Ca, Ti, Mn. From Na/Si ratio, possible compositions 
include: aegerine, rosenbuschite and astrophyllite. 
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in labelling graphs and diagrams, the analyses and sample 
designations are referenced by column numbers. All analyses 
are for phenocrysts unless otherwise noted. In most 
instances, the within grain and grain to grain variations 
were qualitatively examined using short counts and small 
Sweep scans. The actual analyses were done using a twenty 
micron square or larger sweep area on the most abundant 
compositions. The resulting analyses are taken to be more 
representative for crystal fractionation arguments than 
point analyses would be. The large scans and collected 
groups of twenty-five second counts were also done to avoid 
problems of sodium excitation or volatilization. 

In a molecular sense, 98% of the feldspar compositions 
are accounted for by the anorthite, albite and orthoclase 
end members. The most important component beyond those 
mentioned is barium (as celsian). 

Plagioclase compositions from microprobe analyses agree 
well with optical determinations. Phenocrysts from hawaiites 
often have multiply reset or reversed concentric zonation. 
The groundmass plagioclase of these disequilibrium rocks may 
be more sodic than either phenocrysts or whole rocks (see 
analysis from 8/20-2/4307, upper plateau unit). 

Figure 6-3 shows plagioclase-whole rock pairs in the 
atomic system calcium-sodium-potassium. It is evident that 
potassium is only a minor constituent of the plagioclase and 
that the whole rock and the melt are enriched in potassium 


during plagioclase crystallization. The distribution of 
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Figure 6-3. 
Level Mountain lavas and feldspars in Ca, Na, K atomic. 
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whole rocks as seen here could not be simply related by 
fractionation of their phenocryst plagioclase alone because 
the trend of whole rock variation is nearly perpendicular to 
the feldspar vector. plagioclase alone. Presuming the 
basalts to be related by fractionation, the controlling 
phase would have to be a calcium rich phase such as 
clinopyroxene, with only minor effect due to the feldspar. 
Alkali feldspar analyses presented are for phenocrysts 
only. The anorthoclase phenocrysts tend to show patchy 
chemical zonation (in sodium, potassium) on a scale of 30 
microns or less. The zonal differences in Ab and Or are 
Slight and not nearly as pronounced as the difference in 
composition between phenocrysts and groundmass. The 
groundmass almost always has microlites of nearly pure 
Sanidine. Groundmass alkali feldspar is always more 
potassium rich than phenocryst alkali feldspar. Anorthoclase 
phenocrysts may have rims that are either more sodic or more 
potassic than their cores. The presence of two distinct 
groupings is apparent in figure 6-3. The grouping with the 
higher calcium content represent the low volume metaluminous 
Salic lavas. In this set the calcium of content of 
phenocrysts is not significantly lower than in the rock as a 
whole, and the sodium contents are higher. The second group 
are the alkaline and peralkaline salic lavas. With one 
exception, the phenocryst feldspars are more sodic and less 
calcic than their whole rock systems. The two groups could 


possibly be related by the addition of intermediate 
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plagioclase or sodic plagioclase and amphibole. Within 
either group, fractionation of the indicated feldspars 
cannot account for whole rock variations. For further 
elaboration of this and other fractionation arguments refer 
to the subsequent section on whole rock chemistry. The 
peralkaline group variation could be explained by the joint 
effect of alkali feldspar and a calcic phase (hedenbergite). 
The detailed compositional variation of alkali 
feldspar-phenocrysts and whole rocks has been examined by 
mol ratio plots involving sodium, potassium, aluminum, 
calcium and their various combinations. No obvious 
mathematical relations could be discerned. Nicholls and 
Carmichael (1969) made two rather definite statements 
relating alkali feldspar compositions and peralkaline 
residual glasses. Those statements were: (i) K/Na feldspar 
is greater than K/Na whole rock and (ii) feldspar 
compositional variation shows pronounced restriction with 
incresing agpaitic index. The first conclusion is only true 
for Level Mountain sample 25/5e, in fact the converse is 
generally true. Their second conclusion appears to hold for 
Level Mountain lavas. 

Feldspars are the most abundant phenocrysts for all 
Level Mountain lavas but their compositions, if involved in 
fractionation at all, cannot alone account for trends in 
whole rock variation. There are three distinct whole 
rock-feldspar sets: (i) basalt-calcic plagioclase, (ii) 


peralkaline salic-sodic anorthoclase and (iii) intermediate 
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metaluminous-sodic alkali feldspar, with or without sodic 
plagioclase. In all three cases, the whole rocks show trends 
of calcium variation, while the feldspars do not have 
Significantly different calcium levels from whole rocks. The 
most likely phases to account for the calcium variation 


would be calcic pyroxenes or possibly amphiboles. 


IRON TITANIUM OXIDES AND RELATED MINERALS 

Microprobe analyses of iron titanium oxides for 
representative Level Mountain lavas are presented in Table 
6-5 for spinels and Table 6-6 for ilmenites. Analyses were 
obtained using a focussed beam spot, unless otherwise 
indicated. Structure formulae have been calculated for 
ilmenites on the basis of four oxygens. The spinel types 
have additionally been represented as percentage of end 
member molecules. 
Spinels 

The most abundant cations are iron and titanium; with 
essential aluminum and manganese in subordinate amounts. 
Magnesium and zinc are usually present in minor to trace 
amounts. Chromium and vanadium have mutually sympathetic 
presence and distribution and are limited to oxides from 
basaltic samples. Trace amounts of calcium, nickel and 
cobalt may be present but show no sensible variation. In 
ulvospinels from salic lavas, chlorine and sulfur often 


appear in trace amounts. 
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Figure 6-4 shows the distribution of oxide analyses 
(both spinel and ilmenite type), in the system 
Ti0,-Fe0O-0.5Fe,0,. Four points do not fall close to the two 
equilibrium compositional joins, magnetite-ulvospinel and 
hematite-ilmenite. The two points below the 
magnetite-ulvospinel join are the PAO spinel samples 
previously discussed. The single point between the joins 
could be an unexsolved composition. In flow PAC a large 
magnetite phenocryst showed broad lamellae of unusual 
composition. These unusual lamellae plot on the FeO-TiO, 
join. Although it is dominantly a Ti-Fe composition, it is 
not a common mineral. On the basis of five oxygens there are 
2.98 cations; giving (Fe,M)Ti,0O; as a possible formula, 
indicating the low oxidation end member of the 
pseudobrookite series. It 1s relatively common to have 
exsolution relations between magnetite and ilmenite in the 
Level Mountain lavas. This occurrence of a pseudobrookite 
could indicate a peculiar activity for titanium in this 
flow. Another possibility is that these lamellae, with a 
lower density and higher titanium content than ilmenite 
formed at a different (lower) pressure than the 
crystallization of most lavas. 

Ilmenites 

Ilmenites, while not universally present, have been 
found in every lava type at Level Mountain. Essential 
cations are iron, titanium and manganese. Trace and minor 


constituents include zirconium aluminum, and magnesium. The 
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Figure 6-4, 
Level Mountain lavas and iron-titanium oxides Ti‘’-Fe?*-Fe?* 
cation proportions. 
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occurrence of traces of zirconium is most typical for 
ilmenites in the salic lavas, rather than in the basalts; so 
it would seem that the zirconium content in ilmenite 
reflects its content in the rock. 

In salic lavas the ilmenites contain detectable niobium 
which was confirmed by wavelength dispersive analyses. It is 
thought that zirconium and niobium probably substitute for 
titanium as a perovskite type component. The levels of these 
minor elements vary from rock to rock, but seem to be 
present in small concentrations in several coexisting 
mineral phases. 

The importance of these oxide minerals is twofold. 
Depending on the total iron-bearing mineral assemblage, an 
oxide phase can be used as an indication of oxygen fugacity, 
(Buddington and Lindsley, 1964; Carmichael and Nicholls, 
1967). The prevalence of ilmenite, or the absence of oxides 
in the case of the peralkaline trachytes, indicates fairly 
low £0,. The distribution of whole rock compositions with 
respect to Fe?*/(Fe** + Fe*®*) is greater than for their 
oxides. Salic rocks have greater variation in oxidation 
State of iron and are generally lower in titanium than the 
basic or intermediate rocks. While the fractionation of 
magnetite and ilmenite may drive the iron oxidation ratio up 
Or down, it almost invariably depletes the residual melt 
with respect to titanium. The basalts may have two exsolved 
Oxide phases, while the salics tend more commonly to have 


one or none, as is common for the peralkaline compositions. 
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Basalts are usually driven to higher Fe** (lower f0,) by 
oxide fractionation. The basalt PAO has three very different 
oxides. The ilmenite composition can be considered to be an 
ordinary phase for a rock of this composition, but the other 
two cannot. Petrographically these are xenocrysts and they 


are high in spinel content. 


ALKALI FERROMAGNESIAN MINERALS 

Aenigmatite 

Aenigmatite, while an uncommon mineral in igneous rocks 
as a whole, is a characteristic phase in peralkaline lavas 
and their plutonic equivalents. Aenigmatite is a triclinic 
inosilicate, with a deep reddish-brown colour, anda 
composition that does not vary appreciably from its ideal 
formula: Na,Fe,TiSi,.O,, (Deer, Howie and Zussman, 1978). Its 
occurrence, previously reported on Level Mountain (Ostensoe, 
1960), confirms the peralkaline condition and indicates a 
low oxygen fugacity (Marsh, 1975). The best examples of this 
phase on Level Mountain occur on a whaleback outcrop north 
of Meszah Peak. There, striated blades of red-bronze 
coloured aeigmatite up to 1.2cm long were found in 
miarolitic cavities of a trachyte stock, see plate 6-1. An 
XRD pattern was obtained for this material with the 
calculated d spacings presented in Table 6-7, showing good 
agreement with published values (Thompson and Chisholm, 
1969; Kelsey and McKie, 1964; Ernst, 1962). Although 


aenigmatite is common in the peralkaline lavas from Level 
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Plate 6-1. 
S. E. M. photograph of aenigmatite prisms in a miarolitic 
cavity. Blocky pseudohexagonal crystals are sodic pyroxene, 
alkali feldspars are tabular. Note 400 micron scale bar 
below sample number. Operating conditions for Cambridge 150 
Stereoscan S.E.M.: 20KV EHT, 3 amp filament current at 
Saturation level, 0.15 picoamp beam current, prefocussed 
condenser setting of 2.5 amps at 90 volts, focussing 
condenser setting of 3.0 amps at 90 volts, working distance 
12mm. 
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Table 6-7 xRDL Pattern for Aenigmatite 8/25-52/5880 


d(A) I/1j hk] 
8 .088 35 001,010 
7.400 1 01T 
6.366 1 171 
4.820 4 011,11T 
4.412 4 072,027 
4.201 3 221,207 
3.708 7 022 
3.487 2 230,231 
3.324 1 11 
3.24] 4 130, 102 
3.153 100 012,021 
3.062 1 122,321 
2.945 1 013,03T 
2.861 1 331,103 
2.822 1 122 
2.757 4 322, 123 
2.710 18 731,003 
2.657 2 1T3 300 
2.549 17 213,420 
2.460 3 142,202 
2.416 7 022,233 
2.346 3 170,470 
2.307 2 301,421 
2.220 1 323,273 
2.199 3 132,042 
2.120 13 207,333 
2.073 3 252, 350,043,034 
2.009 4 122,250 
1.945 2 221,232 
1.913 2 032,131 
1.814 2 

1.732 2 

1.682 2 

1.635 4 

1.625 8 

1.613 3 

1.590 

1.575 1 

1.559 2 

1.555 2 


1 Pattern made using CuKa radiation, from 2° to 60° 20, with 
slit settings of 0.5°-1.0°-0.5°, with an 8 second time 
constant at a scan rate of 0.5° 20 per minute and a chart 
speed of 30 cm per hour. 
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Mountain, in the polished thin sections prepared for 
microprobe analysis aenigmatite was confined to the 
groundmass. In those samples, aenigmatite occurred in 
patches up to 20 microns in size, but texturally it was 
intergrown with other phases making it impossible to obtain 
good quality microprobe analyses. Energy dispersive spectra 
obtained with a focussed beam and twenty second counting 
times contained peaks for Na, Si, Ti, Mn, and Fe. Sometimes 
a low intensity Zn peak was noted but peaks for Mg, K, and 
Ca were absent, although small amounts of these elements are 
typically reported in the literature. For comparison, 
aenigmatite compositions from other peralkaline lavas are 


presented in table 6-8. 


ACCESSORY AND ALTERATION MINERALS 

Probe analyses of selected accessory and alteration 
phases are presented in Table 6-9. Apatite has been 
recognized as an accessory mineral in all lava types at 
Level Mountain. Usually apatite occurs as acicular 
inclusions in the primary phenocrysts. In the intermediate 
lavas, particularily benmoreites and phonolites however, 
apatite occurs in phenocryst aggregates along with 
hortonolite, andesine, and ulvospinel. Such an analysis from 
flow PBP is given in the table. Chlorine is present in the 
apatites particularily from the salic lavas. Flourine may be 
present as well but detection of the soft XRays from this 


element is beyond the analytical capability of the EDA 
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Table 6-8 Aenigmatites From Other Alkaline Volcanic Centers 
Structural Formulae aS Cations per 20 Oxygens 


Si 
Al 
Fe3* 


Al 
Fe3t 
Ti 

M 
Fe2+ 
Mn 
Ca 


Ca 
Na 
K 


100 Ti 
Fe + Ti 
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Key to Sample Numbers 
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Pantelleria obsidian (Carmichael, 1962) 
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m 
nm 


Nandewar, New South Wales, peralkaline trachyte (Abbott, 1967) 
Rainbow Range, B.C., pantellerite (Bevier, 1978) 
Mt. Edziza, B.C., comendite (Yagi and Souther, 1974) 
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Table 6-9 Microprobe Analyses of Accessory Minerals from Level Mountain 
PBP Chliorapatite 


Si 
Mg 
Fe 
Mn 
Ca 
Na 
P 

Cl 
0 


8/27-63/6788 Zircon 


Element 


Wt .% 


0.24 
0.16 
1.02 
0.16 
39.93 
0.06 
18.57 
0.08 
39.75 


Element 


Wt .% 


15.80 
49.44 
0.06 
0.28 
0.07 
34.37 


PAC Septechlorite 


Si 
Ti 
Al 
V 

Fe 


Mn 
Mg 
Na 
0 


Element 


Wt .% 


10.98 
0.04 
ice 
0.05 

29.28 


0.17 
Toon 
0.22 
49.79 


Structural 
Formula 


0.04 

0.032 

0.090 

0.014 0 
4.87 

0.01 

2.93 

0.01 

13.00 


Structural 
_Formula 


1.01 

0.976 

0.002 0 
0.009 

0.00 

4.000 


Structural 
Formula _ 
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-681 0 
.015 

33 
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-459 12.0 
.744 

147 
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Theoretical Formula clinochlore/cronstedite 
, Fest) 2(Fe3*,Al,V,Ti Si )g029(0H)16 


(Mg,Mn, Feet 


Recast 
Oxide Analysis 

Si02 25.98 
TiO 0.08 
We ke 
2 
Foob3 20.36 
Fed ra a 
MnO 0.23 
MgO 13.88 
Nao0 0.33 
H20 Tsc4) 
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system. The recognition of phosphorous and halogen bearing 
phases is important to understanding the minor element 
geochemistry and gas phase composition of the salic lavas. 

Carbonate aS a veSicle and fracture infilling is one of 
the most common alteration minerals. Carbonates from the 
basalts are usually calcite or aragonite and sometimes have 
noticeable levels of Sr as verified by atomic absorption 
analyses. The groundmass/intersertal carbonate from 
phonolites, tristanites and mixed flows is consistently high 
in Fe and Mn. In a molecular sense these are 85% calcite, 
10% siderite, 5% rhodochrosite. Some of these high Fe-Mn 
calcites have inclusions of pyrite and are probably 
hydrothermal in origin. 

Zirconium minerals are usually found as accessories in 
the peralkaline salic lavas. Baddeleyite microphenocrysts 
have been found in comendite pitchstones. Zircon is the most 
common Zr bearing mineral found as inclusions in pyroxenes 
and alkali feldspars as well as in the mesostasis. Butterman 
and Foster (1967) studied the pure ZrO, - $10, system and 
found Baddeleyite to be stable below 1170°C which is at a 
lower temperature than the stability field for zircon. This 
would seem to indicate that the liquidus temperature for the 
baddeleyite bearing glassy comendites is lower than that for 
the other peralkaline flows which more commonly contain 
zircon. 

The most abundant alteration mineral in basalts was 


petrographically chlorite, but probe analyses showed it to 
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be essentially iron silicate with minor Al, and hence more 


akin to septechlorite or septechamosite. 
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CHAPTER 7. PETROCHEMISTRY 


INTRODUCTION TO THE WHOLE ROCK CHEMISTRY AND ANALYTICAL 
METHODS 

More than one hundred and ten major element analyses 
were performed by electron microprobe EDA of whole rock 
glasses after the method of Schimann and Smith (1980). This 
type of whole rock analysis has previously been reported by 
Rucklidge et al (1970) and by Reed and Ware (1975). Data 
reduction was by the method of Smith and Gold (1979) and 
Smith (1976). Thirty XRF analyses, mostly on sample splits, 
were performed by Dr. J.G. Holland at the University of 
Durham after the method of Brown et al (1977). The XRF data 
additionally include minor and trace element concentrations. 
Other trace element data were collected by atomic absorption 
spectroscopy. Whole rock isotopic compositions for oxygen, 
strontium and lead were surveyed for the entire range of 
petrochemical variation. 

All rock powders intended for major element analyses 
were measured for H,0O- and loss on ignition. Total water 
analyses were made by the Penfield tube method. All whole 
rock powders were tested for the presence of carbonate. 
Those few samples, mostly older basalts, which tested 
positive were accurately measured for CO, in an 
evolution-absorption apparatus modified after Goldich et al 
(1959). The weight percent of ferrous iron was determined by 


acid dissolution of a known weight of sample and titration 
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with a standardized solution of KMn0,. The sample weights 
were corrected for water content before calculating the 
weight percent of Fe?*. 

To prepare glasses for microprobe analysis, the whole 
rock powders after L.O.I. were transferred from crucibles to 
molybdenum pedestals and fused under reduced pressure in an 
image furnace. This apparatus was constructed by K. Schimann 
(Schimann and Smith, 1980). Probe analyses were made against 
anhydrous mineral standards with concentrations of standard 
elements in excess of sample concentrations. Minor elements 
were referenced to default values. The standards used, with 
their weight concentration of critical elements, are 
presented in table 7-1. Whole rock analyses were recast as 
oxides with iron calculated Sample homogeneity was verified 
by 50 second scans on different areas. as ferrous oxide. 
Those analyses with totals greater than 98% were considered 
acceptable. Most analyses had totals greater than 99%. The 
errors associated with these whole rock analyses can be 
limited to three kinds. The first, acciracy,, is related to 
concentration. The lower the concentration of an element in 
the sample the poorer the counting statistics. For instance 
Si has three to four significant figures while Mn has only 
one to two. In the case of the energy dispersive analyses of 
glasses, the concentration of each element in the standard 
is greater than in any sample. Some replicate analyses 
(E.D.A.), and XRF analyses for comparison are presented in 


table 7-2. The variations between these analyses should be 
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taken as representative of errors for all whole rock 


analyses presented. 


MAJOR ELEMENT GEOCHEMISTRY 


Introduction 

Eleven series of histograms have been prepared for the 
prominent lava types, figure 7-1. Along the abscissa are 
plotted the weight ranges for the eleven major chemical 
components. The number of samples in each range is plotted 
along the ordinate. Components (except H,0*) have been 
recast aS oxide percent on a dry basis. Total iron is 
expressed as ferrous oxide. Basaltic and intermediate lava 
types have been classified according to Irvine and Baragar 
(1971). The distribution of Level Mountain lavas on Irvine 
and Baragar's normative classification plot follows the 
typical sodic alkali basalt series, see figure 7-2. 
Peralkaline salic lavas have been classified according to 
MacDonald (1947b). Source derivation tests have been 
performed for all Level Mountain lavas. Partial melting 
calculations for basalts used a pyrolite source (Green and 
Ringwood, 1967) and for salics used a clinopyroxene 
amphibolite from Three Valley gap. Results of these 
calculations appear in the tables of chemical analyses as 
*#MELT. 

For ready comparison of the overall chemical variation 


of Level Mountain lavas with other alkaline volcanic suites 
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Figure 7-1, 
Histograms of chemical composition for the various Level 
Mountain lavas. Number of samples versus oxide content (dry 
basis). Rock types: RHY - rhyolite; CoM - comendite; PANT - 
pantellerite; TRACH - trachyte; BEN - benmoreite; TRIS - 
tristanite; PHO - phonolite; MUG - mugearite, HAW - 
hawaliite; AB - alkali basalt; ANK - ankaramite. 
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Figure 7-2. 
Level Mountain lavas on the normative classification plot of 
Irvine and Baragar (1971). 
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the familiar AFM and CaO - Na,O - K,0 plots, both mol ratio, 
are presented in figures 7-3 and 7-4 respectively. The dense 
clusters of points reflect the bimodal distribution of lava 
types with few intermediate compositions. The bimodal 
distribution and compositional gap is also well expressed in 
the atomic ratio plot (Fe+Mg)/Al versus (Na+K)/Al, figure 
7-5. In this plot the intermediate lavas, which are thought 
to be formed by a combination of magma mixing and crustal 
contamination, are represented by squares, while the primary 
magma types are labeled with dots for salic and triangles 
for basalts. The next pair of figures 7-6 are atomic ratio 
plots of Si/Al versus agpaitic index. Note the similarity in 
overall range and distribution for Level Mountain lavas as 
compared with other major Cenozoic volcanic centres of the 
Intermontane Belt. 

Ankaramites (Olivine Rich Basalts) 

In the case of the ankaramites, which tended to be 
coarsely crystalline and altered, only four analyses were 
made, table 7-3. For the norm calculations on these oxidized 
samples, iron oxidation state has been adjusted according to 
the TiO, content (Irvine and Baragar, 1971). These lavas are 
higher in magnesium than the other Level Mountain basalt 
types, while the concentration of all other major elements 
considered is similar to, or lower than, in the basalts. 
Table 7-4 presents some representative ankaramite analyses 
from other localities. The only published ankaramite 


analyses from the Cenozoic of British Columbia are from 


m doe. eehlyO.s - 


.Vlewvitoagees b-* 


Mee: 
7), ao oe tna? SA 2 


got eerupet fn bg 


tyudéitdatbh Lebomrd eri! {osktea eaiiiog 9 
: 7 
<td of nortf2ogaoes a2 siecaontla = t 4 4 
. ey 
2 {few oeis @: gee ae week 
ats 
[A\ (44eH) evetey [AN (phteen) tolg otand 3gmod, 
' a ie 
w .weval sasibeqmetng ada 401g eins al 
as 
b tin bem 70 anidmas. s a 
oe | 
a 5 a 
: bs inezs7get 935 sn a a 
lee 103 era dziwv beledsl sis aX3 7 
s-f e#o400 }o %i£sq Jxon dT ae — 
¥ 
seal e906 auUatsv iA 2 Ae : 
; ee 
; yi noigudlizeib bas. senaa ff 
: 0232 10f su tento Pr " *’ 
- ‘ = “us J 
ere 
l\e¢tigesd dol® Cat 0} eet: 
7 ai Yaa 
é isidw .ebtimevedne sdd tose: rs 
isos t002 yine ,barstle bas enilis re 
a 4 , a? 
“ no shel jeiuglen Won nd 209. 


168 


Pilourne 7-3. 
AFM plot for Level Mountain lavas. 


Figure 7-4. 
Distribution of Level Mountain lavas in CaO-Na,O-K,0. 
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Figure 7-5. 
Atomic ratio variation diagram for Level Mountain lavas. 
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Figure 7-6. 
Si/Al versus agpaitic index atomic ratio plots for Level 
Mountain lavas and other Late Cenozoic lavas of the 
Intermontane Belt. There are really two distinct trends 
here. For the basalts, AI variation is independent of Si/Al, 
while for salic lavas, Si/Al variation is largely 
independent of AI. 
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Table 7-4. Ankaramite Chemical Analyses From Other Localities 


Spanish 
Weight Creek New New 
Percent QV19&QV20! Georgia’ Hebrides’? Tahiti* Hawaii® Alaska*® 
$i02 45.87 48.68 48.59 43.26 Al 325 49.2 
Ti02 2.2 0.43 0.67 3.40 1.61 One 
A1203 '3g44 12.44 St 9.69 2.07 2.4 
Cr203 0.08 
Fe203 Me ke 3.90 4.81 3.66 1245 2.0 
Fed 9.10 6.49 4..97 8.97 10.41 6.8 
MnO 0.16 OF2Z0 0.19 0.16 Ose Ou 
MgO TOI 11.36 5 he Ge Ie 12.64 19.96 1954 
CaO 10.60 dod 7, 1 Sit 12°10 7.88 18.9 
Na20 3.34 1.93 2216 1.59 1.38 Ose 
H20 1258, 1 ee ia 0.95 1.18 0,35 O.1 
P205 0.46 0.24 0.14 0.61 One) — 
H20 0735 1.68 0.45 2.46 Oo i2 0.8 
TOTAL 100.03 WOO 100.33 99.72 99.82 99.9 


(1) D.W. Fiesinger and J. Nicholls (1975). 
(2) R.L. Stanton and J.D. Bell (1969). 

(3) A.J. Warden (1970). 

(4) A.R. McBirney and K. Aoki (1968). 

(5) G.A. MacDonald (1949). 


(6) J.C. Rudnick and J.A. Noble (1959). 
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Spanish Creek in Well's Gray Park (Fiesinger and Nicholls, 
1977). The Spanish Creek flow is higher in Ti, Na, K and P 
and lower in Mg than the Level Mountain ankaramites. 

At Level Mountain, the older ankaramites differ from 
the other basalt types mainly in amount of phenocryst and 
possibly xenocryst olivine, or in degree of fractional 
crystallization of olivine. For the variation index, 
(CaOt+Na,0+K,0)/Al1,0,; mol, which is a measure of alkalinity, 
the high Al,0O, ankaramites are comparable to the alkali 
basalts at 1.3. The late ankaramites have very high contents 
of clinopyroxene (60%) both in thin section and in their 
norms. Their values of the aforementioned alkalinity index 
at 4.9 and 8.3 are out of character for normal basalts and 
even the strongly alkaline basanitoids which are generally 
less than two. They are spatially and temporally associated 
with peralkaline trachytes that sometimes are laden with 
nodules rich in clinopyroxene. It is interesting to note 
that fractionating a feldspar of any composition from these 
ankaramites would produce peralkaline residua. The levels of 
alumina in these ankaramites (5.67%), are very low for 
alumina in any igneous rock type. Their peculiar chemistry 
is somewhat akin to komatiites which are unknown from the 
Tertiary and not associated with alkaline suites. The 
peculiar circumstances of their occurence at Level Mountain 
in conjunction with their refractory chemical nature suggest 


that the low alumina ankaramites are evolved rocks. 


b 
i 
‘ledot™ a6 re i262 cao afi a 


cull a4 —_ , 4 
+87) 4 « ) ae Serio is es wol? opt) a ah, 


iejineararna noe2nuesM lave 307 nail 
.2 S83 1G THATS r9bf0..2a?t + ane 


reheat Fo Inyons A ¢imiem Py 


1hIROD $132 ao TMBIBANS 


ngtaains 936f° ont 3% 
ate 7 ay 
page ptiat? a # dod | (toa) seein ‘9 


ag 
[saee Lemon + $328 2vado eo sug et 
io s16 doldw abtoetiasasd on niledse’ io 


Vw re 


siege: “935 et 


14.4 


= 
I ° 
wos WOl 23 
2 ime? i6i iW 
— on 2 _ 
pil. Ta a0 
7 
z -@e2¢U 


niasnioM Aeve G3 
epee eivran Lee 


—— 


174 


Alkali Basalts from Level Mountain and a Comparison with 
B.C. and World Averages. 

Chemical analyses of twenty-four primitive basalts from 
Level Mountain are presented in Table 7-5 along with their 
C.I1.P.W. norms. These include both alkali basalts in the 
sense of Irvine and Baragar (1971) and associated lavas 
which are petrographically and chemically similar but 
contain normative hypersthene. Their order of presentation 
in the table is according to stratigraphic succession with 
examples from the plateau units 1, 2, 3, and 4 and also the 
Stratocone units 5b, 6b or younger. Chemically primitive 
basalts are present for most of the eruptive history. For 
the basalts, the only major oxide with a bimodal 
distribution is MgO. Among the minor components, BaO and H,O 
also showed bimodal distributions. Analyses with high and 
low values for MgO were grouped and averaged. The higher MgO 
group was found to contain significantly less Al,0, and Na,0O 
than the low MgO group, but comparable levels of all other 
components. This division on MgO did not correspond to a 
clear division in time or space for occurrence, as the low 
MgO group includes both analyses from the basal plateau unit 
and from the stratocone cap. The high MgO group could be due 
to less olivine fractionation (olivine enrichment) or a 
higher degree of partial melting. Slight olivine enrichment 
is petrographically supported. Similar cluster analyses were 
performed on simple tests of the mean for BaO. The high BaO 


group, again not stratigraphically distinctive, contained 
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marginally higher values of P,0O; and K,O and marginally 
lower values of MnO. Perhaps the high BaO+P,0,+K,0 lavas 
represent smaller degrees of partial melting or residual 
element type enrichment due to fractionation of olivine and 
clinopyroxene. For the H,0 division, the high total water 
group showed greater levels of $10, and MgO with lower 
levels of Na,O. Replacement of 5 to 10% of the groundmass 
glass by a septechlorite could account for the H,0 and MgO 
increase but this would cause an FeO* increase, which is not 
observed, and an $10, decrease rather than the observed 
increase. Also the high H,0O group does not uniformly show 
more extensive groundmass alteration than the low H,0 group. 
Another possibility is that the high H,O group represents 
Slightly more hydrous melting conditions than prevailed for 
the majority of the basalts. For any of the three 
clusterings discussed, it should be emphasized that no 
cluster corresponds to a distinctive group of flows either 
in age or location. If all of the Level Mountain basalts are 
primary, perhaps these slight chemical differences represent 
local variations in source composition or differences in 
melt extraction and derivation. These clustering tests are 
Summarized in Table 7-6. 

Level Mountain basalts are petrographically all of 
alkaline affinity, but chemically some approach the 
Subalkaline type considered to be transitional to 
tholeiites. Some of these lavas in fact have hypersthene as 
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total absense of phenocryst orthopyroxene or groundmass 
subcalcic pyroxenes. Some of the basalts have high enough 
Silica contents for dominant hypersthene, (silica 
Saturation) yet have sufficient alkalis to plot above 
MacDonald and Katsura's (1964) alkaline/tholeiitic line. The 
average alkalinity index, (mol(CaO+Na,0+K,0)/Al1,03), of 1.42 
1s not appreciably different than a typical value for Level 
Mountain or B.C. MgO is high enough and TiO, low enough that 
these rocks do not belong with hawaiites either. The 
problem, given their petrographic characteristics and 
alkaline association on one hand versus their silica 
Saturation seen aS normative hypersthene on the other are 
(i) how to classify them and (ii) are they really 
distinctive chemically from more ordinary alkali basalts. 
The Level Mountain basalts were divided on normative 
criteria into alkaline and transitional types to see what 
the systematic differences were, if any. Additionally, the 
transitional types from Level Mountain were chemically 
compared to similar basalts from central B.C. Again tests of 
mean and standard deviation were applied. The results are 
presented in table 7-7. The transitional basalts from Level 
Mountain and Central B.C. are virtually indentical. Any 
differences in major element chemistry for the bona fide 
alkaline basalts at Level Mountain versus the transitional 
types are within one standard deviation. It is felt that the 
calculation of dominant normative hypersthene for some 


basalts is not statistically significant of any major 
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chemical distinction either in classification or genesis. It 
may relate in part to the oxidation state of the iron in the 
analysis and its effect on calculated normative iron bearing 
Silicates. The alkaline basalt group does have marginally 
higher levels of TiO,, K,0O and P,0O, with lower SiO, than the 
transitional basalts and a slightly higher alkalinity index. 
This may reflect a difference in the degree of partial 
melting, with the transitional type representative of a 
higher percentage of partial melting or melting at higher 
pH,O. Again it should be emphasized that the two types are 
thoroughly interspersed stratigraphically. Even if the two 
types represent differences in mantle melting conditions or 
extent, the chemical variation between them is continous 
rather than discreet. Any variations in the derivation 
process are thus probably continous as well. 

Table 7-8 presents average and representative alkali 
basalt analyses for the Cenozoic of B.C. The lavas and 
averages presented here were selected to represent alkali 
basalts of similar derivation and at a similar stage of 
evolution. Hawaiites have been excluded as these are 
probably derived from alkali basalts via extensive 
fractional crystallization. A range of 14.0<(CaO+MgO)<20.0 
weight percent was used to restrict the effect of 
fractionation processes. A cutoff of 
1.3<(CaO+Na,O+K,0)/A1,0; mol ratio was applied to select 
lavas of undeniable alkaline affinity. Thirty four lavas 


were found to meet these criteria and were hierarchically 
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averaged by area, for a data set of eleven averages. These 
were in turn used to generate a B.C. average and also as a 
test for source chemistry variation and differences in 
degree of partial melting. In a similar fashion a world 
average alkali basalt, exclusive of B.C., was generated from 
140 analyses selected from the published literature 
representing fifteen sub-regions. These are presented in 
Table 7-9. For comparison, this world average is also 
presented in Table 7-8. The Level Mountain average and B.C. 
average are within one standard deviation for all 
components. The B.C. average and world average are within 
one standard deviation for all components. Level Mountain 
appears to be higher in $10, than the world average but if 
more stringent selection criteria are used, the Level 
Mountain average of fourteen alkali basalts from table 7-7 
is within one standard deviation of the world average for 
all oxides. On this basis all of the B.C. basalts are judged 
to be ordinary alkali basalts as they fall within the 
variation of continental and oceanic alkali basalts from 
elsewhere in the world. There is no obvious strong 
correlation between the B.C. analyses and alkali basalts 
from any specific tectonic setting or region of the world. 
Neither the basalts from Level Mountain nor from B.C. as a 
whole, are obviously "rift type" basalts or "continental 
Margin type". 

It is noteworthy that all of the B.C. plateau types 


have above average values for FeO* as do all of the lavas 
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Table 7-9. Comparative analyses of Alkali Basalts. 


Recast *Var World Std Hawaii New Zealand OFB NSW Uganda 
World Avg Avg 140 dev (3) (14) (13) (75) 
$102 47.18 3.06 47.06 +1.44 46.38 49.14 44.80 49.56 47.51 45.2 
T102 cy 22.88 2.36 = .64 2.79 3.17 1.96 1.42 2.82 2.3 
A1203 14.97 8.51 14.93 £1.27 14.74 15.71 13.86 16.09 16.85 16.0 
Feo* 11.48 11.79 11.45 21.35 10.95 13.45 WZ225 2 AOA 10.53 13.95 
MnO cae 11.76 ie / = 702 213 0.15 O.47 == [om he} 0517 
MgO 8.36 2A IO 8.34 +1.76 9.81 6.03 11.07 7.69 5.62 7.60 
Cao 9.96 10.98 9.93 +1.09 10.54 8.38 10.16 11.34 9.28 9.30 
Na20 3.10 11.65 3.09 + .36 2.70 3.69 3.19 2.80 3.47 3.20 
K20 1.08 50.00 1.08 + .54 0.85 7.39 1.09 0.24 1.08 0.81 
P205 0.58 37.93 58 Ee HS 0.37 0.73 Oo. 55 dad 0.93 0.39 
H20 0.75 54.67 On7o =. .41 0.33 ss 0.73 a tes 0.59 1.00 
100.00 99.74 100.18 99.84 99.83 100.00 99.83 939.95 
Pantelleria —. Africa Kurile Boseti Anjovan Y &T Baikal SRP Massif Centra) 
(3) (2) (2) 1 1962 Rift (15) (6) 
(2) 
$102 47.32 48.10 48.12 47.92 46.88 46.53 47.37 46.34 44.72 
Ti02 3.16 1.45 2.07 1.86 2.39 2.28 2.40 be 2.63 
A1203 15.57 13.20 15.31 16.37 12.08 14.31 15.00 14.35 14.45 
Feo* 11.85 10.31 9.36 10.00 12.33 12,65. 11.49 13.73 10.66 
MnO 0.15 0.17 0.18 0.18 0.20 0.18 mh .19 -20 
MgO 5.81 10.13 9.32 7.04 98.50 9 54 7.78 7.81 10.41 
CaO 9.53 12.06 8.86 11.24 10.29 10.32 8.17 9.52 9.90 
Na20 3.64 2.45 3.02 2.92 2.92 2.85 3.13 2.84 3.83 
K20 1.31 0.35 2.30 .66 1.01 0.84 1.67 .82 is ta 
P205 0.84 0.24 0.48 .36 0.48 "<= 48 .80 0.83 
H20 0.81 1.04 nh 1.04 0.54 0.08 1.54 34 1.01 
99.99 99.50 99.14 99.59 98.62 99.58 99.20 99.51 100.06 
1. Hawaii alkali basalts (3) avg Polulu, Haulalai from MacDonald 1949, Mauna Kea from 
MacDonald to Katgura (1964) 
2. New Zealand basalts (14) avg Banks Peninsula Speight 1938 
Auckland (13) (1974) 
3. OFB (75) Ocean Floor Basalts al! MORB 42 Atlantic + (Carlsborg, E Pacific, JdeFuca) 
remainder Pearce (1976) 
4. NSW Nanderwar Abbot (1969) 
5. Uganda Moroto Varne (1968) 
6. Pantelleria (3) Villari (1975) + Romano (1968) 
7. €. Africa Erta’Ale Range (1975) Barberi et al 
8. #1026 Oliutorskii Progi (Koryakskoe Nagorie) ed NB Zolomova (1966) 
9. Boseti (2) 152f 153 avg Brotzu et al (1975) Ethiopia 
10. Anjovan East Indian Ocean Thompson + Flower 1971 
11. Yoder + Tilley (1962). Martin (1975) 
12. Koselev (1977) avg 118AB Baikal Rift 
13. Snake River Plain (15) Stout + Nicholls 1977 
14. Massif Central Paul (1970) nodule bearing basalts (basanites) 
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from the Stikine. This could imply that the mantle source 
region beneath the Stikine (Heart Peaks, Level Mountain, 
Stikine River, Edziza) is enriched in FeO*. This could also 
imply more extensive melting from iron enriched upper 
mantle. If alkali basalt is a primary mantle derived magma, 
and if these lavas in question are related only in type but 
not by genetic sequence, then these chemical variations 
might represent vertical stratification or variation in the 
source region. High Ti, K, and P probably represent 
undepleted mantle as does high Na and Al versus Mg. 
Variation diagrams have been constructed for SiO, 
versus TiO, and Al,0O,;, versus MgO for the B.C. data set, 
figure 7-7. The effect of variable source contribution for 
typical spinel lherzolite constituents is shown as vector 
Subtraction from the B.C. average. Model lherzolite phases 
were forsterite, hypersthene, chrome diopside and spinel 
with a whole rock Mg/(Mg+Fe) atomic ratio between 0.93 and 
0.94. This composition was chosen to approximate the B.C. 
nodule compositions presented in a subsequent section. Kuno 
and Aoki (1970) estimate the bulk Mg/Mg+Fe of the upper 
mantle to be close to 0.89 with olivine equal to Fo,, and 
spinel less than or equal to 0.8. Using these slightly less 
refractory values only affects the spinel vector's direction 
by about 3° without affecting its length. The other vectors 
are unaffected. The distribution of points in both plots is 
Parallel with the clinopyroxene vector. The length of the 


field on the MgO versus Al,0, plot is consistent with either 
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Figure 7-7, 
Oxide variation diagrams for basic lavas of the IMB. Vectors 
demonstrate the magnitude of fractionation effects for 
dispersing the compositions. The point "W" is the world 
average alkali basalt. 
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(i) a difference of 22% in the contribution of clinopyroxene 
to the partial melt or (ii) between 13% and 20% difference 
in the amount of partial melting responsible for the 
variation in compositions. If case two is true the major 
difference in percent of partial melting would still likely 
be related to the contribution of clinopyroxene. Model 
calculations for percent partial melting to derive these 
basalts from pyrolite range from 10% to 20% which agrees 
with the figures above. The variation in the TiO, versus 
Si0, plot is also parallel to the clinopyroxene vector but 
here spinel has much greater leverage such that the entire 
range in SiO, values could be explained within 6% variation 
of spinel contribution to the partial melt. 

Variation diagrams of P,0, versus T10,, Na.,0O and K,0 
have also been constructed for the B.C. data set, shown in 
figure 7-8.These plots are particularily useful in 
addressing the origin of the compositional variation. The 
least squares fit straight lines represent residual element 
models. The slopes (ratios) for K,0 and TiO, versus P,0, are 
very close to unity while Na,0/P,0;, is about 2. If, for 
instance, P,0, and K,0 do not enter into any crystallizing 
phase and their ratio is fixed in the source region, then 
fractional crystallization of any non P,0;, and K,0 bearing 
phase will result in increased levels for both oxides, but 
Still in fixed proportion. The vector corresponding to ten 
percent less clinopyroxene concentration is again 


Subparallel to the trends but it is immediately obvious that 
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Figure 7-8. 
Oxide variation diagrams for basic lavas of the IMB. 
Straight lines are least squares linear regression fits to 
the data. Curved lines are least Squares exponential fits to 
B. C. data plus pyrolite. The point "W" is the world 
average. The point "B" is the B. C. alkali basalt average 
with the end of the arrow indicating the B.C. hawaiite 
average. The chemical dispersion here is probably a source 
Variation and fractional fusion effect rather than 
differentiation. 
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fractionation alone is inadequate to explain the variation. 
To explain the variation in P,0, some 80% fractional 
crystallization would have to occur between compositions 
like Edziza's and those at Heart Peaks or Atlin. The 
variation in P,0O, could be explained in terms of source 
content or contribution of apatite. Then Edziza's partial 
melts would receive about 0.5% apatite while Heart Peaks 
would receive about 1.5% apatite. It is possible that all of 
the B.C. lavas represent the same degree of partial melting 
but that available apatite and other trace minerals are 
variable. 

Another way to explain minor element variation is 
through variable degrees of partial melting of spinel 
lherzolite with a rather ordinary incompatible or minor 
element content. For this type of model, the minerals 
containing the minor elements melt first and further degrees 
of melting dilute their concentration. Na.,O and TiO, are 
also contained in several mineral phases so they will behave 
like diluents rather than pure residuals. For this reason 
expressions of the form 

b 
P,O, = a(Na,0) 
were fit to the B.C. data set and the same hypothetical 
source pyrolite composition (Green and Ringwood, 1967) used 
in the basalt derivation calculations of ROCK. These models 
imply that the Heart Peaks and Atlin represent the lowest 


degree of partial melting and Edziza the highest. In 
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comparing the B.C. minor element distribution to the world 
average point, if we assume both have the same minor element 
source signatures then it can be seen that most of the B.C. 
lavas represent a higher degree of partial melting, implying 
higher heat flow or steeper geotherms for B.C. If however, 
we assume that all of the lavas represent the same degree of 
partial melting, then for their P,0, levels the B.C. lavas 
are enriched in alkalis. This could tie into the occurrence 
of related peralkaline salics in B.C. Another possibility is 
that if melting degree and alkali contents are typical, then 
the mantle source region for B.C. lavas is slightly depleted 
in phosphorous with respect to the world average. Perhaps 
some additional insight can be gained in examining element 
ratio plot,s, figure 7-9. On the Ca/Si versus (Fe+Mg)/Si 
plot the olivine orthopyroxene and spinel vectors parallel 
the (Fe+Mg)/Si axis, while clinopyroxene has a slope of one. 
The scatter of points lies between these two vectors and 
closer to the origin than the world average point. This 
implies either a higher degree of partial melting or greater 
source depletion for the B.C. lavas than alkali basalts in 
general. On the Ca/Si versus P/Si plot, figure 7-10, most of 
the B.C. lavas plot at lower levels on both axes than the 
world average. This could imply, for the same degree of 
partial melting, that the B.C. source is depleted in both 
clinopyroxene and apatite. Alternatively if B.C. mantle is 
normal it could imply higher degrees of melting. The 


interesting thing about the scatter is that it occurs along 
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Figure 7-9. 
Ca/Si versus (Fe+Mg)/Si atomic ratio plots for basic lavas 
of the IMB. Vectors indicate slopes for ideal mineral 
stoichiometry, arrows point to increased contribution, melt 
vector indicates direction of increased partial melting. 
Compared to the world average alkali basalt, the IMB 
compositions appear to have a reduced contribution of 
clinopyroxene and spinel either from source depletion or 
increased partial melting. 
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Figure 7-10, 
Ca/Si versus P/Si atomic ratio plot for basic lavas of the 
IMB, Compared to the world average alkali basalt the IMB 
compoSitions appear to have a reduced contribution of 
clinopyroxene and apatite, either due to Source depletion or 
increased partial melting. 
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Figure 7-17. 
K/Si versus (Fe+Mg)/Si atomic ratio plot for basic lavas of 
the IMB. Note the direction expected for increased partial 
melting and the scatter of IMB values in K/Si both above and 
below the world average. The implication is that B.C. 
basalts represent about the same level of partial melting as 
the world average but that they have considerable source 
heterogeneity. 
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a line of constant sum. This is also perpendicular to the 
melting vector. This line could imply a relatively constant 
degree of partial melting but some source heterogeneity of 
the form highest source clinopyroxene content at the lowest 
apatite content and vice versa. The values of (NatK)/Si 
scatter around the world average, refer to analyses in Table 
7-8. With the low Ca/Si values for B.C. this implies a lower 
contribution of clinopyroxene to the melts at relatively 
Similar source enrichment. On the K/Si versus (Fe+Mg)/Mg 
plot, figure 7-11, the B.C. lavas are shown to have scatter 
around the world K/Si value. This is also true for total 
alkalis. Most of the scatter is perpendicular to the melt 
vector and parallel to amphibole and phlogopite which are 
considered to be possible alkali source minerals (Best, 
1974; Mysen and Boettcher, 1975; Boettcher and O'Neil, 
1980). This plot can be interpreted to mean that if the 
range of partial melting for B.C. is on both sides of the 
world average, then the B.C. source is depleted in minerals 
like phlogopite or amphibole. Further discussion of this 
point is deferred to the section on lherzolites. If this 
plot is really showing variable source depletion for hydrous 
minerals and incompatible or minor elements, then there is 
as much mantle heterogeneity between the Rainbow range and 
the Itchas, which are less as 100km apart, than there is for 
the whole of B.C. Also, it would appear that if this plot 
can show variation in degree of partial melting, then Heart 


Peaks, Edziza and Level Mountain have fairly similar minor 
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mineral content in the source, but reflect as great a range 
in partial melting as is observed for all of B.C. From this 
data analysis it does not seem possible to simultaneously 
retain the concepts of a chemically and mineralogically 
homogeneous upper mantle, a constant degree of partial 
melting implying a narrow range of geothermal gradients and 
an undepleted upper mantle with high incompatible element 
Signature as is typically assumed for alkaline lavas. 

The results of this analysis indicate that most of the 
chemical variation for B.C. alkali basalts can be explained 
in terms of variable source content or contribution of 
clinopyroxene to the melts on the order of twenty-two 
percent. The contribution of clinopyroxene to the B.C. lavas 
is also lower than for the world average alkali basalt. The 
minor element variation is consistent with variable degrees 
of melting, with most of the B.C. basalts being formed by 
more extensive melting than the world average. This could be 
accommodated by either higher geothermal gradients, from a 
minimum of 21°/km by Fujii et al (1981) to 25°/km by Parrish 
(1981), or by a less refractory upper mantle. The latter is 
suspected due to the high iron and silica contents and low 
(Fe+Mg)/Si values. 

Hawaiites and Their Relationship to the Stratigraphic 
Variations in Major Element Chemistry of all Basic Lavas 

The hawaiites have been classified according to Irvine 

and Baragar (1971). They are distinguishable from the more 


primitive basalts and the mugearites on the basis of 
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normative plagioclase composition and normative colour 
index. The analyses are arranged in stratigraphic succession 
and presented in Table 7-10. The hawaiites can be chemically 
distinguished from the basalts on the basis of their higher 
Ti0,, Al.0,, FeO*, Na,O, K,0 and P,0, contents and lower MgO 
and CaO contents. 

Within the hawaiite group at Level Mountain there 
exists a bimodal distribution for: TiO,, Al,0;, MgO, CaO and 
P,0,. None of the high or low clusters for any of The test 
of the mean is presented for two of the clusters, Al.,0O; and 
P.O; in table 7-11. The low Al,0O; group are much more 
strongly alkaline. The Al,0, division not only corresponds 
to the Ti0,., MgO and CaO division previously noted but also 
levels for SiO,, FeO*, Na,O and K,0 which were not initially 
apparent in the normal distributions of these oxides. The 
high P,0, cluster also contained high TiO, and low MgO and 
CaO. This is interpreted to mean that concentration of P,O, 
and TiO, proceeds with fractionation of typical basalt 
phenocryst phases such as olivine, clinopyroxene and 
plagioclase (Anderson and Greenland, 1968). 

Representative hawaiite analyses for other localities 
of the Intermontane Belt are presented in table 7-12. The 
average Level Mountain hawaiite is not statistically 
different from the overall average for the Intermontane 
Belt. However there do appear to be differences from area to 


area. 
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Oxide Variation Plots for Basaltic Lavas and Fractionation 
Arguments 

On the basis of the correlation calculations, nine 
major Oxide variation plots were selected to project most of 
the significant chemical variation. By dividing the lavas 
into basalts and hawaiites, very few clear trends emerged. 
At this point in the data analysis it was found that 
Significant trends existed for time-map units which included 
both basalts and associated hawaiites. Also alkali basalts 
and hawaiites for the same time package were found to be 
chemically distinctive from subsequent and preceeding units. 
Consequently the variation plots were prepared for 
distinctive time and map units. 

Major element variation diagrams for plateau units 1 
and 2 (the basal two eruptive sequences of the plateau 
building stage) are presented in figure 7-12. The effects of 
fractional crystallization are shown for vectors for sample 
PAR which is an alkali basalt with mol ratio 
(CaO+Na,0+K,0)/Al1,0; = 1.72. The lengths of the 
fractionation vectors correspond to the degree of fractional 
crystallization and the directions depend almost entirely on 
the fractionated phase composition and very little on the 
whole rock. The fractionation vectors can be freely 
translated from point to point as long as their orientation 
and length are kept constant. This is particularily true of 
a vector like plagioclase on the FeO-MgO plot or olivine on 


the Na,0-CaO plot. For compositions with one or more 
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Figure 7-12, 
Oxide variation plots (anhydrous weight basis) for lower 
plateau lavas (units 1 and 2). Vectors denote possible 
effects of mineral fractionation. 
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principal components in the plane of the variation diagram, 
such as ilmenite on FeO#-TiO, or olivine on FeO*-MgO, 
translation to another composition involves small changes in 
length (up to 5%) and small rotations in orientation (up to 
10°). The labels correspond to: o - olivine, PAR, Fos3;; c - 
clinopyroxene, PAR, DiscHd,3; p - plagioclase, PAR, ANn¢o. 
All three of these phenocryst vectors correspond to 10% by 
weight subtraction. The vector labelled i corresponds to 
ilmenite, PAR, and is for 3% subtraction. Vectors a and cc 
are alteration vectors corresponding to 10% replacement of 
groundmass glass by a=septechlorite, PAC, and 5% replacement 
by cc=calcite respectively (addition vectors). These 
calculations were made particularly for the oldest lavas to 
distinguish the effects of alteration related to chemical 
variation from primary igneous processes. 

On examining this figure (7-12), it is apparent that 
most of the oxide variation can be accounted for by <10% 
Olivine fractionation, <15% clinopyroxene fractionation and 
<6% ilmenite fractionation. The effect of chlorite 
alteration is in fact away from most of the variation 
trends. One flow of this set, PAW, had amygdules of calcite. 
This point stands apart at 12.46% on the CaO axis and its 
separation from the main lava group can be accounted for by 
5% calcite alteration. No other lavas show this effect. 
Linear trends appear on MgO versus FeO#, FeO* versus Ti0, 
and to a lesser extent on MgO versus Al,O;. It must be 


emphasized that the sequence of points along these trends 
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does not correspond to eruptive sequence. This lack of 
timewise order for trends on the variation diagrams is a 
persistent theme at Level Mountain. It would appear that the 
eruptive sequence samples a related suite of compositions, 
and while they appear to represent a fractionation sequence, 
the points along a given trend are not erupted in order. In 
most flows, appropriate degrees of fractionation of observed 
phenocrysts could account for the compositional variation. 
When sampling was performed the dense phenocryst-poor 
portions of flows were selected. The existence of flowage 
differentiation phenomena has been demonstrated by Komar 
(1972). The source magmas responsible for map units 1 and 2 
could have been relatively uniform alkali basalts, but the 
existence of phenocrysts and the various flow and sampling 
processes, not to mention variable crystallization of the 
sampled flow portions, behave like dispersion vectors to 
string the compositions out. 

Major element variation diagrams for Plateau unit 3 are 
presented in figure 7-13. The composition used in 
fractionation tests was flow PAM. It is a K-poor tholeiite 
in the sense of Irvine and Baragar (1971); however, the flow 
is petrographically an alkali basalt and has an alkalinity 
index of 1.35. Of the twelve compositions plotted for unit 
3, five classify as K-poor tholeiites, one is a hawaiite and 
the other six are alkali basalts. Not only do the lavas 
appear out of stratigraphic order along the variation 


trends, but the alkali basalts and "tholeiites"” are 
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Figure 7-13. 
Oxide variation plots (anhydrous weight percent) for plateau 
basalts from map unit 3. The lines radiating from a common 
point indicate mineral fractionation directions. The two 
lavas joined by a common point are for (29/1f£) an alkali 
basalt flow and its filter pressed derivative. This 
indicates that most of the compositional variation here is 
probably not related to protracted fractional 
crystallization of several phases. 
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thoroughly interspersed. The fractionation vectors used were 
10% for: c - clinopyroxene, PAO, Di,.,Hd,,; o - olivine, PAR, 
FOs3;; and p - plagioclase, PAR, An,, and 2% for m - 
titanomagnetite, PAC, Mt,,Ulv,.. Most of the variation here 
can be accounted for by 15% olivine fractionation with less 
than 5% titanomagnetite, which is in agreement with 
phenocryst mineralogy. In this figure (7-13) there is an 
additional arrow drawn to connect alkali basalt 29/1f (an 
ophitic textured, ponded, columnar jointed flow) with its 
filter pressed derivative hawaiite (that occurred as narrow 
dikelets along the columnar joints). The arrow points 
towards the hawaiite and its length represents greater than 
60% fractional crystallization of olivine, titaniferous 
augite and plagioclase. Crystallization (fractionation) of a 
Single phase can usually account for most of the trend, but 
combined fractionation of several minerals or less produces 
off-trend variation. This distinction 1s very evident for 
the filter pressing case. In comparing the sequential map 
units another interesting twist is seen. The dispersion 
trends are more or less parallel and oriented along parallel 
mineral vectors but there is a shift in level between 
sequential units, compare for instance the FeO* versus S10, 
plots or the MgO versus Al,0O; plots in figures 7-12 and 
7-13. The shift between successive eruptive pulses is like 
the off-trend variation produced by protracted fractional 
crystallization of several phases. Perhaps sequential units 


are produced from different sources altogether or from the 
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Same source which had undergone something on the order of 30 
to 40% fractional crystallization between the eruption of 
‘SuccesSive map units. If this speculation has any merit then 
there would have to be very long lived large volume magma 
chambers at subcrustal depths which persist for hundreds of 
thousands to millions of years. 

Major element variation diagrams for plateau unit 4 are 
presented in figure 7-14. The composition used in 
fractionation tests was flow PAF, a quartz normative 
hawaiite. Of the ten points on these diagrams, four are 
alkali basalts and the remainder are hawaiites. Of the 
hawalites, two are nepheline normative and one is quartz 
normative. Two linear groups of points can be seen on the 
Na,0-Al,0, plot. All of the ordinary hawaiites and alkali 
basalts comprise the trend at the lower Na,O level. The 
quartz and nepheline normative hawaiites together with the 
remaining ordinary hawaiites comprise the trend at high 
Na,0. The lavas of unit 4 erupted nearly in timewise order 
along these trends. On the MgO-Al.0;, MgO-CaO and MgO-FeO* 
plots, the eruptive sequence more or less proceeds from 
hawaiites to alkali basalts with a steady increase in MgO. 
There are also lineations on S$iO,-FeO* and S$10,-Al.,0; with 
no time For this map unit clinopyroxene, and to a lesser 
extent plagioclase and titanomagnetite, appear to be the 
controlling vectors. The vectors are 10% for: o - olivine, 
PAF, Fo;;; p - plagioclase, PAF, An,o; c ~ clinopyroxene, 


PAF, Di,;Hd;,; and 2% for m - titanomagnetite, PAC, 
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Figure 7-14, 
Oxide variation plots (anhydrous weight basis) 
plateau lavas from map unit 4. 
discussed in the text. 
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mt,,;ulv,. These lavas frequently contain large cumulophyric 
clots of black clinopyroxene plus plagioclase, interpreted 
to be cognate. Unlike the previously discussed map units, 
the variation cannot be explained without in excess of 40% 
fractional crystallization, implying longer time intervals 
between eruptions than for earlier units. An eruptive 
sequence from hawaiite going to alkali basalt could imply 
top to bottom emptying of a magma chamber (which had 
differentiated to hawaiites at the top, with related alkali 
basalts deeper). Alternatively, the onset of a new eruptive 
cycle could create a magma chamber which differentiates to 
hawaliite. The repeated injection of fresh related alkali 
basalts releases the hawaiite as injection of fresh magma 
proceeds, until finally the alkali basalts come out as the 
eruptive pace picks up and insufficient time is available 
for further differentation. Whichever of these possible 
mechanisms is responsible, the lavas of unit 4 represent the 
first major outpouring of hawaiites. 

Figure 7-15 portrays the compositional variation for 
the basic lavas of the stratocone units. Because of the 
large variation in lava compositions represented here, the 
test composition used for the fractionation arguments is not 
a specific lava but the average of the Level Mountain 
hawaiites. Fractionation compositions were selected from 
lavas with similar element ratios to the average. The 
vectors used in fractionation tests are: 


p - plagioclase, PBX, An,,, 10% 
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Figure 7-15, 
Oxide variation plots (anhydrous weight basis) for basic 
lavas of the stratocone Stage. For discussion of chemical 


Variation by map unit and explanation of fractionation 
VECTOrs, refer to text. 


- ~~ 
St-\ aug . 54 
: ‘ ‘ “ve ‘Anal zzol 

, ' » as 4 @ 
u SHg I] ow , Se Ro she 
> (gee dri. : <o¢e 

1Oo Ll RSe vet | 7 “2 ‘in 
4) FH ' 4 Qa 
. { > 


210 


- clinopyroxene, PBX, Di.,Hd.,, 10% 

- olivine, 8/20-2/4307, Fos.rim, 10% 
titanomagnetite, PAC, Mt,;, Ulv., 5% 

- ilmenite, PBX, 5% and 

- quartz, 5% addition to simulate the effect of 
crustal syntexis/contamination. 
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For many of the plots the oxide composition and vectors are 
orthogonal. As in the previous figures, the distinction of 
lava compositions between eruptive units is readily made. 
Again the arrangement of points along trends is not ina 
timewise order either within a map unit or for sequential 
map units. 

The alkali basalts and hawaiites of map unit 5b 
petrographically contain phenocrysts of clinopyroxene (black 
augite) and plagioclase, sometimes in clots as in flows from 
map unit 4. The analyses plotted come only from the sequence 
from Meszah Peak Kakuchuya Valley. The unit commences with 
hawalites, progresses to alkali basalts, then back to 
hawalites. The distribution of the unit 5 points is 
consistent with clinopyroxene plus plagioclase fractionation 
each around 10%. When the map unit 5b lavas are compared to 
those of map unit 4b it can be seen that they are chemically 
distinctive and do not represent a timewise progression of 
the unit 4 trends. 

The basic lavas of unit 6b were more widely sampled 
than for 5b. The sequence at Meszah Peak begins with 
hawaiite and progresses through quartz normative hawaiite to 
transitional hawaiite-alkali basalt. The other lavas sampled 
from Kaha-Lost Creek and Dudidontu Creek were all quartz 


normative hawaiites. These lavas contain olivine or 
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iddingsite plus titanomagnetite phenocrysts with infrequent 
plagioclase. There chemical variation trends are consistent 
with olivine plus oxide dominated fractionation, with total 
range of fractionation being less than 20%. Compared to the 
basic lavas of 5b, the 6b lavas are all on the quartz 
enriched side. A vector for alkali feldspar was not plotted 
but for 10% addition this vector would be roughly parallel 
with the plagioclase fractionation vector and points the 
opposite direction. Contamination of basic magma with a 
minimum melting composition granite would be along a vector 
between quartz and alkali feldspar. This is also the 
direction for incorporation of comendite as magma mixing. 
The 6b basic lavas could differ from their 5b predecessors 
by less than 5% crustal contamination or magma mixing with a 
more Salic melt. Both of these possibilities imply the 
existence of crustal level magma reservoirs as did the 
upward basifying trend roughly shown in the basalts of units 
4, 5b, and 6b. 

No basic lavas were analysed from any of the stage 7 
units a, b or c. During this period basic compositions were 
Subordinate and no well developed basic stratigraphic 
sequence was encountered. However in three localities 
marbled lavas (admixed basic and salic) were encountered. 
The admixed salic lava was trachyte but adjacent 
stratigraphic units were metaluminous phonolite ash flows, 
tuffs or plugs. These occurrences demonstrate the 


contemporaneity of basic and salic melts and lend further 
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credence to mixing hypotheses in the derivation of the 
intermediate lavas. 

Only two basic lavas were analysed for unit 8, a 
transitional hawaiite-alkali basalt flow from Meszah Peak 
and an alkali basalt dyke from Egnell Creek on the 
southwestern plateau margin. The hawaiite contained olivine, 
plagioclase and titanomagnetite phenocrysts in addition to 
about 0.1% by volume of xenoliths. The xenoliths were 
dominated by two varieties: (i) a plutonic textured 
troctolite (olivine and plagioclase with minor 
clinopyroxene) and (ii) granitic gneiss. The latter crustal 
fragments showed disaggregation and partial melting at grain 
boundaries. The hawaiite is more silica-rich and iron-poor 
than the alkali basalt (which could be explained by less 
than 5% crustal contamination). The alkali basalt contained 
megacrysts of complexly zoned plagioclase and three or more 
generations of plagioclase phenocrysts. This lava is richer 
in CaO which is consistent with plagioclase fractionation. 

The flows of map unit 9 begin as nepheline normative 
hawaiites and proceed to slightly more silica rich ordinary 
hawaiites. Their order along the trends is not their 
eruptive sequence. These lavas are almost totally aphyric, 
although one locality on the west shoulder of Meszah Peak 
contained sparse xenocrysts or phenocrysts of quartz, 
partially resorbed. The slight variation in the composition 
of these four flows is in keeping with either less than 5% 


Olivine fractionation or less than 2% variable crustal 
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contamination. The MgO-Al.0;, and MgO-CaO plots argue 
strongly for olivine fractionation while MgO-SiO, and 
FeO*-SiO, favour crustal contamination. Some combination of 
the two is considered most likely. As with the unit 8 lavas, 
these are more primitive than the basic lavas of units 5b 
and 6b, and 4. The unit 9 lavas are enriched in TiO,, K,O, 
and P,0;,. The hawalites from map unit 9 could be 
differentiates from a larger volume of more primitive alkali 
basalts having no surface manifestation. 

Selected analyses were tested against a suite of 
hypothetical source compositions chosen from Level Mountain 
and other B.C. inclusions: lherzolite, troctolite, eucrite 
and gabbro, table 7-13. Other peridotites were also tested. 
Pyrolite was used as the hypothetical composition for the 
derivation test because the actual nodule compositions for 
B.C. were too depleted to give rise to most of the basalt 
compositions. The B.C. lherzolites, typically lacking any 
P.O, and low in TiO,, gave possible derivative melts less 
than one percent by volume. Derivation calculations were 
also made using LSPX, a Fortran least squares matrix 
calculation routine written by F. Chayes. For any of these 
calculations, the lava composition is related to a 
hypothetical source such as a peridotite, through 
arbitrarily selected mineral compositional vectors. 
Obviously the assumptions of source and phase compositions 


severely restrict the significance of these calculations. 
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Table 7-13. Model Source Compositions 


SiO, 
TiO; 
AisQ; 
Fe.,0,; 
FeO 
MnO 
MgO 
CaO 
Na.,0 
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H;O+ 


CM 


Class- 
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9568 
vr 
$37.05 
0.16 
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The hawaiites, and lava types transitional between 
hawaiite and alkali basalt, could always be distinguished 
from more primitive basalts on the basis of the value for 
Maximum percent partial melting from pyrolite. The cutoff 
was 10.5%, with hawaiites below and alkali basalts, 
ankaramites and hypersthene normative basalts above. The 
average percent partial melting for all hawaiites was 
7.7£1.9% with a total range of 4.4% to 10.3%. The average 
percent partial melting for all primitive basalts was 
14.042.7% with a total range of 10.6% to 18.4%. The ranges 
Calculated for alkali basalts are similar to the 15 to 20% 
range cited in Fiesinger and Nicholls (1977). Among the 
other source compositions tested, the Telegraph Creek 
lherzolite (Littlejohn and Greenwood, 1974) was so depleted 
that the amount of partial melting was never greater than 4% 
for any derivative composition, and usually the solution was 
zero. From this is concluded that lherzolite inclusions of 
this type may represent depleted or residual mantle but are 
not very likely basalt sources. The troctolite inclusion 
from Meszah Peak gave answers ranging from 18% to 32% for 
the Level Mountain basalts which is considered to be too 
high. 

Fractional crystallization calculations were also 
performed to determine the minimum amount and proportion of 
Subtraction of phenocryst phases which could relate two lava 


compositions. An example calculation was made on the Level 
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Mountain lavas using the average alkali basalt to yield the 
average hawaiite, through fractional crystallization of the 
Same mineral vectors used in figure 7-12. The solution, 
Ppresenced in table 7-14, gdve 7333%9F0,"4)° 18.75% Di, eHd;,., 
feweneans,,, and 3.5% MtE,.,Ulv6 for’ a total fractional 
crystallization of 43.99, all expressed in weight fractions. 
While the values for some oxides of the solved hawaiite may 
be outside the range of one standard deviation, they are 
within range for the hawaiite data set. The misfit is 
attributed to inappropriate choice of mineral compositions. 
Recalling the discussion of the preceeding section, the 
inferred degree of fractionation within map units was 
generally lower than this but in fact the difference between 
SuccesSive eruptive pulses was of this magnitude. If 
fractionation of this magnitude is indeed required for 
hawalite derivation, then the amount of related plutonics 
beneath Level Mountain and the other centers with voluminous 
hawaiite may approach the surficial lava volumes in 
magnitude. Another way of viewing this is with the partial 
melting calculations which demonstrated that roughly half 
(45%) of the source amount was required to yield the 
hawaiites compared to the basalts. 
Intermediate Lavas 

It has been repeatedly stated that Level Mountain has a 
bimodal suite of basalts and salic lavas. There are however 
lavas of intermediate composition including mugearites, 


phonolites, benmoreites and tristanites. While never 
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abundant, they are best represented in the stratocone, units 
Sa,6b,7c,8. Not only are these lavas chemically diverse but, 
with one exception, they are all metaluminous or 
peraluminous and not obviously related (via an ordinary 
differentiation sequence) to the sodic alkali basalt series 
on one hand, or to the peralkaline trachytes and comendites 
on the other. The major element chemistry of these lavas is 
depicted by the histogram in figure 7-1 and tabulated in 
table 7-15. A cursory look at the range of compositions for 
each oxide in this group will show them to be intermediate 
between basalts and salics; thus the name. 

The compositional variation for these lavas was 
presented on a set of eight variation diagrams, figure 7-16. 
Different symbols are used to keep track of eruptive 
sequence. On the variation diagrams, the effects of 
fractionation are shown for a transitional benmoreite, flow 
PBO from unit 8. The fractionation vectors are: 

-TO.1vine, PBO,, Fo,;, 4% 
- clinopyroxene, PBO, Di,,.,Hd,.Acs, 5% 
plagioclase, PBP, An,,Ab,,0re, 10% 


- alkali feldspar, PBO, An,;Ab,,0r,;, 10% 
- ulvospinel, PBO, Mt..,Ulv.., 5% 


cooano 
| 


Linear trends appear on most plots but again there is no 
timewise order along the trends. Unlike in the previous 
diagrams, these vectors do not simply parallel the variation 
trends. In fact the fractionation vectors don't even satisfy 
the variation of unit 8 benmoreites very well. A combination 
of all five phases would be required to describe the 


existing variation for all intermediate lavas for a total of 
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Table 7-15. Miscellaneous Intermediate Lavas. 
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Figure 7-16. 
Oxide variation plots (anhydrous weight basis) for 
intermediate lavas of the stratocone stage. Most chemical 
variation for these disequibrium textured lavas fits a 
mixing origin from basalt plus trachyte. 
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15 to 20% fractionation, with ulvospinel apparently 
dominant. There is an additional problem in attempting to 
translate these vectors to other points on the diagram. The 
phenocryst compositions and populations of these lavas are 
all different. 

The mugearites 26/3c and LMI20d had two phenocryst 
(xenocryst) clinopyroxenes; large partially resorbed 
titaniferous augites and small enhedral sodic 
ferrohedenbergites. The unit 7a series phonolites 
(stratigraphic section LMI20 samples) had cumulophyric clots 
of anorthoclase, with or without green clinopyroxene, but 
always with riebeckites which sometimes occurred as mantles 
on fayalite cores. Most of the intermediate lavas also 
contained resorbed plagioclase phenocrysts of andesine 
composition, a feature also noted for some hawaiites. 
Mantled feldspars and disequilibrium multiply zoned 
feldspars have also been used as evidence for magma mixing, 
(Hibbard, 1981). Most of the phenocryst assemblages do not 
appear to be in equilibrium with their groundmass and 
regardless, they still fail to describe the variation trends 
on the plots. 

The time-map units which contain the intermediate lavas 
always contain some other dominant lava type, usually basalt 
or trachyte. In unit 7c near the headwaters of the south 
fork of the Beatty Creek, basalt and trachyte lavas occur 
together as mixed flows. The other flows of this unit are 


trachyte or the intermediate types. Crustal gneiss xenoliths 
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have been found in trachyte of unit 7c and benmoreites of 8. 
Perhaps, as field evidence suggests, the intermediate lavas 
Originate by magma mixing with variable crustal 
contamination. 

The phonolites of unit 5a were sampled from the north 
wall of the Kakuchuya Valley. All were light coloured rocks 
with trachytic texture. 9/1-86/5421 was the thickest of such 
flows sampled, about 24m, with columns 1.5m across at the 
base that flaired toward the top of the flow. Phenocrysts 
included anorthoclase to 1.5cm, abundant biotite 0.9cm x 1mm 
thick subhedral platelets and 0.5cm prisms of black 
amphibole (barkevikite or basaltic hornblende). Flows PBV 
and PBU also contained microphenocrysts of amphibole and 
biotite. The only other lavas which contained phenocryst 
amphibole and biotite were the tristanites of unit 7. PBU 
and PBV are chemically similar to phonolites of the 
trachyte-phonolite association in general (see table 2-3 of 
Carmichael, Turner, and Verhoogen, 1974). Compositionally 
9/1-86/5421 is distinctive in being the only peralkaline 
phonolite sampled, having A.I. (Na,0+K,0)/A1,0; mol = 1.01 
and ns calculated in the norm. Chemically this sample is 
very similar to phonolites from oceanic islands, see for 
instance the peralkaline phonolites from Tenerife, Canary 
Islands (eg. TC167 Edgar and Parker, 1974) and to 
"undersaturated trachytes" from Cinquo Picos Terceira, 
Azores (Self and Gunn, 1976). Two hypothesis are given for 


the genesis of the Tenerife Suite: fractional 
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crystallization or fractional melting to give two primary 
melts such as discussed by Yoder (1973). If bimodal 
fractional melting describes the genesis of alkaline 
volcanics like Terceira and Level Mountain, then phonolites 
could be the result of magma mixing or they could even be a 
primary melt, (Takahashi and Kushiro, 1981). 

"Petrogeny's residua system" for undersaturated alkalic 
melts is the compositional plane Q-Ne-Ks. The Level Mountain 
phonolites plot along the thermal minimum through between 
the alkali feldspar join and the one atmosphere reaction 
point, see figure 7-17. The phonolites and feldspathodial 
trachytes of Dunedin Volcano, Price and Chappell (1975) show 
a Similar scatter about the thermal valley but a better 
defined trend toward the reaction point. The one atmosphere 
reaction point and thermal minimum are according to Bowen 
and Schairer (1938) and Schairer and Bowen (1956). The one 
kilobar data are from Hamilton and MacKenzie (1965). The 
four metaluminous phonolitesof Level Mountain plot nearer to 
the alkali feldspar join, while the peralklaine phonolite 
approaches the reaction point. Once a melt reaches this 
system, further alkali feldspar fractionation will drive 
residual liquids to the reaction point. These lavas are just 
across the silica saturation barrier and "on trend" for the 
more voluminous trachytes. 

The benmoreités and transitional benmoreite-trachytes 
of unit 8 outcrop over approximately a 100m interval on the 


south face of Meszah Peak. Here alkali feldspar phenocrysts 
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igure 7-17. 
Distribution of Level Mountain phonolites in the 
underSaturated portion of petrogeny's residua system (weight 
basis). Ris the ikbar reaction point after Hamilton and 
MacKenzie (1965). R' is the latm. reaction point after 
Schairer and Bowen (1935). Note that the Level Mountian 
phonolites plot between the one atmosphere reaction point 
and the thermal minimum on the alkali feldspar join. The 
peralkaline example is closest to the reaction point. 
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were present, aS were gneiss inclusions. While the field 
appearance of this scoria was rather different than the 
aforementioned flows, the chemistry of the benmorite was 
virtually indentical, compare analyses PBO, PBP and 
8/25-55/6722. The alternating sequence of benmoreite with 
trachyte is good field evidence of a mixing origin for the 
benmoreites. Some crustal contamination may also be involved 
as evidenced by the gneiss inclusion. Representative 
analyses of intermediate lavas from elsewhere in the IMB are 
presented for comparison to the Level Mountain intermediate 
lavas in table 7-16. The striking feature of all these 
analyses compared to associated basalts and trachytes, are 
the high levels of Al.0;, typically in excess of 17.0%. The 
mugearite from Edziza was originally described as 
trachybasalt. Similar mugearites from the Rainbow Range are 
described by Bevier (1978), as phenocryst rich (10 to 25%) 
thin red, grey and brown flows in sections with abundant 
breccias in excess of 70%. These petrographic and field 
descriptions are virtually identical to the mugearites and 
benmoreites at Level Mountain. These were also originally 
described as trachybasalts by Hamilton and Scarfe (1977). 
The mugearites at all three centres (Level Mountain , 
Edziza and Rainbow Range) occur at high enough statigraphic 
levels to be interspersed with comendite or trachyte acid 
lavas. Nonetheless, a mixing hypothesis for the origin of 
these intermediate lavas has not been commonly invoked. At 


Edziza, Souther and Symons (1974) explained this chemistry 
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to be representative of the "primary magma" for that 
eruptive unit. While for the Rainbows, Bevier (1978) 
explained their derivation by 70.71% fractional 
crystallization of a hawaiite parent. It is doubtful that 
mugearite could be a primary (unfractionated) mantle derived 
magma particularily for the Level Mountain and Edziza 
examples which are devoid of P.O, and very enriched in Si, 
Al, Na and K relative to undeniably primitive basalts. 

The three phonolite analyses presented here, all from 
flows rather than from tuffs, are chemically quite a 
coherent group. These examples from the Itcha Mountains and 
Mount Edziza are most similar to the flows PBU and PBV at 
Level Mountain. They appear to be chemically similar or 
transitional to tristanites. Note the small amount of 
difference between SS5, a corundum normative phonolite and 
SS3, a tristanite from Edziza, (table 7-16). At Level 
Mountain the tristanites have even higher levels of Al.,0,. 
While there was no petrographic evidence, such as 
xenocrysts, to support a magma mixing hypothesis for the 
Origin of the Itchas phonolites, the potential for crustal 
contamination is still there. 

A single stage magma mixing model is presented in table 
7-17 for each of the three intermediate lava types: 
mugearite, benmoreite and phonolite. These mixing models are 
not very sensitive to either the basaltic composition used 
or the trachyte composition, as the proportions work out to 


be almost identical. While the points for average hawaiite 
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and average trachyte have not been plotted on figure 7-16, 
they span the dominant variation trend, with a mixing line 
that 1s nearly to it. Similar mixing relations have been 
discussed for alkaline volcanics by Brooks and Printzlau 
(1978) and by Anderson (1976). 

In the case of the mugearite model, a basaltic magma 
and a trachyte magma must be mixed in subsequal proportions. 
Starting compositions of both alkali basalts and hawaiites 
have been used for a total range of basalt weight fraction 
from 0.4944 to 0.5234. The biggest discrepancies for the 
case shown are for T1i0O.,, P.O, and CaO. The calcium problem 
can be solved by varying the composition of the basaltic 
parent within one standard deviation. The TiO, and P,0O, 
excesses occur pervasively in all mugearite calculations 
attempted. However, they can be explained if in addition to 
mixing of basalt and trachyte in model proportions there is 
somewhere between 0.25 and 1% apatite fractionation and less 
than 3% fractionation of an iron titanium oxide. 

Benmoreites are restricted to map units 7 and 8, soa 
hawaiite from that period is presented to model their 
derivation. Again the choice of hawaiite composition is not 
really that critical to the mixing proportions. The fit for 
the model presented here is quite good. The only 
discrepancies are the Na,O and Al,O; which are small enough 
to be accounted for by alternate choice of parent trachyte 
or by less than 2% of assimilation of granitic material. 


Because these flows contain granitic gneiss xenoliths that 
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are partially melted either solution is considered 
acceptable. 

The phonolites are compositionally quite variable but 
the mixing proportions for basalt and trachyte vary less 
than 2% for any choice of parent lavas. When these 
calculations are made for phonolites PBU and PBV the 
addition of crustal material is required in the 5% range to 
account for their higher levels of SiO,, Al.,0O;, Na,O and 
K,0. These particular flows are nearly as enriched in Al,0O; 
as the tristanites. There are no other magmas at Level 
Mountain with sufficiently high levels of Al,0; to explain 
such products by magma mixing alone. The lack of a suitable 
parent composition exists for Na,0 as well, in the case of 
9/1-86/5421, which is a peralkaline phonolite. Presumably 
this points to open system behavior for alkalis. The mixing 
test for phonolite genesis only partially explains their 
major and trace element chemistry. In addition to magma 
mixing there would have to be a component of crustal 
contamination and open system behavior for alkali and 
incompatible elements. 

The trachyte compositions which satisfy these mixing 
models are all metaluminous. The Al,O; and Na,O levels of 
the peralkaline trachytes are too low to meet the daughter 
lava compositions. This presents a bit of a problem as the 
primary trachyte magmas at Level Mountain are thought to be 
peralkaline. The admixed xenocrysts in most of the phonolite 


lavas such as anorthoclase, sodic ferrohedenbergite and 
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riebeckite are all derived from peralkaline trachyte 
parents. There is another way to visualize the metaluminous 
trachyte average used in these calculations; that is a 
peralkaline trachyte which has experienced volatile loss and 
crustal Saedanidat don. Wuhese complications could accompany 
the magma mixing and eruptive process. 

The best explanation for the low volume compositionally 
variable intermediate lavas is that they lie on a mixing 
line between the more common hawaiites and trachytes, with 
variable amounts of crustal contamination. When fractional 
crystallization models were attempted for the intermediate 
lavas from Level Mountain using a basaltic parent (basalt or 
hawaiite) and analysed phenocryst compositions there were no 
acceptable solutions. Additionally such fractionation models 
fail to explain the disequilibrium phenocryst assemblages, 
and they do not make allowance for crustal contamination 
such as the high Al.,0O;, hydrous minerals and crustal 
xenoliths would indicate. Mixing models presented here can 
account for most of the compositional variation without 
recourse to extensive additional fractionation or 
contamination. 
Major Element Petrochemistry of Salic Lavas From Level 
Mountain and a Comparison of Similar Lavas From Other 
Localities in the Intermontane Belt. 
Introduction 

The volumetrically dominant salic lavas are 


peralkaline. Chemically they span a range of more than 10 
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weight percent in SiO,, including types which have variously 
been described as comendites, comendite trachytes, 
pantelleritic trachytes and pantellerites. Physically they 
range from aphyric pitchstones, to tuffs, ignimbrites, 
stocks, laccolits and holocrystalline eutaxitic trachytes. 
Chemically their hallmark is having molecular 
(Na,0+K,0)/Al1,0,>1. This agpaitic index greater than unity 
implies excess alkalis over the amount required to form 
feldspars and manifests itself as the occurrence of alkali 
ferromagnesian silicates. 

For the purposes of minimizing the chemical variation 
within salic rock classifications, only four groups are 
recognized. Peralkaline trachyte includes the types 
comendite trachyte, pantelleritic trachyte and pantellerite 
of MacDonald (1976). Most of the Level Mountain varieties 
are of his pantelleritic trachyte type. These lavas lie more 
than 98% in the system (SiO,, Al.0;, FeO*, Na,0, K.,0O) with 
Silica content in the mid to upper 60's and more than 12.5% 
for normative colour index. Comendite corresponds to 
MacDonald's (1974) definition and is a weakly peralkaline 
rhyolite with silica content in excess of 70%, and total 
iron less than 5% which gives rise to a normative colour 
index less than 12.5%. In these lavas MgO and CaO are 
reduced to the level of minor or even trace constituents. 
Normative recognition of these peralkaline lavas is based on 
the occurrence of Ac and/or Ns, and depending on the CaO 


content, Wo is generally present as well. 
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The trachyte group includes both alkaline 
(CaOt+Na,0+K,0)/Al1,0;>1, and metaluminous types. They range 
from wollastonite to hypersthene normative and ina few 
cases include normative corundum. The alkaline types 
dominate. Petrographically and mineralogically they include 
types with phenocryst fayalite and hedenbergite and no 
unusual groundmass mineralogy in addition to lavas whose 
groundmass phases include the same alkali ferromagnesian 
minerals characteristic of the peralkaline trachyte group, 
aenigmatite and alkali amphiboles. The chemistry of this 
trachyte group contrasts slightly with the peralkaline 
trachytes in having lower Si0, and FeO*, higher Al.,0O; and 
higher levels of all the minor components CaO, TiO,, P20s, 
MgO and MnO. These chemical characteristics for the alkaline 
trachytes lead to a more diverse and apparently more 
strongly peralkaline groundmass than in the case of some of 
the peralkaline trachytes. 

The rhyolites are chemically quite variable, ranging 
between comendites and trachytes for most oxide 
concentrations but they are distinguished by low Na,O. In 
fact the weight ratio of K,0/Na,0 for this group is 
generally greater than one. This is highly peculiar compared 
to any other lava type at Level Mountain. These rhyolites 
range from metaluminous to peraluminous and have normative 
hypersthene with or without corundum. Rhyolites are 
Spatially and temporally more restricted than comendites. 


The most comon occurrence of rhyolites is as aplitic tuffs 
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and partially devitrified pitchstones. Glassy commendite 
dykes commonly have a crystalline white rhyolite selvedge, 
which would seem to indicate that rhyolites are altered 
forms of more common comendite melts (Noble, 1967). 
Peralkal ine trachytes 

Chemical analyses and norms of the peralkaline 
trachytes are presented in table 7-18. Of the thirteen 
analyses in this group, MacDonald's classification (1974) 
would call six of the rocks peralkaline trachytes of 
pantellerite affinity, three peralkaline trachytes of 
comendite affinity, and four pantellerites. The thirteen 
analyses represent all three salic eruptive cycles (map 
units 5a, 6a, and 7) of the stratocone. Most of the 
compositional variation can be seen in the samples from unit 
6a, SO no Stratigraphic association or genetic sequence is 
inferred. Comparing the averages for the three subgroups 
Shows them to overlap within one standard deviation for all 
ten oxides, table 7-19, so the three groups were 
Subsequently combined. For the entire population, Si0, 
varied inversely as Al,0,, Na,O and P,0O;,, while Al,0; varied 
normally with Na,0. Such variation could be explained by 
alkali feldspar plus apatite fractionation. A system open to 
SiO, could also partially account for these covariances. The 
minor components MgO and MnO vary normally with each other 
and with FeO*. Random variations in content of phenocryst 
Olivine (Fa rich) and possibly clinopyroxene could account 


for the low and highly variable levels of MnO and MgO in 
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Table 7-18. Chemical 
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addition to their covariance with FeO*. their covariance 
with FeO*. For confirmation of this, recall the enrichment 
of MnO in K,O varies inversely with MgO, MnO and CaO, 
reflecting the variation between phenocryst populations 
dominated by anorthoclase, versus those dominated by 
phenocryst fayalite and sodic ferrohedenbergite. TiO, varies 
normally with FeO* and inversely with Al,0O;. Sometimes there 
is phenocryst ulvospinel or ilmenite associated with olivine 
and clinopyroxene, although generally the peralkaline 
magmas, particularily at Level Mountain, seem to have 
crystallized in an oxide absent field as discussed by Ernst 
(1962) and Marsh (1975). The FeO#-TiO, covariance could 
indicate that some iron titanium oxide fractionation may 
have played a role in the genesis of the Level Mountain 
peralkaline trachyte magmas, while the combined effects of 
temperature, oxygen fugacity and peralkalinity preclude 
oxide phases during groundmass crystallization. 
Comend ites 

Chemical analyses for representative comendites are 
presented in table 7-20. In reviewing the figure 7-1 
histograms for comendites, the distributions appear to be 
normal for all oxides except $i0,, K,0 and MnO. The 
population was arbitrarily divided at the means for these 
oxides. Means and standard deviations were calculated for 
high and low splits to see if any statistical differences 
were indicated. The results are presented in table 7-21. The 


high silica set was lower in Al,0, and marginally lower in 
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Table 7-20. Chemical 
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Ti0,, while being more strongly peralkaline and having a 
higher atomic ratio of Na/K. While the high silica, high 
agpaitic index and low Al.,O; could be explained by alkali 
feldspar fractionation of appropriate composition, the 
higher Na/K ratio negates this possibility. The ratio of 
Na/(Nat+K) for feldspar to whole rock is typically about 
1.12, so that any alkali fedlspar fractionation would lower 
Na/K, which is in the opposite sense to the variation 
described. The chemical variation between the high and low 
Silica groups is better explained by open system behavior, 
particularly for Si and Na, coupled to volatile loss during 
eruption and crystallization, or to subsolidus 
devitrification and alteration of groundmass glass (Noble, 
1968). The clusters for low K,0 and high MnO were coincident 
and showed higher agpaitic index and Na/K ratios. These 
differences could equally well be explained by 
crystallization and also by glass alteration processes. 
Noble (1965, 1967, 1968) and Noble and Haffty (1969) have 
extensively described crystallization, devitrification and 
alteration processes in peralkaline rocks from the western 
U.S.A., Canary Islands and Australia. Lowering of S$i0, and 
peralkalinity index (AI) are universally seen. The only 
elements whose levels may be unaffected by these processes 
are Ba, La, Nb, Rb and Zr. 

Low and variable contents of phenocryst acmitic 
pyroxenes and iron titanium oxides have been 


petrographically noted for the comendite pitchstones. These 
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petrographic characteristics can probably account for most 
of the covariances involving FeO#, Cao, MgO, TiO, and MnO as 
discussed above. The direct variation of K,O with Al,O, 
could be explained either by feldspar fractionation, 
alteration or open system behavior. The inverse correlations 
with SiO, could be explained by alteration or open system 
behavior. While quartz fractionation is possible, the 
occurrence of quartz phenocrysts in comendites was rare. 
Trachytes 

Chemical analyses for fourteen trachytes representing 
map units 5a, 6a and 7, are presented in table 7-22. These 
lavas differ chemically from the peralkaline trachytes in 
having lower SiO, and higher Al,0;. They also differ from 
peralkaline trachytes on average in having slightly higher 
levels of minor elements, T1iO,, MgO and CaO, but slightly 
lower FeO*x. They are metaluminous with an A.I. less than 
unity, but on the basis of alkali ratio they are virtually 
indistinguishable from the peralkaline trachytes. These 
lavas are highly variable and were not very amenable to 
analysis of covariance as a Single set. They range from Di 
plus Wo normative (alkaline), to Di plus Hy normative and in 
the extreme Hy plus C normative. On the basis of agpaitic 
index, the corundum and hypersthene normative types cannot 
be distinguished from their diopside plus Wollastonite 
normative types. The more encompassing alkalinity index 
(CaO+Na,0+K,0)/Al,0, must be used. For the trachytes, SiO, 


varied inversely with TiO,, MgO and CaO. Unlike all of the 


n 70] JOUG@9R Bidetew oss anise recom tations eee 
‘<) i i 
Or? ,OPH ,090 .400% onivlown) eeorehearan ener 
42iw GoH 96 Molisliay 3279116 ett levee 5s au 
ISHOLIGRTL teqehie! vd .sedFie cenisigns 9 


) sa3evr! off. 19. vered melzaye neues oat ae 


io 


Go %O Reisseigetia Yao ber nisigqus ed Bluos as 
eis ,sidiesoqg 421 ooisenoriosi2 <3 28Up si tdv sa0m% 
a 


\ 
ne 
igi aw tetfo1smes fi alayroonedq S27aUp id 22097 


—— 
ee ere yt a ay 
131 astyionts asetivol to? aseylens fe2tmed3) 


12 art ; 9icse Si Ssineestg exs Tone Be 6a o22004 


a 


poy 
£ * 
ing73 seq sd? mma? vi feoimero zediib. i 


i 


BOTs 152246 os{G godt ..d, 14 sadodd tks ,Ot2 ier a 


ai 


vedo id tripkle z ad ni sos19ve no eetvdsesd oni bas rr 
werd 
_ aw 


me shies ahi 


iipile sud Os) bis O9M 4 Oi? \, stnemble gone 10.3 9¥e 
ae 
16 z 4 oa agiw euvonin legen eye Yonik =#003 3S8¥¢ 
a 

Yilavduiv sxe ys at iledle 20 ebesd ada ad sas she! evi 


oak 


ont .2e3ydoe7s sailsdAiszeq sad mez? a eit 
O2 sidsueme ~[ysv jon stew bane sldeiaayv vidpid: 

iC @eor! sone acy jee elipnte a es sonsivavor 3 ae 
Ns Soe ovitagion ys eule@ id o3 ont letibeds avi 
2iJisaps Jo alesd sz ao oi Jaqs0N).o we 
Tones) €ecvi svidsmien aha lauagyed een 
: tinovzel low aulq sbhinqoife bess 
esq] etial etic enizesgoosns » 


«diz eaydonrs ade a Scaggs 
ais i9 tie eAitay se 22 ben 


a 


a ae 


— 


243 


OL 


60 


ae 


80° 


9I 
6P 
1S 
66 


6s" 


(<4 
60 


6z" 


ts 
th 
te 
6 


ow 


$3" 


so 
ee 
oo 
th 
ez° 
eh 
SS 
ev 
9% 
be 
8 
ve 
is" 
ce 


BCEellinwi 


oo 


eoo00on0me- 


ev 
Oe 


mo oan ocoo-mmoooo ww 


© 


zo" 
ve 


ve 
ee 
iS 
8p 
60 
ob 
sé" 
30° 
06" 


--H9O OO 


09 8 
9s} 

vO Lv 
es 8c 


vO eh 


60 


rt 


40 
se 
ve” 
90° 
ia° 
es" 
te 
90 
94 
61 
99 
98° 
ve 
te 


mocmmooomwoodd 


VOELIIWI 


vi 


eet 
seo 
e6°0 
ees 
be 1s 
86°0 
92 °OS 
90° O€ 


96'S 


vi 9b 


e1'o 
670 
s0 0 
ozo 
60'S 
ve Ss 
wes 
ve oO 
twzo 
GE'h 
es ¢ 
vo Sh 
860 
€v tg 


Ores 
Nid 


vi 


eco'4 
Sup 


SBE 
ves 
w3°0 
o6't 
4th 
tL °O 
146°0 
Si 6b 
96 O€ 


ees 


© 
Cs 
© 
- 


@ 
© 
Haw~oo-mMmoood ° 


vs'O 
ere 


0919 
wad 


vi 


ec 0 


vo € 


66°C 
cok 
zwe'o 
Svs 
6e'4 
vive 
60°0 
$9 €S 
og TE 


wee 


- 


z 
r) 
TOwVNNOO- OHOOO osx 


“ 


ad 


ve 


9t°O 


tees 
ces 


Le 4 
cso 
se'oO 
pees 
19°4 
e9o°S 
ote 
8p eS 
96 62 


66°15 


- 


Nn 
o 
-~OoVMNOO-OHOOO ow 


OB8rs 
Aad 


vi 


ve 
9t 
St 
is’ 


s-00 


.9° 


4L°0 
68 0 
bhot 
Ov S 
wee 
ti ts 
49°0€ 


coe 


9L° 


» 


18 
ce 
io] 
61 
91 
94 
te" 
sO 
96° 
'S 
S6" 
93° 
ez 


NOwUM-CO-HNMOooSD 
= 


Spi 
ad 


vi 


° 


tv 


tz 
96° 
86° 
ge” 
03° 
ec’ 
16” 
Ov 
to 


monwowwr= 


biG 
98°0 
se°Os 
tO OE 


zee 


te 


© 
- 


° 
v 
mowr-moo-wmoocnd 


9ee9 


96-9%/8 €E-71/8 7S-G7/8 


vi 


2 
N 
oo 


a 
o 
-OorTOMNO 


EL” 


rt 
- 


4° 
os” 
ve 

OV 
€3° 
os” 
tL 
mw 
te° 
80 

vee 
6c 
6S 

ar 


moMmnMoormnwzwOoood 


Orig 


wi 


tO" 
6h 


BL 
Bh 
Ls 
80° 
6y- 
es 
St 
co" 
Gt 


8 
vL” 
Sv 
03° 


ie coocort*mo oo 


or 
aw 


16° 


| 


rt 
- 


O° 
oO: 
vO" 
80 
es 
40° 
86° 
aw 
1 
es" 
ai 
03° 
i 
) 


wowNnmMooowwoodd 


3 


nO01LS 
Oses 


vo 


ec 
Ls 
gs" 
t6 


o OFrmN 


te 


se Os 


se os 
6S GZ 
esc 

Gl tb 


© 
A] 
a 


°o 
Orr mMNCOOM <000 


v/9% 


vs vs vs LINN 
avn 
NOIL 
1 1 4 -vwoldl 
-ssvi19 
vO 0 coo Ad 
vio who 8 806Pt 0 av 
vert 
68 0 vi'O 94°0 VW 
90°¥ ves iw 
Sve (40: a <1 AH 
6se ees 614 (s4) 
9t 0 6h os 0 (Na) 
twO'p an 4 1a 
sé6'4 ons (S4) 
91°0 yO" (Na) 
164 4O'4 (0m) 
OA 
vi OS €6°4S 76 SS Vd 
66 4 e945 ier (Nv) 
Si 8p GtOS ~t9 4G (av) 
O€ Lz €6°6t ZEe'te 4o 
s1 oO ) 
$6°6 7 a e9°t oO 
s6 as 4e°9) SS°Sh 113% 
ovg 
zoo 10°O s 
ty 'O 90°O0 81°O - OCH 
Le-o Gro 93900 +O7H 
900 so;O oOF'0O SOt%d 
es yp vO S te's Ot» 
sos 16°S OD OZ7&N 
ges 68°O OO'8 oe9 
446-0 60°O0 02 Oo o5w 
€2°0 st°O OF'O ouW 
voor Be Pv r6°0 oe4 
6L°¢ w6°O +=96'E eote4 
16 vb 6L St 4G Lt eotiv 
or 0 Bse°O Ovo work 
19°¥9 68°69 vi vo tors 
oees 


96-2/6 IZ/Pe IH /Pe 


"B@PAYOR4, 4OY BOBAL BUY ([BO;WEYD “7ZZ-L @1QRUy 


_ 
_ 7 
- 
- 
a! 
5 
m 
— 
\ 7 
‘ P * 
= eeeseenes é 1 22520¢ 
" ee *=ee@ee a ° eaecbeas 
‘ ¢ 
[s2e4 Gtse s 1 55585 
: a + oss 


ts: ssacedss & : capegrseaas 
° ’ = T| 


9 v~s « 9 een 
ctactsee’ ¢ oa darsestecwaey: 
, & -o G8 0eh~Oownee s: 
122 % 33 %e* re2tocaee 

: 22g re eo949~ ee | 


schsceeee a4 | Oe Seseeeeees Pris af) 
>- #0 PEE Hobs. 
: ~ 
=s g eects: ? 23) 
: eared nbd © , ae 


rr ey ; 


ao 2835? 1243 


244 


other salic lavas, there was no negative SiO, - Al,O, 
correlation. TiO, varied normally with CaO and P,O,. Al,O;, 
varied inversely with MnO. FeO* varied normally with MnO and 
both varied inversely with Na,O and K,0. The absolute levels 
of Na,O and K,0 are actually higher on average for the 
metaluminous trachytes than for peralkaline ones. It is 
difficult to predict a simple fractionation, assimilation or 
alteration model to account for all of these variations. 
Rhyol ites 

The chemical variation in the rhyolites is the most 
extensive of the salic rock types. Chemical analyses are 
presented in table 7-23. All examples are limited in time to 
Statocone unit 7. No rhyolites were noted for other map 
units. In a normative sense, these include Hy plus Di, Hy 
and Hy plus C types. According to Irvine and Baragar's 
(1971) classification these are tholeiitic rhyolites of the 
potassium poor series and calc alkaline rhyolites. For the 
rhyolites, SiO, has a particularily wide range. FeO*, Na,O 
and K,0 may show bimodal distributions. High and low group 
averages were calculated for these four oxides. The high 
Silica group and the low FeO* group were identical. Averages 
are presented in table 7-24. Averaging on the basis of SiO, 
and FeO*x levels additionally revealed differences greater 
than one sigma for TiO,, Al.O; and A.I. such that the high 
SiO, cluster was of higher alkalinity and lower TiO,, Al,0O; 
and FeO*x. The minor elements MnO, MgO and CaO were also 


marginally lower (by their standard error ranges) in the 
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Table 7-23. Chemical Analyses of Rhyolites. 


LM13 LMI7 8/7-19 68/16-43 8/16-44 8/27-66 8/27-63 8/27-62 


5830 5990 6350 5862 6788 68390 
$io2 72.80 66.50 67.57 66.58 68.19 74.13 69.62 71.88 
Ti02 0.15 0.43 0.34 0.43 0.56 Rah, 0.0 0.15 
A1203 13.68 14.79 14.96 15.07 14.36 12.93 14.17 13.86 
Fe203 2.66 6.33 2.35 4.06 6.36 0.96 0.55 0.63 
FeO 0.28 0.20 3.78 2.57 0.74 1.25 2.78 2.413 
MnO 0.08 0.08 0.16 0.18 22 0.05 O.11 0.10 
MgO 0.04 0.07 0.0 0.18 0.15 0.09 0.0 0.0 
Cao 0.0 0.34 0.98 0.93 0.74 O.11 0.53 0.55 
Na20 4.78 5.49 3.90 3.57 2.90 3.91 4.33 4.71 
K20 4.76 4.92 5.81 5.49 4.49 5.14 6.28 5.39 
P205 0.02 0.04 0.0 0.03 0.03 0.03 0.0 0.0 
H20* 0.62 0.40 O.11 0.87 i Pea hed 1.28 1.70 0.52 
H20- OH? 0.33 0.03 0.0 0.0 0.0 0.0 0.08 
s 0.01 0.07 0.05 0.09 0.02 
Y%MELT 17.26 16.63 14.23 14.91 18.11 15.94 12.99 15.30 
° 26872 “48707 18.39 22.44 3207 S157 16.94 21.59 
c 0.66 0.47 1.64 3.50 0.80 
OR 28.13 29.07 34.33 32.74 26.89 30.73 35.87 31.85 
(AB) 40.45 46.45 33.00 30.46 24.88 33.59 35.45 39.85 
(AN) 1.18 4.86 4.47 3.48 0.35 0.65 0.76 
PL 40.45 47.64 37.86 34.93 28.35 33.94 36.10 40.61 
(wo) 0.79 0.82 
(EN) 
(FS) 0.89 0.33 
Ol 1.68 1.76 
(EN) O210° 0.17 0.45 0.37 0.22 
(FS) 4.74 0.89 1.40 3.70 2.39 
HY On10 _O217 4.74 1.34 0.37 1.63 ac 70 2.39 
MT 0.69 3.41 5.92 1.16 1.41 0.77 0.91 
IL 0.28 0.43 0.65 0.84 1.08 0.21 0.28 
HM 2.18 6.33 5.63 
™ 0.17 
RU oO.4g 
AP 0.09 0.07 0.67 0.07 
PY 0.02 0.13 0.09 0.17 0.04 
CLASS- 
IFICA- R R-T R R R R R R 
TION 
MAP 


UNIT 7A 7A 7B 7B 78 7c 7c 7c 
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high SiO, group. These chemical differences could be 
explained if the rhyolite classification here did not 
represent a single continuous related series of lavas, but 
rather two groups which were derived from comendites and 
trachytes respectively. Some of the high silica group have 
either phenocryst green sodic pyroxene or groundmass sodic 
amphiboles suggesting that they are related to comendites. 
The divisions according to Na,0 and K,O were not 
coincident and no significant differences in level for any 
other major components were noted in either case. There are 
two possible explanations for the Na,0, K,0 variation. 
Either they are random, or there is open system behavior for 
these alkalis (which does not affect the other elements the 
Same way). In that the rhyolites have the lowest Na,0O 
contents for salic rocks, and also the lowest Na/(Na+kK) 
ratio, the variation is probably not random. However, the 
lack of correlation between Na,O level and other oxides 
precludes a single mineral fractionation origin or 
alteration trends, as was suggested for the comendites. 
Covariance and correlation coefficients were calculated 
for the eight rhyolites of unit 7. Silica varies inversely 
with Al,O,, FeO*, MnO and CaO. TiO, varies normally with 
Al,O;, FeO*, MnO, MgO, CaO and P,0,. Al.O; varies normally 
with FeO* and CaO. MnO varies directly with MgO and CaO but 
inversely with Na,O. MgO varies inversely with Na,0O and K,0 
but normally with P,0,. CaO varies inversely with Na,0O. K,0 


varies inversely with P,0O,. The very unusual aspect of this 
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analysis is the lack of positive correlation between alkalis 
and any other oxide. The negative correlations all involve 
$i0,, Na,O and K,.0 and may imply open system behavior 
(selective loss) for these components. The positive 
correlations for the minor components could be explained by 
Fe-Ti-oxide and clinopyroxene fractionation or by dilution 
of these components with feldspar and quartz. 
Chemical Variation for all Salic Lavas 

Representative chemical analyses for salic lavas from 
elsewhere in the Intermontane Belt are presented in table 
7-25. Sample R-T from the Rainbows and HPD from the Heart 
Peaks are basically similar to the metaluminous trachytes of 
Level Mountain. For the average metaluminous trachyte from 
the Rainbows, all values are within the range of observed 
variation for Level Mountain. For the Heart Peaks dacite, 
the only values both out of range and beyond one standard 
deviation are the minor elements MgO, CaO and P,0,. The 
rhyolite from Edziza is similar to the high silica rhyolites 
from Level Mountain for all but the trace element MgO, which 
is greater for the Edziza sample. The peralkaline trachytes 
and pantellerites from the Rainbows and the pantellerite 
from Edziza fall within the variation of the peralkaline 
trachyte group from Level Mountain. The Edziza comendite 
compares favorably with the high silica comendites from 
Level Mountain. 

A series of eight oxide variation plots are presented 


for all of the salic lavas from Level Mountain, figure 7-18. 
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Figure 7-18. 


variation plots (anhydrous weight basis) for salic 
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of the stratocone stage. No distinction is made by map 


only by classification. For discussion of trends, 
to the text. The FeO versus Al,0,;, plot has been used 
by MacDonald (1974) for the classification of peralkaline 
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The four lava groups are easily distinguished on the basis 
of SiO., TiO,, Al.,O; and FeO*. are seen between peralkaline 
and metaluminous trachytes and between comendites and 
Continuous trends peralkaline trachytes. The rhyolites plot 
between comendites and trachytes, generally separating into 
two distinct clusters. The plot of FeO* versus Al,.0O,; is the 
basis for MacDonald's classification (1976) of peralkaline 
lavas. His division boundaries and labels for comendite C 
pantellerite P and their related trachyte types CT and PT 
are also shown. 

Fractionation models are presented for a typical 
comendite, peralkaline trachyte and metaluminous trachyte in 
figure 7-19. Mineral vectors correspond in direction to the 
removal of analysed phenocryst compositions and in length to 
the amount removed. The metaluminous trachyte used is flow 
8/26-56/6336 with its own phenocryst analyses for 
clinopyroxene (5%) Di,,Hd;;, olivine (5%) Fo., and iron 
titanium oxide (1%) Ulv,,;. Plagioclase (10%) An,,. and alkali 
feldspar (10%) Or... come from 8/30-82/6150, another unit 7a 
trachyte flow. Quartz is shown as a 5% fractionation vector 
(rather than assimilation) for plotting ease, and for 
consistency so that all vectors represent the shift of 
composition for fractionation. Quartz does not occur as 
phenocryst phase in any metaluminous trachyte. Assimilation 
of quartz would be parallel to this direction plotted but on 
the opposite side of the whole rock point. The main effects 


of any crystal fractionation and quartz feldspar 
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Figure 7-19, 
Mineral fractionation vectors for selected Salic 
compositions (anhydrous oxide basis). Symbols for lava types 
are aS in figure 7-18. 
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assimilation are in the direction of increasing SiO, and 
decreasing FeO*. The FeO - Al,0O;, plot is the only one 
showing a reasonable lineation for metaluminous trachytes. 
This variation is parallel to the clinopyroxene and feldspar 
vectors, which matches the dominant phenocryst phases. The 
trend on FeO*-Ti0, could be fit with a combination of 
feldspar, clinopyroxene and ulvospinel 73 as modelled, or by 
a more iron rich Fe-Ti oxide phase. The other variation 
diagrams show only a cluster of points for metaluminous 
trachytes. This cluster on most variation plots could be 
Spanned in two ways, either by more than 30% combined 
crystal fractionation of the phases indicated or by up to 
20% crustal contamination. A combination of the 
fractionation and contamination is also possible. For all of 
the variation diagrams not involving alkalis, the 
metaluminous trachytes are removed from the peralkaline 
trachytes in a direction which is inconsistent with minor 
amounts of crystal fractionation or crustal contamination. 
As little as 10% fractionation of any feldspar moves this 
trachyte from slightly metaluminous to slightly peralklaine 
and a stronger effect would be expected for the wollastonite 
normative members of this group. This classic "feldspar 
effect" has been alluded to in differentiation arguments for 
the origin of peralkaline rocks (Bailey and Schairer, 1964, 
1966). The problem here is that feldspar fractionation alone 
can neither explain the variation in the metaluminous 


trachytes nor generate the composition of any known 
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peralkaline trachyte (or other peralkaline lava) for Level 
Mountain. 

The peralkaline trachyte fractionation models used the 
average of the thirteen Level Mountain compositions in this 
group. Fractionated phases used were: alkali feldspar 25/5d 
(10%) Or;., quartz (5%), aenigmatite (5%), 
riebeckite-arfvedsonite alkali amphibole (5%), olivine Fa,; 
(4%), and aegerine Ac,,Hd,(5%). All of the ferromagnesian 
minerals were shown to demonstrate how little difference in 
direction and magnitude the variation in the ferromagnesian 
fractionating phase can make. The trends on the variation 
diagrams are consistent with alkali feldspar (anorthoclase) 
plus alkali pyroxene or aenigmatite fractionation. This is a 
good match to petrography. The lineation on S$i0,-Al,0,; is 
particularly inconsistent with any significant crustal 
contamination. Any quartz plus feldspar assimilation would 
be across trend. 

The only reasonable ways to reduce peralkalinity appear 
to be contamination by aluminous rocks or open system alkali 
loss. Two possible chemical statements of such a process are 
given below. 

19°tNa, S15 0,.+2HC1le22 2510 .a 4 NaClo+)H50 

2) 2NaAlSi,O, + 2HCl = Al,Si0O, + 2NaC1+H,0 
In either of these cases a higher pH.,O in the adjacent wall 
rock would favor preservation of the peralkaline condition 
while higher pH,0O in the melt would result in sodium (brine) 


and possibly silica metasomatism at the expense of the 
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magma's peralkalinity. These relations also apply to 
crystallization at the surface where pH,0 and pTotal are 
low. 

Two successful fractionation tests were found for 
deriving average comendite from average peralkaline trachyte 
using analysed phase compositions from Level Mountain flows. 
The low £0, case gave 40% fractional cyrstallization as (4%) 
Fa,s, (1%) clinopyroxene Di,,Hd;,, (28%) alkali feldspar 
Or,, and (7%) aenigmatite. The high £0, case gave 44% 
fractional crystallization as (5%) olivine, (0.1%) 
clinopyroxene, (30%) alkali feldspar and (7%) alkali 
amphibole. Both of these models provide good fits to the 
five major oxides Si0O,, Al,0O;, FeO*, and Na,O and K,0O which 
describe 99%+ of the beginning and ending compositions. MgO 
has the worst fit on both models, probably due to the choice 
of an olivine composition which was too magnesium rich. Use 
of an olivine composition with iron content greater than 
Fa,; would presumably improve the fit. This lack of real 
solutions for models involving fractionation of iron 
titanium oxide minerals matches the petrography, in that the 
peralkaline flows are generally devoid of opaques. In both 
modes the generated, P.O; and CaO are slightly high for both 
acceptable models, but not outside the range of compositions 
for such peralkaline lavas at Level Mountain. The small 
discrepancy for these oxides could be satisfied by less than 


0.2% of apatite fractionation. 
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An interesting feature of the fractionation models is 
the very small change in agpaitic index which accompained 
the 40% plus crystallization. The starting peralkaline 
trachyte had A.I. = 1.11 and the average comendite had A.I1. 
= 1.08. The first (low £0,) model predicted A.I. = 1.11 and 
the second A.I. = 1.12. These small changes in peralkalinity 
index match the overall distribution of Level Mountain's 
peralkaline compositions. Larger increases in A.I. would 
occur for feldspar fractionation alone and this in fact is 
observed for many of the peralkaline suites from elsewhere 
in the world for example at Pantelleria, Fant'ale, etc. At 
Level Mountain it would appear that the peralkalinity, 
(activity of ns in the melt), is effectively buffered by the 
crystallization of suitable alkali ferromagnesian minerals 
particularly sodic pyroxene, aenigmatite and alkali 
amphibole. 

For the comendite fractionation vectors rock 9/2-99 was 
reasonably close to the average so that it was used as model 
composition. The clinopyroxene Ac,,Hd,.(5%) and alkali 
feldspar Or;,.(10%) compositions belong to the same rock 
sample while ilmenite (0.25%) was from 8/27-62/6890 another 
glassy comendite. Quartz was 5% in length. The reason for 
the small ilmenite vector is the low titanium level in the 
rock which limits the amount of fractionation. It is 
interesting to note that the "feldspar effect" here is very 
subdued as 10% alkali feldspar fractionation produces no 


change in peralkalinity for the rock. The overall scatter of 
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the comendite compositions implies a combined fractionation 
of alkali feldspar, clinopyroxene and possibly quartz. 
Neither Fe-Ti oxide fractionation nor assimilation of quartz 
plus feldspar is consistent with observed variation trends. 

Attempts to derive rhyolites by fractional 
crystallization of metaluminous trachytes failed to produce 
the distinctive high K/Na ratios and the metaluminous to 
peraluminous character. Using the average rhyolite, and an 
alkali feldspar 8/10-43/5990 Or,, from a rhyolite 
ignimbrite, 10% feldspar fractionation reduced the A.I. from 
0.90 to 0.89. Feldspar fractionation in this case is 
obviously not a very strong lever for further reducing 
alkalinity. 

Three special projections were constructed, in the 
system Si0,-Al,0,-Na,0-K,0 (molecular) as a further aid to 
understanding the petrogenesis. All of the salic lavas lie 
more than 90% within this tetrahedral composition space, 
with the major non tetrahedral component being FeO*. The 
first of this series of plots, figure 7-20, is oriented 
perpendicular to the quartz-albite-orthoclase plane 
conventionally shown as 'Petrogeny's Residua System'. This 
plot was originally introduced by Bailey and MacDonald 
(1969) and has been used extensively to understand the 
genesis of peralkaline rocks. The compositions of the Level 
Mountain salic lavas diverge from the quartz-feldspar plane 
in such a way that the most peralkaline (most silica-rich), 


and most metaluminous (least silica-rich) compositions are 
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Figure 7-20. 
Level Mountain salic lavas in special projection of 
Na,O0-K,0-Al,0,-Si0, tetrahedron (molecular basis), after 
Bailey and MacDonald (1969). Horizontal axis corresponds to 
peralkalinity index, vertical axis corresponds to increasing 
Silica saturation. Midline is trace of petrogeny's residua 
plane (Q-Ab-Or). The lines correspond to the quartz-feldspar 
cotectics. 
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the furthest removed from the Q-Ab-Or plane. The Level 
Mountain lavas describe a continuous series parallel to the 
plane defined by the feldspar join and the chord 
(Na,0+K,0)=5, SiO0,=95 on the alkali silica face. Notice that 
the lavas are parallel to this plane on the silica-rich 
Side. This means that feldspar does not lie within the plane 
of lava variation, so that any alkali feldspar fracitonation 
will not lie in the plane of the variation. Unlike many 
peralkaline suites which have been studied, the rock 
compositions here do not stem from the feldspar join. 
Traditionally feldspar fractionation is invoked as a 
dominant differentiation mechanism and means of enhancing 
peralkalinity. In this case, even if the variation trends 
were Originally linear, alkali feldspar fractionation would 
introduce a curvature convex towards the silica apex (a 
positive radius of curvature about S$i0,). Such a curvature 
Or dispersion is not seen. The location of the individual 
qQuartz-feldspar cotectics is shown as they extend from the 
qQuartz-feldspar plane into the peralkaline volume, after 
Bailey and MacDonald (1969). The distribution of lava 
compositions probably defines the two feldspar minimum as it 
extends toward the Q-Ab-Or eutectic. 

The next figure 7-21 shows the projections of the 
peralkaline lavas and some of their accompanying feldspar 
phenocrysts. Again, this alkali ratio trapezoidal plot was 
introduced by Bailey and MacDonald (1969). Here the figure 


has been constructed for the plane parallel to the Level 
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Figure 7-21. 
Level Mountain peralkaline salic lavas projected onto their 
limiting plane in the peralkaline volume of the 
Na,O-K,0-Al,0;-Si0, tetrahedron (molecular basis) after 
Bailey and MacDonald (1969). Light lines join whole rocks to 
their alkali feldspar phenocrysts. 
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Mountain compositions. The purpose of these projections is 
to remove the distortion inherent in projecting compositions 
into the quartz-feldspar plane. The base of this trapezoid 
is the Ab-Or join at 12.5 mol percent Al.,0;, but variable 
alkali ratio. The top of the trapezoid is in this case the 
5% (Na,O+K,0), 95% SiO, chord on the alkali-silica face and 
thus has a value of 0% Al.0O;. Each peralkaline suite really 
requires the construction of its own diagram. The important 
use of this figure is not to compare with other rock suites, 
but to compare a given suite with the quartz-feldspar 
water-Saturated liquidus surface. Here the topology of the 
liquidus surface and feldspar minimum are shown as a 
function of pH,O. On the 5kbar surface the quartz-feldspar 
minimum corresponds to a ternary eutectic, data from 
Carmichael et al (1974) and Tuttle and Bowen (1958). The 
alkali ratios, (Na 100)/(Nat+K), for Level Mountain 
peralkaline lavas generally lie between 60 and 70. With the 
exception of one comendite, all of the salic lavas contain 
feldspars (anorthoclase) which are more sodic yet. The 
majority of the lava compositions plot to the Or side of the 
Quartz-feldspar minimum at any presSure. A Similar effect is 
seen for granites in the Q-Ab-Or traingle. This shift 
towards the Or composition presumably relates to projection 
away from other components such as FeO, MgO and CaO, or to 
an actual shift of the minimum toward Or in the more complex 
Natural systems (Naney and Swanson, 1980). The similarity in 


position and orientation of the Level Mountain silica 
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Saturated peralkaline rocks to the high pressure liquidus 
topology and the feldspar thermal minimum in Q-Ab-Or is 
great enough to propose that these may be a series of 
equilibrium eutectic-like liquids derived by a common 
fractional melting process, rather than fractional 
crystallization of feldspars. 

Figure 7-22 is an analog of Bailey and MacDonald's 
(1969) alkali ratio trapezoid plot as a projection into the 
aluminous part of the system. The plane shown here is 
defined by the point 10%A1.,0, - 90% SiO, on the silica 
alumina join with the alkali feldspar join as the base of 
the triangle. Note that metaluminous trachytes fall between 
60 and 70 for alkali ratio (100 Na)/(Na+K), as did the 
peralkaline trachytes across the Q-Ab-Or plane. The 
rhyolites show much more scatter, with reduced alkali level 
and reduced Na/(NatK) values. The alkali feldspar 
compositions are all more sodic than their whole rocks, and 
just as in the peralkaline part of the system, here too, 
they represent a strong lever. However, as before, the 
alkali feldspar does not appear to cause the variation trend 
for the majority of the lavas. 

The plot most commonly used to assess the variation of 
alkalis and alumina is (Na+K)/Al versus Na/(Na+K) (atomic), 
which is analogous to a Slice through the salic tetrahedron 
at a constant SiO, level. While it would be preferable to 
construct this plot at a specific SiO, level, it would be 


difficult to choose a particular $i0, value because of the 
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Figure 7-22, 
Distribution of metaluminous trachytes and rhyolites with 
alkali ratio projected onto their limiting plane in the 


metaluminous volume of Na,0O-K,;0-Al,0,;-Si0, (molecular 
basis). 
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variation (8% mol or 13% weight) in, the salic lavas. For 
this reason the A.I. versus alkali ratio plot is retained. 
Plotted on this figure 7-23, are all of the salic lavas and 
corresponding data for alkali feldspar phenocrysts. The 
feldspars and whole rocks do not both lie in this plane. 
This figure emphasizes the leverage that the highly sodic 
anorthoclase phenocrysts have in any fractionation scheme. 
Feldspar vectors are nearly orthogonal to the trends for 
peralkaline trachytes and rhyolites, and in the case of the 
comendites they appear to generate dispersion rather than 
any clear trend. The only case of a feldspar vector being 
parallel to a rock composition trend is for the metaluminous 
trachytes. The straight line is the least squares fit to the 
1kbar water saturated liquidus minimum given by: A.I. = 9.48 
Na/(Nat+K) - 4.97; r?=.79 using data from: 1) Tuttle and 
Bowen (1958), 2) Thompson and MacKenzie (1967) and 3) 
Carmichael and MacKenzie (1963). The compositions from 
Carmichael and MacKenzie are (4.5%Ac + 4.5% Ns) at A.I. = 
1.46 and for (8.3%Ac + 8.3% Ns) at A.I. = 1.98. The effect 
of indicated feldspar fractionation on metaluminous 
trachytes is toward the liquidus minimum plane. Apparently 
the peralkaline trachytes and some of the comendites 
Straddle this plane for the liquidus minimum, which was the 
Same deduction from the previous plots. Their dispersion 
along this direction is again not obviously related to 
anorthoclase fractionation, although such a path is possible 


from peralkaline trachyte to comendite. This path length was 
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Figure 7-23, 
Plot of alkali ratio versus agpaitic index for Level 
Mountain salic lavas. Liquidus minimum and alkali feldspar 
fractionation trends after Bailey and MacDonald (1974). The 
Level Mountain lavas appear to originate near the liquidus 
minimum. 
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equivalent to 40% to 44% net fractionation, of which 30% was 
anorthoclase. The Level Mountain lavas do not lie along a 
Single smooth curve for extensive feldspar fractionation as 
do most of the peralkaline suites reported in the 
literature. Typically, when alkali feldspar fractionation is 
the dominant genetic process, lava compositions should lie 
along smooth curves such as the examples for Pantelleria 
(Brotzu et al, 1974) or Fant'ale (Gibson, 1974). Note that 
the latter portions of these trends, at high peralkalinity 
index, coincide with the trace of the liquidus minimum. The 
Level Mountain lavas apparently attained this liquidus 
minimum position at a lower peralkalinity and sodium level 
than 1S typical for peralkaline suites. The Level Mountain 
lavas do not lie along such a singular feldspar 
fractionation trend. Level Mountain's peralkaline trachytes 
may represent a primary melt in the peralkaline part of the 
Salic tetrahedron. Their distribution seems to match the 
liquidus topology between one and five kilobars pH,O. 
Comendites may be derived from the peralkaline trachytes by 
fractional crystallization involving significant amounts of 
at least two mafic minerals in addition to anorthoclase. 
The rhyolites lie along a trend which extrapolates to 
O=Na/(NatK) at O=(Nat+K)/Al, which would represent either 
contamination by aluminous rock or open system alkali loss. 
FeO* is their only significant non-tetrahedral component. 
Its effect as seen in the distribution of natural granites 


is to shift the composition of the thermal minimum towards 
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Or. This is the direction seen for the Level Mountain 
rhyolite trend. The strong reduction in peralkalinity seen 
here 1S not common and is not explained by the iron 
component (3 to 4% mol FeO* should not have such pronounced 
effect). The rhyolites are also not easily related to the 
metaluminous trachytes by feldspar fractionation and they do 
not follow the variation trends for the metaluminous 


trachytes. 


TRACE ELEMENT CHEMISTRY 


Introduction 

Trace element concentrations for a related suite of 
igneous rocks can provide added insight into the controlling 
petrogenetic processes. The first set of models attempting 
to quantify the behavior of trace elements during a magmatic 
process was presented by Neuman, Mead and Vitaliano (1954). 
This work treated the process of fractional crystallization 
uSing the Nernst distribution law. The current state of 
quantitative trace element models for magmatic processes is 
reviewed by Allegre and Minster (1978). Arth's (1976) review 
is mathematically similar and presented in a more concise 
and easily understood fashion. He also presents 
representative and average partition coefficients for the 
major rock forming minerals for the range of volcanic 
compositions from basalt to rhyolite. More recent data of 


Partition coefficients is given by Pearce and Norry(1979). 
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Interpretation of trace element patterns for suites of 
igneous rocks is still an "inverse problem" as discussed by 
Minster et al (1977) and as such does not have a single 
correct solution. 

Aside from modelling, a great deal can still be done in 
a qualitative sense with trace element data. Trace element 
Signatures provide useful comparative tests between suites 
with similar genetic processes. Estimates of source 
heterogeneity can be made in addition to using trace 
elements to distinguish comagmatic suites from unrelated 
ones (Frey, 1980). One can also compare primitive lavas 
between different tectonic areas to get at the fundamental 
differences in the upper mantle, for example White et al 
(1975, 1979) in their comparison of the Azores Plateau to 
the mid Atlantic Ridge. Comparisons of trace elements 
between geochemically similar magmatic suites can be used to 
assess the state of mantle fertility or of relative 
depletion from place to place (Nixon et al, 1981; Hervig et 
al, 1980). 
Data and Analytical Methods 

Trace element analyses for selected Level Mountain 
lavas are presented in table 7-26. Most of the analyses were 
performed by XRF (J.G. Holland, analyst) on whole rock 
powders using previously calibrated rock standards after the 
method of Brown et al (1977). This includes most of the 
values reported for Cr, Ni, Cu, Zn, Rb, Sr, Y, Zr, Nb and 


Ba. Additional elements and some repeats in the above set 
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Table 7-26. Trace Elements. 


Element-Atomic 01 cB h PAP 
No. 
3 
4 
24 132 378 288 319 
28 77 260 231 232 
29 43 48 66 69 
30 108 105 98 89 
37 21 10 15 14 
38 576 507 540 561 
39 31 19 18 16 
40 221 116 #130 156 
41 38 22 25 29 
42 
56 711 #330 347 344 
79 
82 4 9 
930 3 1 
92 ce) 


Element-Atomic 8/21- 29/1f 25/3b 


No. 6/3222 late 
3 26 
4 2 
24 60 87 82 
28 67 48 49 
29 44 78 43 
30 188 132 93 
37 324 30 13 
38 809 633 
39 30 28 
40 259 198 
41 56 31 
42 6 
56 651 319 
79 
82 25 
90 
92 


29/11 
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13-K KD-1 16/19 


° 6 
2 3 
363 246 38 
195 251 29 
96 56 36 
128 98 133 
9 19 22 
510 688 862 
21 20 32 
131 145 287 
24 30 49 
2390 433 398 
a 5 
5 
12) 

8/28- 8/25- 
58/5895 50/6397 

Ss S 

2 2 

74 188 

41 112 

41 153 

133 119 

12 

382 464 

30 

185 

32 

6 
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19 6 

3 
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2 
218 129 
138 129 
61 BT 
113 127 
3 24 
441 934 
20 20 
106 252 
12 47 
174 472 
4 4 
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Table 7-26. Continued 


Element-Atomic 24/1Hi LMI7 9/2- LMI1130b 8/25- 8/12- 8/26- 8/10- 8/5- 8/T- 24/2C 
No 96/5330 52/5880 33/6140 56/6336 1/4915 8/5150 18/5720 
Li 3 20 10 14 20 25 32 
Be 4 | 5 4 9 10 9 
Cr 24 7 14 8 10 1 5 5 4 3 4 10 
Ni 28 5 8 5 3 9 9 6 5 23 5 6 
Cu 23 5 1 71 8 19 34 7 31 28 5 20 
zn 30 119 240 233 178 107 161 163 281 501 340 207 
Rb 37 107 109 75 98 115 89 149 123 
Sr 38 14 4 14 64 14 29 26 3 Ss 6 31 
Y 39 46 60 96 54 45 54 90 B4 
Zr 40 895 965 740 784 637 620 1085 996 
Nb 41 108 121 108 90 70 89 155 122 
Mo 42 6 6 15 13 
Ba 56 323 176 330 966 723 1079 62 377 
Au 73 4 3 
Pb 82 23 22 32 16 293 
Th 30 13 
U 92 3 
Element-Atomic 8/6- 8/6- LMIII30a 25/5d 25/5e 8/16- 9/2- 8/16- 8/25- 9/2- 9/2- 
No. 12-5290 12n 43ccpp 95/5625 44/5960 54/6345 ss 99/6600 
Li 3 16 7 34 69 64 54 61 
Be 4 8 6 11 13 13 11 13 
Cr 24 5 3 13 15 7 5 15 3 10 3 3 
Ni 28 14 3 2 5 7 8 7 1 
Cu 29 14 46 vA Y 11 10 8 16 11 16 14 
Zn 30 323 436 188 356 216 225 315 261 
Rb 37 124 117 116 168 154 310 275 276 
Sr 38 4 33 2 19 16 ° 1 ie} <i <1 <1 
yv 39 63 155 73 128 90 116 128 ta 
zr 40 71212 1122 903 1338 1306 580 | 1592 1617 
Nb 41 153 165 127 176 133 133 166 164 
Mo 42 
Ba 56 16 29 1085 147 127 4 23 6 
Au 73 3 3 3 4 3 3 3 
Pb 82 
Th 90 
U 92 
Element-Atomic LMI3 8/27- 8/16- 8/16- 8/27- 8/27- 
No. 66/5862 43/5990 44/6350 63/6788 62/6890 
ee EN. i ae RS ST ae ae 
ui 3 35 7 8 30 43 
Be 4 11 6 6 8 9 
Cr 24 5 11 3 3 5 3 
Ni 28 4 9 10 9 26 
Cu 29 5 19 12 19 30 8 
Zn 30 211 223 262 221 258 
Rod 37 296 248 162 125 217 
Sr 38 1 5 22 1 ° 
y 39 62 113 WA | 67 109 
zr 40 771 566 834 766 1148 
Nb 41 126 132 116 105 138 
Mo 42 12 
Ba 56 16 3 99g 945 22 
Au 79 3 4 3 3 3 
Pb 82 41 
Th 90 25 
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were analysed by AA under the guidance of A. Stelmach. 
Standards were multi-element solutions. Two U.S.G.S. 
international geochemical powdered rock standards were 
analysed along with the Level Mountain samples to check 
laboratory methods and data reduction procedures. These were 
AGV-1 andesite and BCR-1 basalt referenced in Flanagan 
(1972). One or both of these samples had elements in the 
range of interest for the Level Mountain suite. Elements 
analysed by this method included: Li, Be, Cr, Cu, Sr and Au. 
For samples analysed isotopically, Rb and Sr were determined 
by isotope dilution. Where these values duplicated XRF or AA 
data, the isotope dilution measurements were reported due to 
their higher levels of accuracy and precision. Neutron 
activation analyses were performed on six selected Level 
Mountain samples, which spanned the compositional range of 
lava types, by G. Goles of the University of Oregon. Data 
include the rare earth elements (La, Ce, Sm, Eu, Tbh, Yb and 
Lu), the naturally radioactive elements (Hf, Th, U, K), and 
analyses for Fe, Sc, Co, Cs, Sb, Br. 
Comparison with Other Suites 

A comparison for certain trace element abundances and 
element ratios in basalts is presented in Table 7-27. This 
compilation is intended to place Level Mountain and the 
other B.C. lavas in a global context. The world examples 
include oceanic tholeiites (MORB) and alkaline basalts, both 
from continental and oceanic settings. The B.C. examples are 


demonstrably alkaline by their overall trace element 
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Signature, particularly for the diagnostic ratios Nb/Y 
(Pearce and Norry, 1979) and Nb/Zr (Langmuir et al, 1977). 
If one can consider these basalts at roughly the same 
Stage of evolution or derivation, then their trace element 
abundances and ratios should reflect general characteristics 
of their mantle source regions. Note that Level Mountain has 
relatively typical levels of K, Ni, Sr, but low values for 
Rb/K, Rb/Sr, Nb/Zr, Nb/Y, Th/Ba when compared to the other 
alkaline basalt examples (Langmuir et al, 1977; Heinrichs et 
al, 1980). This probably indicates a relatively depleted 
mantle source with respect to the residual and incompatible 
elements for Level Mountain. Similar low Rb levels are seen 
at Heart Peaks (Casey, 1980) and in the Rainbow Range 
(Bevier, 1978). The only trace elements whose abundances 
seem disproportionately high for the Stikine are Ba and Pb. 
Discussion of the Level Mountain Trace Element Variation 
Selected trace element plots are shown in figures 7-24 
and 7-25. The ratio Zr/Nb has been shown to be relatively 
constant within comagmatic suites (Weaver et al, 1972). For 
all of the East African volcanoes the entire lava suite 
basalt-benmoreite-peralkaline trachyte plots as a straight 
line of positive slope passing through the origin (Ferrara 
and Treuil, 1974). The slope varies from centre to centre. 
That this line passes through the origin is a prerequisite 
of residual element behavior. Bailey and MacDonald (1975) 
have shown evidence for F, Zr, Rb, Cl, Nb, Yb and Zn from 


Eburru of the Kenya Rift, that these are not truly residual 
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Figure 7-24, 
Trace element variation plots for Level Mountain (Zr versus 
Zn, Rb and Ba). Salic compositions show considerable 
dispersion compared to basalts. A distinctive source 
composition and/or open system derivation is implied for the 
Salic lavas. 
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Figure 7-25, 
Trace element variation plots for Level Mountain. Salic 
lavas and basalts generally plot in separate clusters or 
along distinct trends. Basalts show residual element 
behavior for many elements while salic lavas do not. Salic 
lavas do not lie along fractionation trends from basalts. A 
separate origin is implied. 
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elements but evolve in an open system. When pairs of these 
elements are plotted, either they do not pass through the 
Origin or they show an inflection such that the basaltic 
segment includes the origin but the salic does not. In the 
Zr versus Nb example from Level Mountain, the basic and 
Salic segments neither extrapolate to the origin, nor do 
they intersect in positive concentration space. The 
implication using the criteria of Weaver et al (1972) is 
that there are two different source regions for the Level 
Mountain salic and basic magma series. The interpretation 
according to Bailey and MacDonald (1975) is that the salic 
magmas may require open system, perhaps with a Zr, Rb, F, Cl 
rich gas phase. On the Rb versus Zr plot for Level Mountain 
the basalts make one cluster at low values and the trachytes 
string out along a trend at high values. The comendites and 
rhyolites are removed on the high Rb side of this trend. A 
line can be fit to the basalt and trachyte data that passes 
through the origin, given by: 

Bireee). 125 azricze 3490 

where Rb and Zr values are in ppm. This line has a much 
higher Rb/Zr ratio than that from any reported alkalic 
volcano. If the trachytes were indeed derived by a classical 
alkali feldspar fractionation process, this trend would have 
to be flatter because Rb has a finite, positive bulk 
distribution coefficient for alkali feldspar, while that for 
Zr is zero (Bailey and MacDonald, 1975). The calculated 


series of possible alkali feldspar fractionation curves 
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(Bailey and MacDonald, 1975) have lower slopes Rb/Zr than 
the Eburru trend, which also makes them lower than the Level 
Mountain trend. They argue that no crystal-liquid process 
involving alkali feldspar can explain the observed trend 
(either melting or crystallization). This interpretation 
applies even more strongly to the Level Mountain data. The 
rhyolites and comendites were not included in the fit due to 
their high off-trend Rb values. This compares well with 
their high K/Na ratios indicating that Rb and alkalis 
require an open system model. On the Zr versus Zn plot, the 
best fit line to all data has a positive Zn intercept of 100 
while for the comendites and rhyolites a parallel trend at 
higher Zn is seen (intercept 182 on Zn). This time the high 
zn off-trend scatter is due to peralkaline trachytes. The Rb 
versus Zn plot shows a Similar positive Zn intercept near 
100 again implying open system behavior where Zn is 
associated with some volatile rather than crystal or liquid 
phases. The Li and Be data also imply an open system with Be 
carried in a gas phase. To fit all of the data on the Y 
versus Nb plot either requires two lines or an intercept 
that misses the origin, implying open system behavior or 
distinct sources or both. Most of the trace element pairs 
such as Sr versus Y, Sr versus Cr, sr versus Rb, etc. Show 
two perpendicular trends parallel to the respective element 
axes. This implies both separate populations and independent 
controls on the trace element variation. The strongly 


independent nature of the Rb and Sr variation at Level 
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Mountain contrasts markedly for centres with feldspar 
fractionation or fractional melting as the dominant 
petrogenetic process, compare for instance with Aden (Cox et 
al, 1970). The Sr versus Ba plot for Level Mountain is 
unique among alkali volcanoes in having a plagioclase-type 
Ba versus Sr dependence for the basalts but totally 
independent behavior for the salics. There is no possible 
Single stage fractionation or fractional melting argument to 
explain the derivation of the salic lavas from the basalts 
and still satisfy the Ba and Sr data. The Ba and Sr data 
also deny crustal contamination of basalts as a mechanism 
for trachyte genesis because they fail to fit the hyperbolic 
mixing model of Langmuir et al (1977). 

Considered as a set, the K,Rb and Ba concentrations for 
Level Mountain lavas portray the distinction between basalts 
and salics. The basalts and hawaiites lie along a single 
line in K-Rb-Ba space with Rb and Ba showing typical 
residual element behavior expected for either a fractional 
melting or fractional crystallization processes The 
projections of this line into Rb-Ba space is given by 
Rb = 030598"Ba <1'65454,cr?c#¢0.92 
and into K-Ba space by 
KRo#ens 1.7078) Ba: -© 2340.97, «r2%=0 596 
(where all element values are in ppm). By contrast the salic 
lavas have higher Rb and K values which are also too high 
for the projection of the basic trend. For the salic lavas 


K, Rb and Ba are variable and no longer correlated in any 
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Simple fashion. The salic lava types, trachyte to rhyolite 
and peralkaline to metaluminous, are totally interspersed 
along any apparent trends. This overlap of fields precludes 
Simple differentiation trends even within the salic group. 
Considering the large distribution coefficient of 6.12 for 
Ba into alkali feldspar versus whole rock (Philpotts and 
Schnetzler, 1970), a pronounced trend of Ba versus K and Ba 
versus Rb with onset of alkali feldspar fractionation is 
inevitable. The absence of such a regular trend for the 
Level Mountain salic lavas precludes their derivation by 
alkali feldspar fractionation from a basaltic parent. The 
Scattered K-Rb-Ba data for the salic rocks suggests either 
open system behavior for these elements or derivation of the 
diverse salic lavas from source regions with variable 
K-Rb-Ba. This might be expected either for very small 
degrees of partial melting in an inhomogeneous metasomatized 
upper mantle or for fractional fusion of crustal rocks with 
variable alkali feldspar K-Rb-Ba compositions, see figure 
e20; 
Rare Earth Elements in Level Mountain Lavas 

Six lavas chosen to cover the compositional spectrum 
were analysed for rare earths, table 7-28. The chondrite 
normalized abundance patterns for these lavas (basalt, 
phonolite, tristanite, peralkaline trachyte, comendite and 
rhyolite) are given in figure 7-27. The enrichment for light 
rare earth elements (REE) is typical for alkaline lavas, but 


the low degree of light versus heavy enrichment is not. For 
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Figure 7-26. 
Plot of the (K, Rb, Ba) variation for Level Mountain lavas. 
The K, Rb and Ba abundances vary Sympathetically for the 
Level Mountain basalts but the salics do not follow this 
Simple residual element behavior. When K and Rb values are 
normalized to Ba, all Level Mountain lavas plot on a single 
trend which may be due to a fractional fusion process, or to 
some peculiarity of evolution of salic magmas in an open 
System. Alkali feldspar fractionation or assimilation may 
cause dispersion along the trend in either direction. 
Comparative trends are shown for peralkaline suites from the 
Rainbow Range (Bevier, 1978) and Fant'ale (Gibson, 1978). 
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Figure 7-27, 
Rare earth element patterns for selected Level Mountain 
lavas. Normalized abundance as (1.2xsample)/LEEDY CHONDRITE. 


(1) Pantellerite 8/5-8/5150 
(2) Rhyolite 8/27-63/6788 
(3) Comendite 9/2-99/6600 
(4) Tristanite 8/16-44/6250c 
(5) Phonolite LMI 20i 
(6) Basalt 13-k 
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comparison see trends given by Cox et al (1969), Ferrara and 
Treull (1974), and White et al (1979). Also the small 
difference in REE abundance levels for the range of lava 
types at Level Mountain is unusual. Normally the peralkaline 
lavas show REE abundances up by a factor of 6+ 3 from the 
basalts whereas here the patterns cross. Compared to the 
other types, comendite and rhyolite have markedly lower 
heavy rare earth levels than the other lavas. Heavy rare 
earth depletion in general is characteristic of a source 
depleted in clinopyroxene and garnet or of suites derived by 
Substantial fractionation of such phases. Here it probably 
reflects some preferential uptake of heavy rare earths by 
sodic clinopyroxene or other sodic ferromagneSian minerals. 

The single most important feature of the rare earth 
patterns for the salic rocks is the total absence of a 
negative europium anomaly. Peralkaline salic lavas typically 
show strong europium depletion arising from alkali feldspar 
fractionation, in fact Eu depletion is documented as a 
characteristic of peralkaline lavas, (Ferrara and Treuil, 
1974). Sometimes the concurrent or related intermediate 
lavas show positive europium anomalies. Here again this is 
not seen. 
Summary 

There are many distinctive features of the Level 
Mountain trace element signatures for salic lavas. While 
they show enrichment in alkalis, REE, Zn, Be, Zr, Nb, Th, Hf 


and depletion in V, Cr, Ni, the degree of enrichment 
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compared to the East African Rift, Azores and Pantelleria is 
subdued. In most suites U and Mo are enriched, while Sc, Cu 
and Eu are depleted. These aspects of the typical 
peralkaline trace element signature are not observed at 
Level Mountain. Three possible explanations for these 
differences are proposed: i) Level Mountain may have a 
partially depleted mantle source region, ii) There is the 
possibility for open system behavior for many of the 
residual elements and possibly a distinctive gas phase 
associated with the Level Mountain peralkaline salics and 
iii) Lack of pronounced alkali feldspar fractionation and 
development of strong peralkalinity has failed to generate 
the pronounced differences in trace element behavior seen 


for other centres. 


OXYGEN ISOTOPE STUDY FOR LEVEL MOUNTAIN LAVAS AND MINERAL 


SEPARATES 


Analytical Methods 

Oxygen was extracted by fluorination reaction with 
bromine pentafluoride as described by Clayton and Mayeda 
(1963). From silicate samples O, gas is evolved with a 
residue of BrF, and fluoride salts. Samples were air dried 
and weighed, generally a 20mg sample was loaded into the 
nickel reaction vessels. A two hour heating to 300°C in 
vacuum was performed to remove absorbed water and air. The 


BrF, extraction took a nominal twelve hours at 650°C. Oxygen 
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was purified over in a cold trap then combusted with carbon 
of known isotopic composition. The yield of 0, (CO,) gas was 
determined from the pressure differential on a mercury 
manometer. Gas samples were transferred to glass sample 
tubes and and their isotopic ratios were measured. The mass 
Spectrometer used for these determinations was a Micromass 
Model 602D located in the department of Soil Science. 
Isotope ratios were reduced on line with an HP programmable 
calculator. Continuous correction for instrumental drift was 
made by comparing each sample to preceeding and subsequent 
standards. The data are reported with respect to the SMOW 
Standard assuming alpha (CO,-H,0) = 1.0407. Values were only 
accepted if their oxygen yields were within 5% of 
theoretically predicted amounts. Minerals were checked 
against ideal stoichiometry while lava yields were 
calculated from ten oxide anhydrous chemical analyses. Most 
yields are within 3%. 
Level Mountain Oxygen Isotope Study 

Analysis of the oxygen isotope composition for Level 
Mountain whole rocks and mineral separates was begun in 
1976, with the assistance of Dr. Karlis Muehlenbachs. At 
that time, none of the Late Cenozoic volcanics in B.C. had 
been analysed for oxygen. The primary impetus for the study 
was to see if the upper mantle for B.C. and the volcanic 
rocks derived from it were isotopically normal or exhibited 
low '*O character as discovered for Icelandic basalts by 


Muehlenbachs et al (1974). It was also anticipated that the 
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oxygen data would place constraints on petrogenetic and 
evolutionary processes within the Level Mountain suite. In a 
study of thirty-five samples from four Japanese Islands, ~ 
oxygen isotope data were shown to monotonically increase 
with differentiation by fractional crystallization 
(Matsuhisa et al, 1973). In an another study on andesite 
from the Banda arc, IndoneSia , oxygen and strontium isotope 
systems were positively correlated (Magaritz et al, 1978). 
For this latter example a two component mixing model was 
constructed for typically low ‘'*O and low *’Sr/**Sr mantle 
derived magma and subducted sediments. Aside from detailed 
interpretations such as those mentioned above, the Level 
Mountain samples could also be compared to other alkaline 
lavas of known isotopic composition. Some of these data have 
been previously reported by Hamilton, Baadsgard and Scarfe 
(1978). 

A tabulation of the Level Mountain oxygen isotope data 
is given in table 7-29. Values are reported as 6'*O relative 
to the SMOW standard. The order of presentation is 
Stratigraphic with most of the eruptive history having been 
sampled from the base of the Plateau, map unit 1, to unit 8 
of the stratocone. Descriptive notes plus map location are 
given to the side. All of the whole rock values reported 
here are thought to be fresh and unaltered. For a majority 
of samples this is demonstrated by mineral separates. Whole 
rock basaltic values range from 5.7 to 6.5°/,,. for Level 


Mountain. 
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The Level Mountain basalt analyses and the chrome 
diopside from the Telegraph Creek lherzolite nodule indicate 
that the mantle under the Stikine is normal with respect to 
its oxygen isotope composition. The basalts of the plateau 
are all around 6.5+0.05°/,, while those of the stratocone 
are lighter at 5.8+0.05°/,,. This difference probably 
reflects two compositionally distinct sources, or 
alternatively fundamentally different processes such as the 
influence of gas phase or degree of partial melting. The 
different oxygen signature for plateau and stratocone basic 
lavas 1S not associated with "differentiation" because 
alkali basalts and hawiites are sampled in each case and are 
identical within each set. This difference would also not 
follow from fractional crystallization because in the 
basaltic compositions olivine is the dominant fractionating 
phase along with smaller amounts of magnetite, clinopyroxene 
and plagioclase. Starting from a typical plateau basalt of 
6.5°/.. such fractionation would drive the residual liquids 
to heavier oxygen compositions due to the lighter nature of 
most of the fractionating phases. Any fractionation which 
may have occurred for the plateau basalts seems to have a 
negligible effect in dispersing the oxygen signatures. The 
stratocone basalts are lower in Mg and higher in normative 
plagioclase. If they were plagioclase-enriched couterparts 
of the earlier plateau magma they would be isotopically 
heavier, not lighter as they turn out to be. The same 


increase would be expected for crustal contamination. 
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The salic lavas show a total range of compositions from 
5.6 to 10.2°/,,. This range encompasses the basalts and is 
quite large compared to any volcanic suite which has been 
explained by fractional crystallization. It also spans the 
low, intermediate and high '*O salic rock subgroups of 
Taylor (1968). By comparison other peralkaline glassy rocks 
fall in the middle of this range; see Fant'ale, Timor and 
Nevada values in table 7-29. The fractionation between glass 
and alkali feldspar inreases as the whole rock 34'*O 
increases. The largest fractionation (2.0°/,,.) is for an 
anhydrous comendite, 9/2-95/5625. The 6'*O variation in 
these glassy rocks is not related to the water content. Ten 
hydrated obsidians, welded tuffs and perlites with 
corresponding feldspar values from Taylor (1968) show 
fractionations ranging from 2°/,, to 9°/,., with whole rock 
values all greater than 10°/,.. The 6'*O rich comendites, 
trachytes and rhyolites from Level Mountain clearly do not 
fall into this severly altered league. The pitchstones in 
particular are dense and optically clear. An alternate 
explanation of the Level Mountain salic rock variation must 
be sought. 

The alkali ferromagnesian minerals (both riebeckite - 
arfvedsonite type amphiboles, aenigmatites and sodic 
pyroxenes) are isotopically lighter than whole rocks 
(glasses) and alkali feldspars. This is not just an artifact 
of the Level Mountain peralkaline pitchstones because the 


Same relationship is seen for the Fant'ale pantellerite. The 
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alkali amphiboles fall in the general range for plutonic 
amphiboles of more ordinary chemistry (without peralkaline 
affinity), as does the sodic pyroxene fall in the range of 
other igneous pyroxenes. The two aenigmatite analyses, while 
predictably lighter than their whole rock and alkali 
feldspar, must stand on their own as there are no previously 
published values. Since all of the alkali ferromagnesian 
minerals and alkali feldspars are lighter than their whole 
rocks, it could be argued that ordinary fractional 
crystallization would have the effect of driving up the 
residual 6'*O value. Even for large fractionations, this 
effect is sufficient to explain the great variation. 

From major element chemistry the derivation of the 
average comendite from the average peralkaline trachyte 
required 45% by weight fractional crystallization as 30% 
alkali feldspar, 7% olivine, 7% amphibole or aenigmatite and 
1% clinopyroxene. This satisfies bulk oxygen balance to + 
1.25°/,,. When these proportions are used with respective 
oxygen contents and '*O/'‘O ratios, the mass balance can 
also be performed for isotopes. The following model was 


calculated: 


PTrach Ofe.e fore Hd, Aas (AGNI Rieb Com. 


O wt.% 46.49 47.81 32.80 37.15 41.28 41.02 48.8 
8'*0°/oo 6.9 6.7 5 6.0 5.6 6.0 7.30 
fraction 1 =,'36 =e07 a ~307 =.07 £55 
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For this example with 45% fraction crystallization, even 
with all of the minerals lighter than residual glass, the 
shift is only 0.4°/,.. If instead of fractionating, the 
alkali feldspar and other minerals crystallize at isotopic 
equilibrium the effect would be truly negligible. Another 
example was calculated for the fractional crystallization of 
40% alkali feldspar from comendite 9/2-99, using tabulated 
values to give a shift from 8.00°/,, to 8.47°/... Clearly 
fractional crystallization is too limited to explain the 
natural variation in isotopic ratios for the Level Mountain 
Salic lavas. In the case of the basalts and hawaiites, which 
have minerals on both sides of the whole rock ratio, 
calculated fractionation effects are negligible. 
Traditionally the '*O values for volcanic suites have 
been compared to chemical differentiation indices. Matsuhisa 
et al (1973) used Kuno's solidification index, (SI = 100 
MgO/(MgO+FeO+Na,0+K,0)), for the Japanese Islands while 
Taylor (1968) used Larsen's index, (L.1I. = 1/3Si0, + K.0 - 
CaO - MgO - FeO*) for volcanics from all over the world. 
When these indices are calculated for the Level Mountain 
samples and plotted against 6'*O, the basalts and salics 
plot in different clusters that are not obviously related to 
mixing lines or normal differentiation, see figure 7-28. The 
oxygen whole rock values were also plotted against percent 


partial melting, which was calculated from chemical analyses 
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Figure 7-28. 
Variation plots for oxygen isotope whole rock composition of 
Level Mountain lavas: (1) versus Kuno's solidification index 
showing the trend for fractional crystallization in Oki Dogo 
lavs after Matsuhisa (1973), (2) versus Larsen's 
differentiation index showing field of all volcanic rocks 
after Taylor (1968), and (3) versus maximum percent partial 
melts for basic lavas with an amphibolite source. The trends 
for basic and salic lavas are distinct in all cases. 
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using a pyrolite source for basalts and an amphibolite 
source for the salics. As with the other two plots, any 
inferred trends for the basalts produce little or no shift, 
while the gross 6'*O variation for the salics is independent 
of the chemical variation index. 

In attempting to find some chemical variable that 
correlated with 6'*O for the peralkaline lavas, oxide 
variation diagrams and element ratio plots were constructed. 
A weak trend was obtained for sum of alkalis with a slight 
improvement when reciprocal sum of the alkalis was used. 
This plot is shown in figure 7-29. The best fit line to the 
Salic rock data is: 

tAtNa sO +5K30)5 = 43 7635410 6720 2 6Or 07048 
This is a weak trend, nearly horizontal, and the r? for all 
of the salic rocks is only 0.51. The basalts, as on the 
other variation diagrams, plot as a separate field. The 
Salic line is not obviously a mixing trend with any known 
reservoir in mantle, crust or meteoric water. Local meteoric 
waters from Whitehorse and Fort Smith are -23°/,. and 
-22°/,,. respectively. Meteoric water from Level Mountain 
should not be appreciably different since it is of similar 
latitude. For fresh water with less than 500ppm total 
dissolved solids the reciprocal sum of alkalis cannot be 
lower than 1300. Local mantle is probably 5 to 6°/,, with 
0.7% total alkalis at a maximum. While actual crustal 
compositions are unknown, they are probably higher in 3'*0 


and lower in alkalis than the peralkaline magmas, if they 
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Figure 7-29, 

Variation plot for oxygen isotope whole rock composition 
versus reciprocal sum alkalis. Salic lavas from the 
Stratocone stage define a line which may imply open system 
behavior. This line is not related to a crustal component, 
to derivation from basalt or to meteoric water. Basalts plot 


as a cluster showing slightly heavier oxygen for the plateau 
lavas than for those of the stratocone. 
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are similar to calc alkaline plutonics reported by Taylor 
(1968) for the western U.S.A. All of these points are well 
removed in their location (in 6'*O versus reciprocal sum 
alkali space) and direction from the salic lava trend for 
Level Mountain. Two pantellerites from the Rainbow Range 
have been isotopically analysed and they also fall on the 
trend. 

There is also a linear relationship (r?=0.51) for the 
Salic lavas and Na/K ratio given by: 

(Na/K) = -0.16683'*O + 2.95 

A plot of this relationship is given in figure 7-30. These 
relationships with oxygen and alkalis could possibly 
indicate an open system gaseous exchange process between 
peralkaline magmas and granodiorite crust. Assuming the 
crust to be calc-alkaline mesozoic plutonics with 6'*O of 
13°/,. and the two Level Mountain alkali lines to represent 
gaseous transfer type mixing, this crust could be predicted 
to have sum alkalis equal to 8.50% with an Na/K ratio of 
0.78 or weights of Na,O =4.49% and K,0 =4.01%. This 
prediction is reasonable for a typical granodiorite and the 
mesozoic plutonics which outcrop in the Intermontane Belt. 


The model would then be 


Na, 1 6 a 
Salic Melt [ (CO6- Cine) 1% Granodiorite Crust 
5.840.3°/eeo mI Re ote ] 13.0+70.50°/,, 
+ (H,O) 


This type of gas-coupled alkali exchange could explain the 
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Figure 7-30. 
Variation of oxygen isotope composition of salic rocks with 
Na/K ratio. Circled points are for pitchstones (comendites 
and rhyolites). Trend could possibly indicate open system 
exchange of alkalis and oxygen isotopes between a primitive 
comendite melt and crustal granodiorite. 
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open system behaviour for Na,0 and K,0 discussed in the 
major element variation of the rhyolites and comendites. For 
granodiorites to be that heavy in oxygen implies that they 
are second cycle or S type granites. If an oxygen survey of 
the B.C. calc-alkaline plutonics shows them to be lower in 
3'80, in the range of normal primary plutonic rocks (6°/o. 
10°/,0; Taylor, 1968), then the model fails. The major and 
trace element data also support this type of open system 
alkali behavior, independent of the oxygen data. 

The oxygen isotope data for Level Mountain's 
peralkaline lavas is consistent with a mantle derivation and 
subsequent modification by some open system higher level 
gaseous alkali exchange process. The oxygen data cannot 
eliminate either of two principal hypotheses. These are 
direct derivation of peralkaline melts from the mantle and 
derivation by protracted fractional crystallization from 
alkali basalts. However such derivation trends have been 
shown by Matsuhisa for Oki Dogo Island with a normal alkali 
basalt differentiation series (1973) and they do not follow 
the distribution of the Level Mountain data. The range of 
3'*O values for the Level Mountain salic lavas also suggests 


some other process must be involved. 
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STRONTIUM ISOTOPE STUDY FOR LEVEL MOUNTAIN LAVAS 


Introduction 

A strontium isotope study was conducted during 
1976-1977 with the purpose of determining the *’Sr/**Sr 
initial ratio for Level Mountain lavas. It was hoped that 
the strontium isotope ratios would verify a mantle 
derivation for the Level Mountain lavas and constrain the 
participation of crustal components in petrogenesis. 

In the course of the Level Mountain study the 
peralkaline lavas with low levels of Sr, typically below 
5ppm, and high Rb/Sr ratios were found to have a reasonably 
high radiogenic Sr component. The study was expanded to 
include accurate isotope dilution determinations of Rb and 
an attempt to generate Rb-Sr isochrons for these very young 
rocks. Due to the long half life of *’Rb (1.42x10 ~'' 
year~') and the low Rb/Sr ratios for common igneous and 
metamorphic rocks, the accurate application of this isotope 
clock is usually only applicable to rocks older than 100MY. 
The unusual peralkaline nature of the Level Mountain salic 
rocks seemed amenable to such an isochron attempt despite 
their young age (1.2 to 4.6MY from fission track dates on 
glass and paleomagnetic correlation) (see discussion in 
Subsequent sections). 

The Level Mountain Rb-Sr Data 
The compilation of reduced data from the Rb-Sr study is 


presented in table 7-30. Samples are presented in descending 
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Stratigraphic order with a brief description of rock or 
mineral type and any qualifying comments about the mass 
spectrometer run. The error quoted for *’Sr/**Sr is the 
Standard error of the mean for n blocks (10 cycles per 
block) of data. The value of n is reported in the adjacent 
column so that sigma (the standard deviation of that 
measurement) may be calculated and used in place of the 
standard error. The Rb values in parentheses are from XRF 
determinations by Dr. J.G. Holland. In comparing the XRF and 
isotope dilution measurements on different sample aliquots, 
the XRF determination was generally within 10%. For these 
cases the *’Rb/**Sr was calculated from the natural 
abundance of *’Rb, the XRF Rb value and the measured *°*Sr 
value all converted to atomic ratio. These examples were 
mostly basalts which have negligible radiogenic Sr gains 
since their time of eruption due to low Rb/Sr ratios. The 
basalts do not really figure in the isochron portion of the 
study, except aS a comparative estimate of initial ratio. 
The Level Mountain basalts have an average *’Sr/**Sr of 
0.70334+.00059 which compares within the standard deviation 
to Dr. N. Green's determinations (unpublished) from the 
Gorda Ridge of 0.7028+.0002 and Armstrong's unpublished data 
on the Cordilleran basalts. Mount Edziza has basalt values 
in the range of 0.7028+.0001 Green (unpublished). The low 
*7Sr/**Sr values indicate the Level Mountain basalts to be 
primitive and mantle-derived. While the Level Mountain 


average is a little high for the Intermontane Belt (others 
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are typically 0.7031) sample PAR is as low as N. Green's 
Edziza sample #2902a65 at 0.7026 and is in fact the most 
primitive value hitherto reported for the Intermontane Belt. 
The diversity in *’Sr/**Sr ratios for Level Mountain basalts 
aeief rom -0:..7025 .to 10 +7045. 

The intermediate lavas from Level Mountain, including 
the mugearite, tristanite and phonolite samples, show a 
Similar range in *’Sr/**Sr, but extend to higher values than 
for the basalts (0.70309 to 0.70506). These intermediate 
lavas also have primitive strontium values and are 
ultimately mantle-derived. Possible explanations for the 
higher Sr isotope ratios in phonolites would be: more 
radiogenic Sr in source, contamination by radiogenic crustal 
Sr, disequilibrium melting or a different degree of melting 
from basalts but same source, mixing between basaltic magma 
(primitive Sr) and a more radiogenic magma. 

The trachytes, peralkaline trachytes, rhyolites and 
comendites have lower Sr levels than the magmas previously 
discussed and Rb/Sr ratios (weight) that range from 5.0 to 
216.0. This lends itself to appreciable radiogenic Sr 
accumulation, even in the 4.5MY since the eruption of the 
oldest stratocone lavas at Level Mountain. The *’Sr/**Sr 
ratios range from 0.7041 to 0.7076 and even 0.7553 in the 
case of the one extremely radiogenic glassy rhyolite dyke. 
Assuming the initial *’Sr/**Sr ratio for this dyke to be 
0.7033 a “point” age may be calculated from the slope 


formula of 5.886MY. Even accounting for all possible errors 
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this age estimate should be good to +1.5MY. To accurately 
correct for this age effect and obtain true initial ratios 
for the Level Mountain salic lavas, isochrons were obviously 
needed. Before any correction for age there are rhyolites 
and trachytes with *’Sr/*‘Sr ratios in the range 0.704 to 
0.705. Such low ratios could be taken to indicate: mantle 
Origin, derivation from basalts, contamination of (low Sr) 
crustal melts with (high Sr) primitive basalts, partial 
melting of young cruStal rocks (only slightly radiogenic), 
etc. The low values by themselves are hardly conclusive. A 
combination of Sr with other information is needed. 

The effects of crustal contamination on the Sr isotope 
ratio and Sr concentration has been shown for andesite lavas 
of the Rio Grande Rift (Zimmerman and Kudo, 1979). On plots 
of Sr versus *’Sr/**Sr, the contamination trend shows higher 
isotopic ratios at lower Sr abundance, and at higher SiO, 
content. Such a plot has been constructed for the Level 
Mountain data, figure 7-31. Considerable variation exists in 
the basalt data that could be explained by two mixing 
models. Line (1) is a least squares fit to the high-ratio 
basalts given by: 
eee se "Gr at.) > =) 1.181x10"2 se (ppm) + 0.7103 
with a value of r? = 0.78. This trend could indicate 
assimilation of crust with 20ppm Sr and 0.710 ratio by 
primitive basaltic magma of 650 ppm Sr and 0.7025 ratio. For 
this linear model, sample PAR is primitive and 29/1L 


represents about 25% crustal contamination. Line 2 has ten 
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Kigure 7-37 . 
Variation plot of strontium versus strontium isotope ratio. 
Trachytes with the lowest Sr content are the easiest rocks 
to alter the isotope ratio either by decay of Rb, 
contamination with radiogenic crust or contamination from Sr 
rich basalts. The basalts may show 3 trends, one for crustal 
contaminatio and one for magma mixing with a primitive 
trachyte to yield phonolites. Analytical errors are plotted 
with bars except where they are less than the dimension of 
the point. 
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points ranging from basalts to phonolites and tristanite. 
The expression is given by: 

87Sr/86Sr(at.) = -2.27x10-*Sr(ppm) + 0.7041 

with a value of r?=0.60. The Sr ratio of the intercept is 
certainly not a normal crustal value. This model is thought 
to represent magma mixing primitive basalt with trachyte 
having 15ppm Sr and an initial ratio 0.704. According to 
this model, the phonolites scatter between 6% and 17% basalt 
contribution to a normal trachyte. Mugearite by this model 
represents about half and half basalt and trachyte. The 
"parent" basaltic magma may actually be hawaiite, or the 
hawalites may represent contamination of alkali basalt by 
trachyte. This is precisely the magma mixing origin 
attributed to the intermediate lavas on the basis of major 
element chemistry. The increase in SiO, along with *’Sr/**Sr 
is considered to be coincidental, since silica as well as 
the other major elements fit a trachyte component just as 
easily as a Sialic crust. The only way that this mixing line 
could involve crust would be for the crust to chemically 
resemble a trachyte and have a primitive isotope ratio. 
Partially melted granitic gneiss xenoliths were found in 
both hawaiite unit 8 from Meszah Peak and in the tristanite 
dyke, 8/16-44/6250. In both of these cases the isotope ratio 
of the inclusions were similar enough to host magmas to 
assume that the Sr had been partially equilibrated along 
with the melting. The alternative is that these granitic 


gneiss xenoliths with *’Sr/**Sr= 0.7038 to 0.7051, could be 
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the primitive crustal component indicated by mixing line 2. 
The original magma mixing hypothesis, however, is still 
favored. The trachytes, peralkaline trachytes, rhyolites and 
comendites scatter along the isotope ratio axis at low Sr 
levels. No particularly good fit can be made to all of these 
points. It is thought that they represent variable 
contributions of crustal contamination or disequilibrium 
partial melting together with the overprint of radiogenic Sr 
formed since magmatism began. 
Isochron Attempts and Age Estimates 

Tf ithe atomic ratyYos *’Sr/*’Sr and *’Rb/**Sr are known 
for contemporaneous rocks or minerals, slope ages 
(isochrons) can be obtained with the zero *’Rb/**Sr 
intercept giving the *’Sr/**Sr initial ratio. The slope 
formula is: 
mest o*Sr) /(A*tBb/* *Sr-t = t 
where t, the half life of *’Rb, and t, the age, are in 
years. The slope of the isochron in the simplest case is 
determined by two points such as minerals from the same rock 
Or one mineral and the whole rock, coeval rocks of differing 
Rb and Sr contents. The slope can also be the least Squares 
best fit to a series of coeval minerals or rocks. An 
isochron plot has been constructed for salic and 
intermediate samples from stratocone map units 7, 6a, and 
5a, figure 7-32. Various fits were attempted for units 7 and 
5. Apparent ages increase upwards in the section. All three 


map units seem to have different initial ratios. Slope ages 
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Figure 7-32. 
Rb Sr isochron plot for whole rock and mineral samples from 
Level Mountain stratocone. Two point isochrons give 
Scattered ages which are generally too old compared to 
fission track dates. Initial ratios are quite variable. 
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are generally in the range of 5.1+1.2MY which agree with 
most of the age estimates obtained from single rocks, 
assuming the average basalt initial ratio of *’Sr/**Sr= 
0.7034. The problem is that all of the ages are too old 
compared to age estimates from fission track dates, 
paleomagnetism and geomorphology. Two single flows with 
mineral separates from unit 5a give such divergent answers 
as 3.9 and 23.8MY. Lumping all of the samples from the 
Stratocone gives a scattered apparent isochron of slope 
equivalent to 9.4+1.0MY, figure 7-33. The interpretation of 
all these "too old" ages is problematic and is not even 
restricted to Level Mountain. Souther and Symons (1974) have 
published K-Ar dates radiocarbon and fission track ages 
which tie Edziza's magnetostratigraphy to a post 7MY age. 
Some of the Rb-Sr data on comendites and trachytes from 
Edziza have *’Sr/**Sr as high as 0.715 (Green, samples 
#SE2902a66 and SE381567a, unpublished) implying ages ranging 
from 2MY to 24MY. Pantellerite flows from the Rainbow Range 
have *’Sr/**Sr as high as 0.7293 (R-88, Bevier, 1978) 
implying an age of 8.55MY, which is about 1.5MY older than 
the oldest K-Ar dates for the Rainbow Range. These old Rb-Sr 
ages are not restricted to B.C. Similar peralkaline suites 
from Aden, South Arabia have been reported by Cox et al 
(1970) and Dickinson et al (1969) to give isochron ages 
between 20 and 30MY older than the actual eruptive age of 
the Aden volcanoes. These isochrons are constructed by 


combining samples from six volcanoes with no attention to 
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Figure 7-33, 
Rb-Sr isochron plot showing whole rock values for salic and 
intermediate lavas of the sStratocone. From fission track 
dating this sequence of lavas is younger than 4.5MY yet they 
give an apparent isochron age of 9.4+1.0MY. This could 


possibly indicate the age of a melting event in the mantle 
source region. 
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actual eruptive sequence or stratigraphy. As at Level 
Mountain, there is considerable variation of *’Sr/**Sr 
within each of the Arabian volcanic centres. In the Arabian 
instance, the ages are thought to reflect the previous 
volcanic episode of the region (a mantle melting event). In 
the Intermontane Belt there are no conveniently timed 
earlier volcanic episodes. For the Stikine (Level Mountain 
and Edziza) the preceeding volcanic episode was the Sloko 
formation at approximately 50MY. There are also no 
conveniently placed tectonic events, the last being the 
cessation of subduction off the west coast, 200 km away, at 
about 20MY (Grow and Atwater, 1970). Souther (1970) stressed 
the relation of post Miocene alkalic volcanism in the 
Intermontane Belt to be subsequent to the cessation of this 
Subduction activity. The only remaining explanation is that 
the anomalous Rb-Sr ages for Level Mountain reflect some 
thermal event(s) which preceed any surface manifestation of 
volcanic activity. The 9.4MY lineation possibly reflects the 
mantle melting event responsible for voluminous plateau 


volcanism. 


LEAD ISOTOPE GEOCHEMISTRY FOR LEVEL MOUNTAIN 


Introduction 
The Pb isotope systems have been used to compare rocks 
derived from different Pb reservoirs (Schilling, 1973; 


Tatsumoto, 1978) in disucssing variable depletion of the 
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mantle source regions for oceanic island (alkaline) versus 
ocean floor (tholeiitic) basalts. U and Th posess low decay 
rates and low diffusivities (Seitz, 1973) with respect to 
thermal and tectonic events. As a result, the three 
radiogenic Pb clocks keep time slowly and the lead isotope 
ratios are sensitive indicators of geochemical 
characteristics which persist for hundreds of millions of 
years or longer. Consequently, they have been extensively 
used in modelling eoyuerndee turnover and mixing times for 
the mantle (mesosphere-asthenosphere) (Dupre and Allegre, 
1980; Hoffman and Hart, 1978; Sun and Hanson, 1976; and 
Armstrong, 1968). 

Of the four stable lead isotopes, ?°*Pb is the least 
abundant (1.4%) and cosmogenic while 7°*Pb, ?°’Pb and ?°*Pb 
are radiogenic. The radiogenic leads are normalized to *°*Pb 
to provide a common base for comparison. The radiogenic 
leads are the ultimate stable daughter isotopes of uranium 
and thorium decay. The lead isotope information is commonly 
displayed on two atomic ratio plots: 
2°&Db/2°*Pb versus 7°*Pb/?°*Pb and 
2°7Dbh/2°4*Pbh versus 7°*Pb/?°*Pb. 

The distribition of points on the first plot reflects the 
Th/U ratio (initial atomic isotope ratio 7°?Th/?°*U = K) of 
the source region(s). The distribution of points on the 
second plot depends on the U/Pb ratio (initial atomic 
isotope ratio ?°*U/*°*Pb = u) in the source region(s) and 


the joint operation of the two lead clocks. With knowledge 
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of, or assumptions about u and K, single or multiple stage 
model Pb ages may be calculated. The Pb isotope 
characteristics of abyssal tholeiites and oceanic islands 
have been reviewed by Tatsumoto (1978). He notes that all of 
these primitive volcanics have lead isotope ratios which 
correspond to negative model ages (future) for single stage 
growth. This feature and the wide scatter in ?°*Pb/?°*Pb he 
attributes to progressive U and Th enrichment with respect 
to Pb in the mantle source regions for basalts during much 
of geologic time. 
Level Mountain Lead Isotope Data 

The lead isotope ratio analyses are presented along 
with pertinent trace element data in table 7-31. There is no 
Simple correlation of the lead isotope data with "degree" of 
differentiation", major and trace element abundances 
(including lead), *’Sr/**Sr or ‘*O SMOW °/,.. The values and 
overall range in the Level Mountain lead isotope ratios is 
comparable to abyssal tholeiites, oceanic islands, Juan be 
Fuca - Gorda Ridge, East Pacific Seamount Chains 
(Kodiak-Bowie-Emperor), (Tatsumoto, 1978; Church and 
Tatsumoto, 1978), and the Mid Atlantic Ridge (Dupre and 
Allegre, 1980; Sun (1973). That the Level Mountain values 
compare to those other mantle derived suites, while failing 
to follow other geochemical variation, lends itself to 
mantle heterogeniety explanations rather than to 


contamination or mixing models. 
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The lead isotope data for Level Mountain have been 
plotted in figure 7-34. The distribution of the Level 
Mountain whole rock leads is consistent with a single 
source. The compositions of two conformable modern 
Cordilleran ore leads are shown on the figure by T for the 
Yukon Treadwell deposit S for the B.C Sullivan mine, data 
from Russel (1972). Leeman et al (1979) analysed leads from 
the iron rich lavas from Craters of the Moon, (C.O.M.) 
Idaho, which are crustally contaminated, and from the lower 
crustal charnocitic xenoliths which they bear. The C.O.M. 
ferrolatite plots midway between the Level Mountain group 
and the Cordilleran ore leads, while the gneisses are 
extremely primitive plotting well down the concordant trend 
mreache wacinity of 20 6Dh/2°4ph = 14.5+1.5. The purpose of 
pointing out these examples is to indicate the direction of 
crustal contamination. For the Level Mountain whole rock 
lead trend to represent crustal contamination, the virgin 
lead would have to be even more radiogenic than indicated. 
The only reported leads more radiogenic than those for Level 
Mountain are from Ross (Antartica) and St. Helena 
(Tatsumoto, 1978). 

The bold lines through the Level Mountain data are the 
best fit trends for whole rock leads from the Juan de Fuca - 
Gorda Ridge System and the Eastern Pacific Seamount Chain 
after Tatsumoto (1978). The trends for MORB (Dupre and 
Allegre, 1980) and the Hawaiian Islands (Tatsumoto, 1978; 


Sun et al, 1975) are very similar. Two of the Level Mountain 
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Figure 7-34. 
isotope ratio plots for Level Mountain lavas. Lava 


values plot near primitive MORB trend of Tatsumoto (1978). 
Variations are not obviously related to differentiation or 


crustal contamination. 


Implication is fractional melting 


derivation from a heterogeneous but primitive mantle source. 
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samples fall well to the right of the trend on the 

Seer bs-- “PD, versus 7°*Pbh/?**pb" plot, indicating a~ lower 
value of kappa which is to say that there is a low Th/U 
ratio for the source of the hawaiite and tristanite (crustal 
contamination ?). Despite the fact that the Th abundances in 
the Level Mountain lavas are high, other incompatible 
elements such as U and Hf are disproportionately higher 
compared to MORB and abyssal basalts (Dupre and Allegre, 
1980). If the Level Mountain 2°’Pb/?°*Pb versus ?°*Pb/?° ‘Pb 
data are fit with a single line, the slope is more nearly 
horizontal than for the Oceanic volcanics. This implies that 
the mantle source for the Level Mountain Leads is very 
young, i.e. it has been depleted and reset relatively 
recently. 

The abundance data for the Pb, Th, and U from Level 
Mountain whole rocks is presented in figure 7-35, (data 
tabulated in preceeding section on trace elements). The 
plots for U-Pb, Th-Pb and U-Th indicate a single homogeneous 
source region for all of the Level Mountain lavas, which 
compares favorably with the interpretation of a Single 
source for the Pb isotopes. This places constraints on the 
composition of the source region for Level Mountain. 
Tatsumoto (1978) has analysed the U and Th contents of 
mineral separates of a four phase periodite nodule (spinel 
lherzolite) from the strongly alkaline cinder cone at Boss 
Mountain (Takomkane) B.C. The spinel and Cr-diopside are 


enriched in U and Th relative to olivine and orthopyroxene. 
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Figure 7-35. 
Trace element variation plots for (U, Th, Pb) of Level 
Mountain lavas. Here the U/Pb and Th/Pb are lower than for 
alkaline volcanic suites of oceanic islands and the east 
African rift. The U/Th ratios for lherzolite minerals are 
from Takomkane Mtn, B.C. (Tatsumoto, 1978). 
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The Th/U ratios for spinel and Cr-diopside are 0.4 and 0.8 
respectively. Note that the Level Mountain Th/U ratios range 
from about 1.0 to 5.0 but mainly fall along a line of slope 
2.95. This range of Th/U ratios is typical for seamounts and 
oceanic islands as well. None of the aforementioned lavas, 
including those from Level Mountain could have been 
generated by melting of such depleted material as the 
Takomkane lherzolite. The synopsis of the lead isotope and 
related trace element data for Level Mountain is that the 
lavas were derived from a single inhomogeneous source for 
lead isotopes, but relatively constant in U/Pb, Th/U and 
Th/Pb. This source is more radiogenic than that for most 
primitive oceanic basalts and more heterogeneous in 
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